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Abstract 


The purpose of this research is to investigate the use of commercial, off-the-shelf wireless net- 
working technology in planetary exploration applications involving rovers and sensor webs. 
The three objectives of this research project axe to: 1) simulate the radio frequency envi- 
ronment of proposed landing sites on Mars using actual topographic data, 2) analyze the 
performance of current wireless networking standards in the simulated radio frequency en- 
vironment, and 3) propose modifications to the standards for more efficient utilization. In 
this annual report, we present our results for the second year of research. During this year, 
the effort has focussed on the second objective of analyzing the performance of the IEEE 
802.11a and IEEE 802.11b wireless networking standards in the simulated radio frequency 
environment of Mars. The approach builds upon our previous results which deterministically 
modelled the RF environment at selected sites on Mars using high-resolution topographical 
data. These results provide critical information regarding antenna coverage patterns, max- 
imum link distances, effects of surface clutter, and multipath effects. Using these previous 
results, the physical layer of these wireless networking standards has now been simulated 
and analyzed in the Martian environment. We are looking to extending these results to 
the and medium access layer next. Our results give us critical information regarding the 
performance (data rates, packet error rates, link distances, etc.) of IEEE 802.11a/b wireless 
networks. This information enables a critical examination of how these wireless networks 
may be utilized in future Mars missions and how they may be possibly modified for more 
optimal usage. 
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1 SUMMARY 


1 Summary 

Effective utilization of commercial wireless networking technology such as that based on 
the IEEE 802.11a and IEEE 802.11b standards, could bring the vision of proximity wireless 
networks used in future exploration of Mars to reality. However, in its current form, this 
technology may fall short of expectations due to its power requirements, limited tr ansmiss ion 
range, significant attenuation in hilly terrain, and possible poor bit error rate performance 
due to long multipath delays and narsn environments. The three objectives of this three-year 
research program are to 1) simulate the Radio Frequency (RF) environment of a planetary 
surface using high-resolution, Digital Elevation Maps of Mars, 2) analyze the performance 
of current wireless networking standards in the simulated RF environment, and 3) propose 
modifications to the wireless networking standards for better utilization in the planetary 
environment. 

In year two, we have successfully accomplished the second research objective through the 
completion of the following tasks: 

• Using the Year 1 results regarding antenna coverage patterns, ma xim um link distances, 
effects of surface clutter, and multipath effects, we have simulated the physical layers 
of IEEE 802.11a and IEEE 802.11b wireless networking standards in the Martian en- 
vironment. The simulations were conducted in MATLAB using the Comm Access 
Technologies mWLAN toolbox. These results will be used as the basis to begin the 
study of the behavior of the medium access layer for these protocols for the planetary 
environment. 

• We have measured key parameters (RF signal strength, delay spread, data rates, and 
packet error rates) in an outdoor IEEE 802.11b wireless network. These measurements 
were conducted at sites near Tortugas Mountain and Dripping Springs due their sim- 
ilarity to the Mars surface (free of man-made objects, little vegetation, mostly flat 
with some terrain variation and rocks, etc.). These measurements were compared with 
expected performance based upon the DEMs for the sites. We judge the agreement 
to be very good, especially considering the wide variation in measurements that can 
occur by moving the transmitter or receiver antenna by a small distance. 

In this annual report, we present detailed results for the second objective of analyzing 
the performance of current wireless networking standards in the simulated RF environment. 
Our results give us critical information regarding the performance (data rates, packet error 
rates, link distances, etc.) of IEEE 802.11a/b wireless networks. This information enables 
careful examination of how these wireless networks may be utilized in future Mars missions 
and how they may be possibly modified for optimal usage. Our work suggests that: 

• The results of the physical layer simulations for the Martian surface show that success- 
ful communications are possible within a few hundred meters of the transmit antenna 
when the transmit power is more than a few milliwatts and the antenna heights are 
fixed at more than 1 meter above the ground. The packet error rate performance of 
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802.11b without a RAKE receiver seems to be more adversely affected by the multi- 
path conditions than 802.11a. Further, the lowest data rate mode of 802.11a provides 
the best bit error performance. 

• Using higher power in the communications system does not always help the perfor- 
mance of the system. This result is known for terrestrial environments where vegeta- 
tion and atmospheric phenomena are important. It is also true in sparse environments 
with negligible atmospheric attenuation as well. 

• Transmission power and antenna height can be traded to a certain extent. The desired 
link Quality of Service and data rate may be more of a driver in link design and 
performance than transmission power and antenna height. 

• Use of a RAKE-type of receiver can significantly improve performance with 802.11b 
protocols. 

In addition, we include results from Year 1 regarding the validation of the power delay 
profiles. These results were not published in the 2004 annual report. 
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2 Introduction 

NASA’s long-term goals for the exploration of Mars include the use of rovers and sensors 
which intercommunicate through proximity wireless networks. Elements of the network have 
a short transmission range, low power requirements, low cost, and a relatively short-life span 
[1]. The performance of any such wireless network depends fundamentally on the Radio 
Frequency (RF) environment. In order to evaluate and optimize the performance of such 
a wireless network, a basic understanding or model of the channel is important. With 
such a model, better choices for the modulation and coding schemes, equalizer design, and 
positioning of network access point antennas can be made [2]. These choices can also affect 
the overall design and operations of a rover or a sensor web for a planetary environment. 

In the first year of the research, we simulated the RF environment of selected sites on Mars 
using commercial RF planing and propagation software (with appropriate modifications for 
the Mars environment) and high-resolution Digital Elevation Maps (DEMs) from the Mars 
Global Surveyor (MGS) mission. The results of the simulation provided critical information 
regarding antenna coverage patterns, maximum link distances, effects of surface clutter, and 
multipath effects. Using these results, we now have simulated and analyzed the physical 
layers of the IEEE 802.11a and IEEE 802.11b wireless networking standards in the Martian 
environment. Our results give us new information regarding the performance (data rates, 
packet error rates, link distances, etc.) of the wireless networks. This information enables 
careful examination of how these wireless networks may be utilized for proximity networks in 
future Mars missions and how they may be possibly modified for optimal usage. In addition, 
these results could assist mission planners in network design should such wireless networking 
technology be adopted. These studies form the basis for our examination of the medium 
access layers for these protocols that we have now begun and how the protocols may be 
expected to act in the outdoor planetary environment. 

The major accomplishments of the first year’s work are detailed in [3] and they include: 

• Identified a RF planning and propagation software packages (ATDI’s HerTZ Mapper 
and ICS Telecom) that are used in the cellular telephone industry for network planning 
and optimization. The software, with modifications developed at NMSU, allows us to 
import DEMs of Mars and tune key parameters of the propagation model for accurate 
simulations of RF coverage patterns and multipath. 

• Obtained high-resolution (10.4m/pixel) DEMs for the Gusev Crater and Meridiani 
Planum regions and successfully reformatted this data for import into the software. 
Within these regions, sites were chosen for their range and topographic features which 
'were expected to yield different RF propagation effects in simulation. 

• Computed RF coverage patterns and power delay profiles (multipath) for various com- 
munication scenarios at the selected sites. These scenarios include various transmitter 
and receiver locations, transmission power, surface clutter, and atmospheric phenom- 
ena. The basic reference configuration was a 2.4GHz carrier frequency, omnidirectional 
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antennas at a lm height above local ground, and 1W radiated power. The clutter layer 
was synthesized based on the statistical distribution of rock sizes and spacing found in 
earlier Mars missions. 

• Validated the simulation method by modelling a local site and performing actual RF 
field measurements. 

• Utilized the simulation results to answer questions that mission planners might wish 
to ask such as percentage RF coverage in a region or percent RF coverage as a function 
of distance from a network access point. 

The major recommendations from the first year’s work are 

• Investigate how the coverage results and power delay profiles may affect communica- 
tions protocol selection and/ or modifications. These investigations are expected to give 
further insights into the design of proximity network configurations and potential ways 
to modify the protocols for better performance. 

• Utilizing the RF models developed in this research, evaluate and improve IEEE 802.11a 
and 802.11b protocols for use in proximity wireless networks. 

• Investigate how NASA mission planners might partner with New Mexico State Univer- 
sity to make the RF modelling techniques developed in this research more widely known 
and utilized in the planetary exploration community. The methods developed in this 
research could be used to model RF propagation at other sites of interest on the moon, 
Mars, and even in structures such as the International Space Station. 

In this research regarding wireless network simulation on the Martian surface, we assume 
a proximity network model composed of multiple surface rovers and sensors that communi- 
cate via IEEE 802.11a or IEEE 802.1b to a network access point placed on the lander. At the 
lander, data is filtered, aggregated, and packaged for transmission back to earth via an or- 
biting relay station [4], Other assumptions in the network model are 1W transmitted power 
between the network access point and surface nodes (rovers or sensors), omnidirectional 
antennas placed at a l-2m heights, and 2.4GHz or 5GHz, as appropriate for the standard, 
carrier frequency. 

We assume that the standard commercial 802.11 wireless networks are being used. Rel- 
evant parameters for the 802.11a and b standards axe given in [5] and are summarized in 
Table 2.1. The 802.11 protocols contain functions related to connection management, link 
reliability management, and power management. The protocol developers are also looking 
into functions related to link security as well. The protocol supports two basic modes for 
mobile terminal operation: infrastructure mode and ad hoc mode. In the infrastructure 
mode used in this research, the mobile terminals access a wired network via one or more 
access points. In the ad hoc mode, the mobile terminals can communicate directly with each 
other without the need for the access point. 
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Table 2.1: Important 802.11 protocol parameters. 


Parameter 

802.11b 

802.11a 

Operating Frequency 
Maximum Data Rate 
Slot Time 

Minimum Contention Window 
Maximum Contention Window 
Modulation Format 

2.4 GHz 
11 Mbps 
20 fis 
31 

1023 

Complimentary Code 
Keying (CCK) 

Direct Sequence Binary 
Phase Shift Keying 
(DS-BPSK) 
Quadrature Phase 
Shift Keying (QPSK) 
Packet Binary Convolutional 
Code (PBCC) 

5 GHz 
54 Mbps 
9 /xs 
15 

1023 

Orthogonal Frequency 
Division Multiplexing 
(OFDM) 


The wireless protocol needs to be configured so that the medium access layer works to 
avoid collisions and not just recover from channel collisions. The basic method is to use a 
modified version of the IEEE truncated binary exponential backoff algorithm. If the channel 
is busy, the sender backs off an integer multiple of the slot time. The backoff algorithm 
has variations from the strict application to allow a more fair access of the channel by 
those stations currently experiencing a backoff wait. In order to avoid the hidden terminal 
problem, the 802.11 networks can use a request-to-send/clear-to-send mechanism for channel 
access control and collision avoidance. Because the wireless channel can have relatively poor 
performance, e.g. a channel BER of 10 -4 , there is a provision in the 802.11 protocols to 
allow for fragmentation of the packets. A shorter packet will have a higher probability of 
successful transmission so packet fragmentation can result in higher throughput in poor 
channels. Wireless protocols allow for power management through placing devices into sleep 
modes. When a device is in sleep mode, it cannot send or receive data. The device can have 
its sleep/awake duty cycle adjusted for specified power management goals. When a node 
needs to send data to a sleeping mode, the data sender must buffer the data until it detects 
the sleeping node has reactivated and is ready to receive. 
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3 Methods, Assumptions, and Procedures 

3 . 1 Introduction 

In this section, we first describe the methods, assumptions, and procedures used in validating 
the power delay profile results using local sites (this work was not included in the 2004 annual 
report but it did continue into Year 2 of the effort). Next, we describe the main work of 
Year 2 of the research project, namely simulating and analyzing the physical and medium 
access layers of IEEE 802.11a/b. 

3.2 Validation of Local Power Delay Profile Simulations 

In order to gain confidence in the simulation of a Power Delay Profile (PDP) for any environ- 
ment, we attempted to measure a detailed channel impulse response at 2.4 GHz at a location 
where the map data was available and where detailed measurements could be made. How- 
ever , due to the high-speed instrumentation required to capture a complete, high-resolution 
response, we instead decided to measure the RF power for a simple reflection off a large-scale 
topographical feature referenced to a direct path power level. This reflection could then be 
checked against simulation for proper time delay and power level as a validation. In order 
that the channel response would not be too complex, we selected local sites which were very 
flat, had minimal vegetation, and which had a single large sloping or vertical rock face. 

3.2.1 Local Site Locations for Reflection Measurements 

We chose two local sites for reflection measurements: 1) on the east side of Tortugas Moun- 
tain (approximately 10 minutes from the NMSU campus) and 2) Dripping Springs State Park 
(approximately 30 minutes from the NMSU campus). The “reflector” at the Tortugas Moun- 
tain site is a portion of the 200 m sloping, east face of the mountain (see Fig. 3.1 (a))and 
at the Dripping Springs site is a near vertical rock formation north of the visitor center 
(see Fig. 3.1 (b)). DEMs for these local sites are shown in Fig. 3.2. For these two sites, 
the coordinates for the transmitter (TX) and receiver (RX) positions are given in Table 3.1 
and the basic arrangement of transmitter, receiver, and reflector is shown in Fig. 3.3. For 
the Tortugas Mountain site, the distance travelled along the reflection path is 424 m which 
translates to a travel time of 1.4 /xs at the speed of light. 

3.2.2 Measurement Procedure 

In order to measure the reflected RF power off of Tortugas Mountain, we used directional 
antennas (Cushcraft 24012P) to transmit and receive the 2.4-GHz signal. Directional anten- 
nas were required in order to ensure that the received signal was strictly from the reflection 
off of Tortugas Mountain and not directly from a side lobe of the transmitting antenna. The 
actual antenna patterns were measured in the NMSU anechoic chamber and this pattern 
was later used in the simulation. 


6 


I 

> 

I 

I 


3 METHODS, ASSUMPTIONS, AND PROCEDURES 





3 METHODS, ASSUMPTIONS, AND PROCEDURES 




(a) (b) 

Figure 3.2: Digital elevation maps for (a) Tortugas Mountain and (b) Dripping Springs. 
Approximate locations of TX, RX, and reflection points are denoted on DEMs. 


Table 3.1: GPS coordinates for TX and RX positions at local sites. 
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Reflector 



Figure 3.3: Arrangement of transmitter (TX), receiver (RX), and reflector at (a) Tortugas 
Mountain and (b) Dripping Springs. 
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For the measurements at Tortugas Mountain, the transmitter antenna was mounted on 
top of the NMSU telemetry van and the receiver antenna was mounted on a university truck 
as shown in Figure 3.4 in order to have them each at roughly the same elevations. For the 
measurements at Dripping Springs, both transmit and receive antenna were mounted on 
poles approximately 3 m off the ground. The receive antenna is shown in as shown in Figure 
3.5. A Microdyne (TSS-2000) signal generator was used as the 2.4-GHz, 100-mW source 
and an Agilent spectrum analyzer (1544 03 B) was used to measure power levels at the receive 
side. 

First, the noise floor at the site was found with the spectrum analyzer. Then the trans- 
mitter was turned on generate a continuous wave signal with a transmission power of 20 dBm 
(100 mW) at 2.4 GHz. Second, the direct path measurements were made by pointing the 
transmit and receive antenna towards each other. Third, the transmit and receive antennas 
were pointed towards a small target on Tortugas Mountain and RF power measurements 
were made. In order to determine the level of sidelobe power propagating from transmit 
to receive antenna, both antennas were pointed up into the sky and sidelobe power was 
measured. 

3.2.3 Simulation of PDPs 

The DEMs for the Tortugas Mountain and Dripping Springs sites were converted to ATDI’s 
proprietary format for use with ICS Telecom. The procedure for generating the PDPs using 
ICS Telecom was given in Appendix C of the 2004 Annual Report. 

3.3 Measurement of Power Delay Profiles with the YellowJacket 
802.11b PLUS 

In June 2004, we acquired a YellowJacket 802.11b PLUS measurement instrument from 
Berkeley Varitronics Systems (BVS). In addition to measuring many IEEE 802. lib-specific 
parameters, it is also able to measure a power delay profile from an access point to the 
location of the YellowJacket. The PDP measurement can be geocoded using the on-board 
GPS unit and logged to a file. Such measurements provide us a capability to develop real- 
world, space-time RF channel models at 2.4GHz. In this section, we describe how that 
measurement is made. 


3.3.1 Background Information 

Consider a channel with impulse response h(n) with an input (transmitted) random process 
u(n) and an output (received) random process y{n). Compute the cross-correlation of the 
received process with the transmitted (£/[*] is the expectation operator) 


r uy (n) = E[u(k)y*(k -n)] 


= E 


u(k) ( h* (l)u(k — n — l) 

\Z=— oo i 
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Figure 3.4: Mobile antenna placement, (a) Transmit antenna on the top of the NMSU 
telemetry van and (b) Receive antenna on university truck at the Tortugas Mountain (A- 
Mountain) site. 
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mem. 


Figure 3.5: Receive antenna at Dripping Springs. 


= E u(fc) ^2 h(l)u* (k — n — l) 

l 1——OC 

= i? [ h(l)u(k)u* (k — n — l) 


= h(l)E [u(k)u*(k — n — l)] . (3.1) 

1— — OC 

All 802.11b packet transmissions start with several repetitions of the unmodulated Barker 
code. The 802.11b Barker code is 11-bits long and it has excellent correlation properties, i.e. 

r u (n) = 6(n). (3.2) 

The YellowJacket correlates the received 802.11b packet, y(n) against the Barker code at 
the start of the received packet. Thus substituting (3.2) into (3.3) shows that the cross- 
correlations are in fact the impulse response values 


r„,(n) = £ h(l)6(n - () 

Z=— 00 

= h(n). 


The PDP is readily computed from (3.3) by 


P'(n ) = max|/i(n)| 2 


Finally, YellowJacket normalizes (3.4) so that 


T 
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3.3.2 Relative Correlation Window 

Equation (3.4) is computed in the YellowJacket at a time spacing of 1/4 chip (91 ns), for a 
span of about 10 chips, i.e. and displayed in the relative correlation window. YellowJacket 
centers the PDP so that the peak (assumed to be the direct path) is at 0 (hence the relative 
correlation). Thus for convenience, the relative correlations start at 2 chips and end around 
+8 chips. This is illustrated in Figure 3.6 for the YellowJacket unit where the display is set 
to show signal correlation and received signal strength (RSSI). 



Figure 3.6: YellowJacket display showing the correlation display window plus a Received 
Signal Strength Indicator (RSSI) window. 


3.3.3 Multipath Values in Log File 

The first N = 22 relative correlations which are displayed are logged to a file. Presumably 
not all 40 values are stored in order to reduce memory. Thus the time extent of the correlation 
function is 


r = 22 x 91ns 
= 2 [is. 


(3.6) 
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3.3.4 RMS Delay Spread 


The rms delay spread, cr is defined as 

<7 = (3-7) 

where E is the expectation operator and P(n) is the received power at time n. YellowJacket 
computes and displays the rms delay spread using the 22 PDP values as follows. Since we 
have only a finite amount of data the expectations must be estimated. The second moment, 
E [ P 2 (n)] is estimated with 


E[P 2 (n)} 


E^P(n)n 2 

El-o 1 P (n) n 


(3.8) 


where N = 22 is the number of PDP samples YellowJacket has stored. The first moment or 
mean delay , E [P (n)] is estimated with 


E [P (»)] 


S£o P(n)n 
E^P(n)n 


(3.9) 


Using these estimates, the rms delay spread is calculated as 

a = T ] /I[P 2 (n)]- (P[P(n)]) 2 . (3.10) 


where T = 1/22 MHz (1/2 x 1/4 chip period?). 


3.3.5 Final Notes 

The YellowJacket provides a convenient method of measuring real-world PDPs provided the 
delay spread is not too large so that the values alias. Note that if the propagation delay 
between the strongest path and multipath reflections is larger that 9 chips, the multipath 
can “alias” to the start (-2) of the multipath display. 

3.4 Simulation of IEEE 802.11a/b Physical Layer using MATLAB 

The work in this section is based upon the paper “Performance Evaluation of the IEEE 
802.11a and b WLAN Physical Layer on the Martian Surface” by Deva K. Borah, Anirudh 
Daga, Gaylon Lovelace and Phillip DeLeon. This paper was presented at 2005 IEEE 
Aerospace Conference. 


3.4.1 Summary 

The performance of IEEE 802.11a and b WLAN standards on the Martian surface is studied 
in this task area. The Gusev Crater region and the Meridiani Planum (Hematite) region 
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are chosen as example sites based on the mission science and mission success criteria. The 
radio frequency multipath environment is obtained using digital elevation maps from the 
Mars Global Surveyor mission, taking into account the atmosphere and other factors on 
the Martian surface. It is observed that IEEE 802.11a performs well in terms of packet 
error rates at distances up to a few hundred meters from the transmit antenna when the 
transmit power is 1 W and the antennas are located 1.5 m above the ground. Although the 
performance of IEEE 802.11b is found to be more adversely affected, its performance too 
can be improved significantly using a RAKE receiver. It is observed that the lower data 
rate modes of 802.11a show much better results in terms of bit error rates. However, both 
802.11a and b appear to provide effective communications within a few hundred meters of 
the transmitter in the selected sites considered. 

3.4.2 Introduction 

Future space missions on the Martian surface may involve multiple rovers collecting data at 
different locations, and communicating wirelessly with common access points. Such com- 
munications have to be reliable, robust and power efficient. Development and testing of 
such communication technologies from scratch is an expensive proposition. A more cost 
effective approach would be to adapt existing technology with appropriate modifications. 
Towards this objective, this paper investigates the physical layer performance of two well 
known Wireless Local Area Network (WLAN) standards IEEE 802.11a and IEEE802.11b 
under the Martian environment, and identifies the issues that need to be addressed. 

The IEEE 802.11b standard [6] provides data rate options of 1, 2, 5.5 and 11 Mbit/s in 
the 2.4 GHz band. The modulation options include direct sequence spread spectrum using 
differential Binary Phase Shift Keying (BPSK), Quadrature Phase Shift Keying (QPSK), 
Complementary Code Keying (CCK), and Packet Binary Convolutional Code (PBCC). Al- 
though primarily designed for indoor office environments, recent studies have shown good 
performance of 802.11b in outdoor environments [7], [8], However, the performance with low 
height rover antennas on the Martian surface, and the performance comparison of 802.11b 
with respect to 802.11a in the Martian environment are important issues that have not been 
addressed before and need investigation. 

The IEEE 802.11a standard [9] operates in the 5-GHz band and uses the Orthogonal 
Frequency Division Multiplexing (OFDM) technology. It can support data rates of 6, 9, 
12, 18, 24, 36, 48, and 54 Mbit/s. The standard employs convolutional encoder, and uses 
cyclic prefix of 0.8 micro second duration. This enables it to handle the multipath problem 
more successfully [10], [11]. However, longer delay spreads, which can happen on the Martian 
environment with low height antennas and longer transmitter/receiver distances, can severely 
affect its performance. Therefore, the effects of such delay spreads on the Martian surface 
require investigation. 

In previous work [12], the RF environment of the Martian surface has been extensively 
studied. In particular, the RF coverage patterns produced from a 2.4 GHz transmitter with 1 
W radiated power and 1 m antenna height within Gusev Crater and Meridiani Planum have 
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been investigated. These simulations use 11 m/pixel digital elevation maps from the Mars 
Global Surveyor mission. The software used in this study takes into account the propagation 
factors such as planetary radius, atmospheric density and composition, soil chemistry, etc. 
The impact of surface clutter (rocks) on RF propagation has also been examined. It has 
been observed in that study that while significant terrain variation can have a major impact 
on the coverage, sufficient RF signal power for an IEEE 802.11b link is possible at these sites 
over several kilometer distances even with low antenna heights. 

This analysis uses the received power results from [12] and recent results regarding the 
simulation of the multipath environment in the performance evaluation of the 802.11 a and b 
standards. Here, we will study the performance of different data rates for different transmit 
and receive antenna locations and several sites on Mars. It is observed m the results that 
multipath environments can severely affect the performance of 802.11b. The use of receivers 
that mitigate multipath effects, such as RAKE-type receivers, is found to provide significant 
improvement. The performance of 802.11a is also found to be affected by the multipath 
environment, especially in the absence of clear line-of-sight. In particular, the higher bit 
rate modes of IEEE 802.11a are found to be more affected by the multipath effects. Further, 
when the delay spread exceeds the 0.8 fis cyclic prefix duration, the performance drops 

rapidly. , 

T his study section is organized as follows. In Section 3.4.3, we provide an overview of 
the 802.11a and b physical layer (PHY) specifications. The packet structures described in 
this section are faithfully simulated in our simulation results section. In Section 3.4.5, the 
site selection criteria on the Martian surface, and the radio frequency multipath calculation 
are described. In Section 3.4.4, the methodology for preparing MATLAB and the product 
mLAN for the physical layer simulations is described. Section 4.6 presents the simulation 
results in terms of packet error rates and bit error rates. A discussion on the interpretations 
of the simulation results is presented in Section 4.6.6. 

3.4.3 IEEE 802.11 WLAN 

In this section, we present a brief overview of the physical layer specifications of IEEE 802.11a 
and b standards. 

IEEE 802.11 

The IEEE 802.11a physical layer is based on orthogonal frequency division multiplexing 
and operates in the 5-GHz band providing data payload capabilities of 6, 9, 12, 18, 24, 36, 48 
and 54 Mbit/s. The different transmission rates are obtained by varying the modulation type 
and/or the channel coding rates. The system uses 52 subcarriers that are modulated using 
BPSK, QPSK, 16- or 64- Quadrature Amplitude Modulation (QAM). The error correction 
coding uses a convolutional encoder with a coding rate of 1/2, 2/3 or 3/4. 

The Physical layer Protocol Data Unit (PPDU) format is shown in Fig. 3.7. The Physical 
Layer Convergence Procedure (PLCP) preamble field consists of 10 repetitions of a short 
training sequence, and two repetitions of a long training sequence preceded by a Guard 
Interval (GI). A single BPSK encoded OFDM symbol follows. It contains a 4-bit RATE 
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PLCP Preamble 

SIGNAL 

DATA 

1 2 symbols 

One OFDM symbol 

OFDM symbols 


(a) IFFF 802.1 ] a 


— 
PLCP Preamble 

— 
PLCP Header 

PSDU 

144 bit or 72 bit 

48 bit 



(b) IEEE 802.11b 


Figure 3.7: PPDU Format for IEEE 802.11a and b. 


field, a 12 bit LENGTH field, one reserved bit, one parity bit and 6 ‘zero’ tail bits encoded 
with a rate 1/2 convolutional code. The DATA portion contains a 16 bit SERVICE field, 
a Physical Sublayer service Data Unit (PSDU), 6 ‘zero’ tail bits and pad bits, and may 
constitute multiple OFDM symbols. 

The data to be transmitted are scrambled to remove any spectral line from the data. 
They are then convolutionally encoded with a rate-1/2 encoder with generator polynomials 
g 0 = 133 g , gi = 171s, and puncturing is performed if necessary. All encoded data bits 
are interleaved using two steps. First, consecutive coded bits are mapped to non-adjacent 
subcarriers. The second step maps consecutive coded bits onto the less and more significant 
bits of the constellation. 

The OFDM symbols are transmitted using a relatively long cyclic prefix of duration 
T G i = T F ft/ 4, where T fft is the duration of an OFDM symbol. The duration T fft equals 
3.2 /zs. Thus the symbol interval is 4.0 /zs. The PLCP preamble duration is 16 /zs, and the 
SIGNAL symbol lasts 4.0 /zs. 

IEEE 802 . 11 b 

The IEEE 802.11b Direct Sequence Spread Spectrum (DSSS) can provide data rates of 1, 
2, 5.5 and 11 Mbit/s in the 2.4 GHz band. The basic data rate of 1 Mbit/s is provided using 
Differential Binary Phase Shift Keying (DBPSK) while the 2 Mbit/s rate uses Differential 
Quadrature Phase Shift Keying (DQPSK). The above two data rates employ 11 chip long 
Barker sequences for spreading with a chip rate of 11 MHz. 

Higher data rates of 5.5 Mbit/s and 11 Mbit/s are available in 802.11b through the use 
of complementary code keying at the same chipping rate of 11 Mchips/s. Each CCK symbol 
consists of 8 complex chips: e^ 1+ ^ 3+< H e^ + ^>, e j(4n + <h^)^ 

e j(4i+<f>3) , _gj(<Ai+<fo/ j n case Q f 5.5 Mbit/s, 4 bits are transmitted per symbol while 
in the case of 11 Mbit/s, the number of bits transmitted per symbol is 8. The first two bits 
are used to compute a phase change for <f) 1 with respect to phase <f>i of the preceding symbol 
or the phase of the preceding DQPSK symbol if there is a header to PSDU transition. In the 
case of 5.5 Mbit/s, the remaining two bits are used to derive the phase fo and (f> 4 , while 
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the 11 Mbit/s mode uses the remaining 6 bits to compute <j> 2 , <f> 3 and (j) 4 based on QPSK. 
An optional mode replacing CCK modulation with packet binary convolutional coding with 
a 64-state encoder is also available. 

The PPDU format for TF F.F, 802.11b is also shown in Fig. 3.7. Two different preamble 
and headers are defined: long PLCP PPDU format and short PLCP PPDU format. The long 
format cont ains a 144-bit preamble and a 48-bit header, while the short format contains a 72- 
bit preamble and a 48-bit header. The preamble contains two fields: synchronization (Sync) 
and Start Frame Delimiter (SFD). The Sync field is provided to enable the receiver per- 
form necessary synchronization operations. The SFD indicates the start of PHY-dependent 
parameters within the PLCP preamble. The header consists of signal, service, length and 
Cyclic Redundancy Code (CRC) fields. The signal field indicates the data rate that is used 
for the transmission and reception of the PSDU. The service field contains 8 bits, and they 
carry some information about modulation, symbol clock etc. The length field indicates the 
number of microseconds required to transmit the PSDU. Finally, a 16 bit CRC protects the 
signal, service and the length fields. The long PLCP preamble and header axe both trans- 
mitted using 1 Mbit/s DBPSK modulation. In the case of a short PLCP, the preamble is 
transmitted using 1 Mbit/s while the header is transmitted using 2 Mbit/s. The transmitted 
data bits are scrambled at the transmitter and descrambled at the receiver. 

3.4.4 MATLAB Simulation Development 

We use the capabilities of MATLAB to model the physical layer in this investigation. MAT- 
LAB does not have the required tools in its manufacturer’s distribution to perform all of 
the analysis. To augment MATLAB, we acquired the mLAN tools from CommAccess. This 
section describes the effort required to male these tools usable in the analysis and simulation. 

IEEE 802.11a 

The IEEE 802.11a physical layer simulation software was developed around the mWLAN 
Toolbox by CommAccess Technologies. This package provides a large set of MATLAB func- 
tions, most of which are sub-functions that need not be called directly. They provide a 
reasonably complete implementation of the 802.11a OFDM signaling and packet construc- 
tion. In brief, a block of random data is generated and passed to mWLAN’s PLCP_lla() 
function, which returns a complex baseband transmit vector, sampled at 20MHz, represent- 
ing one packet. After channel effects (noise, fading, multipath, etc) are imposed, the vector 
is passed to parse_signal_lla(), which extracts and checks the transmission parameters 
encoded in the packet header SIGNAL symbol. An FFT translates back from the time to the 
frequency domain, and the demodulation and data decoding are performed by the mWLAN 
functions DeMod_lla() and DeMod_post_lla(). 

A significant amount of “glue logic” code is required to piece together these mWLAN 
functions. To simplify the effort, ensure repeatability, and guard against mistakes, a top- 
level simulation script (sim_ER_lla) was developed. It reads experiment parameters (e.g. 
packet size, data rate, E b /N 0 , number of packets to simulate, etc) from a file, so that many 
different simulation experiments could be run without modifying the simulation source code 
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itself. Simulation results were logged into result files, and additional scripts were developed 
to process the results and produce figures and tables. 

In addition to the mechanics of simulation (looping through variations and iterations, 
calling mWLAN functions, storing results, etc), our code implements whatever channel ef- 
fects are specified in the parameter file. The mWLAN toolbox does not provide any channel 
modelling functionality. Our software reads in a power delay profile (calculated from Martian 
terrain maps by ATDI’s ICS Telecom RF modelling software). Rayleigh fading is applied, 
and the result is convolved into the complex baseband transmit vector for the packet. The 
receiver end of the simulation estimates the channel by a linear least-squares method, and 
weights the output of the FFT of each received OFDM symbol. 

To validate the behavior of these tools, the results of our simulations were compared 
to published results. Doufexi [lljpresented BER-vs-Et/Ao results (their Figure 11) for a 
channel with AWGN plus a 50ns exponentially decaying multipath power delay profile. We 
ran our simulation with this channel model for the 6, 12, 36 and 54 Mb/s data rates, and 
obtained very good agreement with the published curves. Comm Access Technologies publish 
pure AWGN results in the User Manual for their eWWoB product (a GUI workbench based 
on mWLAN). Their results are consistently 1.25dB (i.e. a factor of 64/48) better than ours. 
After extended discussions with Dr Song An, who developed mWLAN, we found that the 
difference is exactly attributable to the treatment of noise power in the non-data OFDM 
sub-bands, in our respective definitions for Eb/No ■ The distinction is somewhat arbitrary, 
and in the end, we decided to follow the convention used in Doufexi’s published results, 
which derives from what seems a more intuitive handling of noise power by the DFT of the 
OFDM receiver structure. 

IEEE 802 . 11 b 

The physical layer simulations for IEEE 802.11b were performed using the mWLAN 
toolbox developed by CommAccess Technologies, Inc. The mWLAN toolbox provides several 
modules associated with the transmission and reception of data at the physical layer. The 
modules available in mWLAN include: 

• Bit Generation module to generate random data bits, 

• Square root raised cosine filtering module for pulse shaping (oversampling) the raw 
data for a given rolloff factor, 

• CRC encoding module to encode certain fields in the physical layer header with the 
CCITT-CRC 16 polynomial, 

• Scrambler and descrambler modules to scramble and descramble the data bits in ac- 
cordance with IEEE 802.11b standard specification, 

• Modulation and demodulation modules for the basic data rates of 1 Mbps(DBPSK) 
and 2 Mbps(DQPSK), 
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• Modulation and demodulation modules for the high data rates of 5.5 Mbps and 11 
Mbps(CCK), 

• Modulation and demodulation modules for the PBCC mode (all data rates), and 

• CRC check module to verify the validity of the received data. 

All these modules are available in the form of Windows Dynamic Linked Libraries(DLLs) 
and can be accessed using their associated MATLAB functions. Bit error rate and Packet 
error rate simulations for all the data rates can be performed by writing MATLAB programs 
which use these functions in a logical order. 

The original module for 1 and 2 Mbps modulation had to be changed because it carried 
out DBPSK and DQPSK encoding after spreading the data with the Barker sequence. The 
correct way to go was to Barker spread the DBPSK and DQPSK encoded data to preserve 
the autocorrelation properties of the Barker sequence. Corresponding changes were made in 
the demodulator module as well. Also, the original module provided for only sub-optimal 
demodulation of the CCK based 5.5 and 11 Mbps modes. New modules were developed to 
perform near-optimal demodulation of the 5.5 and 11 Mbps modes. 

After making these changes and additions, BER simulations were carried out in AWGN 
channel to verify the correctness of the modules. The results for the 1 and 2 Mbps schemes 
agreed very closely with the results published in [13]. The results for the 5.5 Mbps and 11 
Mbps schemes agreed very closely with the results published in [14]. 

Several new modules were developed to incorporate the Martian multipath power delay 
profiles(PDPs) generated by ICS Telecom. The PDPs obtained from ICS Telecom had 
a 10ns resolution and had to be fractionally resampled to match the sample time of the 
oversampled (square root raised cosine filtered) data. The power delay profiles were converted 
to normalized channel impulse responses to observe the effect of multipath on BER vs Eb/No 
for various data rates. Modifications were also made to observe the BER and PER as 
a function of distance. Least squares channel estimation was carried out to estimate the 
multipath channel impulse response at the receiver. RAKE receiver was implemented, which 
significantly improved performance. 

3.4.5 RF Environment on the Martian Surface 

We have used the ICS Telecom software from ATDI [15] to obtain the multipath environment 
on the Martian surface. DEM files are converted to ATDI’s format for the Martian sites 
(llm/pixel resolution), and are loaded into the software. 

Site Selection 

We have selected the Gusev Crater and the Meridiani Planum (Hematite) regions [16] 
as example sites for our study. These two regions are chosen considering the mission science 
and mission success criteria [16], [17]. The mission science criteria included evidence of water 
on the Martian surface in the past. The Gusev Crater appears to have been a lake fed by a 
river at one t im e The Meridiani Planum region shows the chemical signature of Hematite 
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minerals associated with ancient water locations. For mission success, the sites are chosen 
“near the equator, low in elevation, not too steep, not too rocky, and not too dusty in 
addition to other factors. The locations of the selected sites are shown in Table 3.2. 


Table 3.2: Sites for WLAN Performance Study. 


Site 

Mars Latitude 

Mars Longitude 

Gusevl - Site 1 

14° 47' 39.35" S 

176° T 29.18" E 

Gusevl - Site 2 

14° 58' 41.95" S 

176° 2' 53.51" E 

Gusevl - Site 3 

15° 11' 35.66" S 

176° 4' 31.23" E 

Hematite4 - Site 1 

2° 11' 0.69" S 

-5° 53' 5.16" E 

Hematite5 - Site 1 

1° 52' 29.16" S 

-5° 25' 39.59" E 


RF Model 

The Irregular Terrain Model (ITM) has been used. It is a general-purpose propagation 
model for frequencies between 20MHz and 20GHz. This model predicts the median attenu- 
ation of a radio signal as a function of distance and the variability of the signal in time and 
space. The predictions are based on electromagnetic theory and statistical analysis of both 
terrain features and radio measurements. 

The ITM source code has been modified for Martian parameters. Atmospheric attenua- 
tion is negligible — actual calculations for a horizontal path on Mars’ surface yield attenua- 
tion of approximately 10“ 6 dB/Km at 2.5GHz [18]. The ITM source code for propagation 
on Earth accounts for atmospheric refraction by introducing an “effective radius” multiplier 
of K = 1.33. The effective radius used for Earth is K times Earth’s physical radius. Mars’ 
atmosphere is so diffuse, even at the planet’s surface, as to resemble a vacuum compared to 
Earth’s. Thus we assume atmospheric refraction is negligible in our study [18], [19]. We set 
K = 1, and use an effective radius equal to Mars’ physical radius. We note that in some 
implementations, an effective curvature (inverse of the effective radius) is used. 


3.5 Simulation of IEEE 802.11a/b MAC Layer using OPNET 

During the summer of 2004, we considered alternatives for evaluating the MAC-layer perfor- 
mance of wireless LAN protocols on the Martian surface. Ideally, such an evaluation would 
be tightly coupled with the physical-layer model developed during the initial phase of the 
research. Simulation software tools which we considered included OPNET, QualNet, NS2 
and various MATLAB toolboxes and packages. We chose OPNET because of its ease of use 
and wide user base. We obtained three licenses and began training during the first summer 
session. 

We discovered that OPNET has limited capabilities in incorporating a physical-layer 
model in a MAC-layer simulation, so we investigated alternatives for extending its capability 
in this area. Two basic approaches were considered: 
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1. Find a reasonable, parameterized physical-layer model, based on previous studies, and 
incorporate this into the OPNET simulator. 

2. Program the physical-layer model using MATLAB, and then provide an interface be- 
tween this tool and OPNET in order to perform a cosimulation. 

The first approach would be preferable if an appropriate physical-layer model could be 
characterized using a sufficiently small number of parameters. The standard software mou- 
ules used by OPNET for simulating a wireless channel employ a very simple model, where 
the probability of packet error is related straightforwardly to the signal-to-noise ratio [20, 
Ch. 11]. These modules would have to be modified, so that significant effects such as mul- 
tipath are included in the packet-error computation along with signal strength and noise 
level. The modification would be accomplished by providing user-defined process modules, 
written in C, according to a standard procedure supported by OPNET [21], This approach 
would be simpler than a full-up cosimulation, and the simulations would probably run much 
faster. The reason for this is that OPNET views the transmission or reception of a packet as 
a discrete event. In a cosimulation, each of these events involving a wireless link would give 
rise to a physical-layer simulation requiring a much finer time granularity, thus multiplying 
the simulation time. 

If a parameterized model proves to be insufficient, then some form of cosimulation will 
have to be used. To provide an appropriate interface, two options were considered: 

1. Use OPNET External-Domain cosimulation tools to exchange data between OPNET 
and MATLAB [22], 

2. Use MATLAB tools for calling the MATLAB engine from the C code used to program 
the OPNET process modules [23]. 

The OPNET External-Domain tools proved to be complicated, difficult to use, and not 
well documented. Further, we were not able to find any literature referring to an OPNET- 
MATLAB cosimulation that was carried out this way. By contrast, procedures for calling the 
MATLAB engine from a C-language program are very well described in literature provided 
by Mathworks, and this approach for linking MATLAB and OPNET has been successfully 
used by other investigators studying the interaction between MAC and physical layers (see, 
e.g. Dham [24]). 

The approach makes use of the MATLAB Engine Library which contains a set of routines 
that can call MATLAB from a C program. We can write C programs called MATLAB 
Engine Programs that can open and close MATLAB, send data to and from MATLAB, and 
send commands to be processed in MATLAB. Thus, we may use the Engine Library to 
pass information between OPNET and MATLAB via the C-language user-provided process 
modules discussed above. 

Here is a detailed procedure for preparing an OPNET-MATLAB cosimulation using this 
approach on Windows XP: 
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1. Ensure that Microsoft Visual C++ is installed on the computer 

2. Open the OPNET Modeler software. Click on Edit— ^Preferences 

3. The preference bind_shobj_prog should have the value bindjso_msvc 

4. The preference bind_static_prog should have the value bind_msvc 

5. Specify the path where the MATLAB Engine Library files are located. This can be 
done by editing the bind jshobj -flags preference. Change the value to 

/LIBPATH::C:\MATLAB6p5\extern\lib\win32\microsoft\msvc60 

6. To use the MATLAB engine the following library files have to specified. 

• libeng.lib 

• libmex.lib 

• libmx.lib 

• libmat.lib 

T his can be done by changing the value of the bind_shobj Jibs preference to libeng.lib 
libmex.lib libmx.lib libmat.lib 

7. Close the preferences dialog box 

8. Open the Process Model by clicking on File— ►New— ^Process Model 

9. Define the state variables by clicking on the SV icon and the temporary variables by 
clicking on the TV icon. 

10. Click on the HB icon and specify the relevant header files. 

11. Click on the Create State icon (the icon at the top left corner) and place a state in 
the workspace 

12. Enter the C program that will call MATLAB by clicking either on the top half of the 
state (Enter execs) or the the bottom half of the state(Exit execs) 

13. Compile the program to check for compilation errors 

14. Click on Intefaces— > Process Interfaces. Change the initial value of begsim intrpt 
attribute to enabled 

15. Associate this process with a node and the node with a network. 
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3.6 Outdoor IEEE 802.11b Measurements 

In order to validate the IEEE 802.11b physical and MAC-layer simulations, actual measure- 
ments of IEEE 802.11b must be made in the outdoor environment. These measurements 
include data rates, bit error rates (BERs), packet error rates (PERs), power delay profiles, 
signal strength, etc. For the validation work, we considered three environments: ffee-field, 
simple reflection, complex terrain. For the free field environment, we conducted experiments 
at the football stadium parking lot on the NIvISU campus. The lot has a gravel surface, 
is very flat and has virtually no vegetation or man-made objects. The nearest man-made 
objects axe at least 200 m away from our staging area. The simple reflection environment 
consists of a large (9' x 6') steel reflector located on the football stadium parking lot. Fi- 
nally, we have chosen the Dripping Springs site (La Cueva rock formations) for measuring 
the IEEE 802.11b protocol performance in a complex terrain. Figure 3.8 shows the measure- 
ment equipment on the parking lot and Figure 3.9 shows the reflector. Figure 3.10 shows 
the topographic map of the setup at La Cueva (rock formations) at Dripping Springs. The 
complete field test procedure is given in Appendix E. 

The measurement setup for the free-field experiments includes a basic IEEE 802.11b 
wireless access point/router elevated to either 1 m or 2 m height with a laptop hardwired 
to it. In addition both a laptop (with an IEEE 802.11b wireless card) and the Yellow- 
Jacket are located at various distances along a radial path from the access point. The entire 
setup is portable and battery powered. A dataflow from either Transmission Control Proto- 
col/Internet Protocol (TCP/IP) or User Datagram Protocol (UDP) based source is initiated 
between the laptops using the network bandwidth measurement software application Iperf 
that is described below. During the dataflow, Iperf can report the data rate and, in the case 
of UDP, datagram errors. The YellowJacket (YJ) is able to log data such as relative signal 
strength and relative correlations (power delay profile). In addition, YellowJacket logfiles 
can be geocoded with location and time information. Finally, a set of MATLAB programs 
were developed to aid in the visualization of YellowJacket logfile data. A short User’s Guide 
is given in Appendix D. 

Considerable effort was made to disable all power-saving features on the laptops and to 
eliminate all unnecessary processes (both operating system and applications) so that sus- 
tained and stable data transfer could be made. Without the power-saving features disabled, 
the results found with Iperf indicated that the data transfer rates were very erratic because 
the wireless card’s RF interface would be cycled on and off. 

The use of commercial 802.1 lx hardware was an expedient and relatively inexpensive way 
to verify the performance of this protocol. The problems such a wireless network faces in a 
terrestrial environment can in some cases be translated to a simulated Martian environment. 
The primary application of 802.1 lx wireless, in buildings and public places, is not directly 
applicable and some of the aspects of the design are not best suited for field use (low power, 
omnidirectional antennas, sealed units, and AC power needs). 

The choice of Dell PC laptops was driven by the ATDI HerTZMapper software for Win- 
dows. This enabled the experimenters to use a single computer to run simulation models 
in the field and use the same computer for the data throughput measurements. Otherwise, 
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Figure 3.8: Measurements on the NMSU football stadium parking lot. In the foreground is 
a laptop hardwired to the wireless access point/router elevated to 1 m. In the background 
are personnel measuring RF data with the YellowJacket and data rates with another laptop. 
Flags are spaced every 20 m to a distance of 100 m. 
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Figure 3.9: (a) Reflector made of steel bookshelves and (b) setup for measurements in the 
simple reflection environment. 
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Figure 3.10: Arrangement of IEEE 802.11 access point (TX) and node (RX) at La Cueva 
(rock formations) in Dripping Springs. 

a Mac PowerBook G4 with built-in “Airport” was available. During our testing, we found 
that the performance of the Apple Airport card (firmware version 9.42) was not as good 
as the Intel PRO/Wireless 2100 card in the Dell Latitude D600. But the burden of the 
Microsoft Windows XP operating system on Mac users was a major factor in the time it 
took to establish reliable, repeatable measurements with our wireless network. 

Eventually, we had to establish an independent Wireless Local Area Network (WLAN) 
with no connection to the rest of the world or internet. We set up Wired Equivalent Privacy 
(WEP) to keep others from easily accessing our WLAN and adding mystery to our mea- 
surements. The Mac OSX has a utility called “Location Manager” that was a great help in 
changing the connections to the wireless Access Points (AP) or campus internet as needed. 
The XP OS did not have a similar capability (even on a laptop) and it was difficult, indeed, 
to switch from a wireless to ethernet link or return. A reboot was usually necessary. Note: 
This may have been a Windows issue that has since been resolved in Service Pack 2. 

The necessity of an independent, self-contained, WLAN seems obvious, but our initial 
field tests brought out even more insidious elements that required tighter control. Our first 
field test to establish the distance limits of performance with 802.11b showed wild variations 
in performance, not related to distance. The appearance of periodic deteriorations in transfer 
rate caused us to suspect that an internal operating system function on the Dell was taking 
CPU cycles. A suspicion of “spyware” or background network access software created a need 
to investigate the Windows XP OS in depth, (see, for example, [25] and [26] for extensive 
additional information on these topics.) 

There is a function of the XP OS that checks periodically for the availability of any type 
of network access, be it Modem, Wireless or Ethernet. This is because the ^Vindows OS back 
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to Microsoft.com for software updates or other internal requirements. Spyware will attempt 
to upload information gleaned from a user’s activities and report them on a periodic basis to 
an unannounced host. These effects can range from occasionally annoying to compromising 
the integrity of the computer. In either case, from the point of view of our experiment, 
they represent an uncontrolled network access and can affect the measurements made over 
the network by being a background noise source. Much effort was put in to disable these 
functions, and then to keep the PC laptop “sanitized” by never letting it connect to the 
general internet again. This complicated the exchange of data from the laptop to desktop 
computers in the lab, but we used a USB flash disk to transfer the necessary files and keep 
the laptop from being violated by connecting over an open Internet. 

Further study of how Windows implements the Internet Protocols led to the measurement 
of Bandwidth Delay Product which is a measure of network performance. This Bandwidth 
Delay Product is the product of multiplying the “bottleneck bandwidth” (slowest hardware 
link speed) by the Round Trip Time (RTT). RTT is most readily measured by a network 
utility called ping. With ping , an Internet Protocol (IP) address is probed and the response 
time (in milliseconds) is reported. From the NMSU campus internet, a ping to NASA Glenn 
takes about 80 ms. We found that a ping over our WLAN to the Dell, one meter away, took 
about 300 ms. Clearly, there was something in our link that was amiss. Testing all possible 
configurations of wired and wireless using the D-Link AP focused the bottleneck at the Dell 
laptop wireless card only. It seemed a daunting task to determine what, if anything, could 
be done to locate the level at which this link was failing. There are (at least) two setup 
or configuration locations for the wireless card: the Windows Network Connections and the 
Intel ProSet application (see Figure 3.11). Two settings are made available for control - the 
RF power and the device power management configuration. The default RF power is set to 
highest, but the power/performance default is set for lowest power (appropriate on a laptop). 
What is not immediately clear is that this setting is not merely battery power usage, but 
Intel’s implementation of Power Save Polling, an option in 802.11 to “turn off” a node when 
network usage is minimal. Not all AP’s recognize Intel’s implementation and apparently 
we were caught in that situation. It was fortuitous that a simple change of this setting 
improved the wireless performance of the Dell/Intel laptop from intermittent and unreliable 
to functionally usable, (see, for example, [27] for further information on this topic.) 

In order to create a stable, realtime, testbed for IEEE 802.11 hardware it was necessary 
to severely restrict the configuration of the Dell laptop, and to a lesser extent, the Mac 
PowerBook. For best performance, the only non-operating system application running would 
be the command window where Iperf executes. Simultaneously executing network monitors 
or Secure Shell functions would allow enough CPU resources to keep the wireless data stream 
uninterrupted in any non real-time operating system (Windows, Linux, or Unix) . 

We also found that the UDP protocol in Iperf takes all the available CPU cycles on the 
Client computer it is running from. This is addressed in a forum response from ncsa. In 
order to guarantee a continuous data rate, the Iperf program uses “busy waiting” to provide 
accurate timing and subsequently absorbs all CPU resources to provide the stream. As Iperf 
is a performance testing tool, it should properly be used temporarily and as the only running 
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Figure 3.11: Intel Preset WLAN configuration window. 


application [28]. 

Before finding Iperf, we had tried to transfer big data files over 802.1.b using AOL Instant 
Messenger, iTunes, and Quicktime. While the free Iperf tool might leave us wondering what 
we’re seeing, it is, at least, quite controllable. The first realization we came to was the use 
of the terms “client” and “server”; they seemed counter-intuitive. Whatever the history 
of networking terminology, in Iperf the “client” sends the UDP data TO the “server” and, 
therefore, the client should be connected at the wireless link under test. 

3.6.1 The Iperf Network Measurement Bandwidth Tool 

Complete information regarding the installation and use of Iperf can be found at 

http : //dast .nlanr.net/Projects/Iperf/ 

Iperf must be installed on two computers connected (wired or wirelessly) to a LAN or the 
Internet. One computer is designated as the server and the other as the client. Dataflow 
originates at the client and terminates at the server. For our setup, the laptop hardwired to 
the access point is the Iperf client while other laptop is the Iperf server. On UNIX systems, 
Iperf is executed by typing Iperf commands in a terminal. On Windows systems, Iperf is 
executed by typing commands in the command line window. 

An example of a TCP/IP data flow using Iperf is as follows. On one PC we type: 

iperf -s -i 1 

to denote this PC is the Iperf server and to update measurement information every one 
second. On the other PC we type: 
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iperf -c 10.0.0.1 -i 1 

to denote this PC is the Iperf client which is connecting to a server whose IP address is 
assumed to be 10.0.0.1, and to update measurement information every one second. 

An e xamp le of a UDP data flow using Iperf is as follows. On one PC we type: 

iperf -s -u -i 1 

to denote this PC is the Iperf server, UDP dataflow, and to update measurement information 
every one second. On the other PC we type: 

iperf -c 10.0.0.1 -u -b 5.5M -t 60 -i 1 

to denote this PC is the Iperf client whose IP address is assumed to be 10.0.0.1, a UDP flow 
at 5.5 Mbps for 60 seconds, and to update measurement information every one second. In 
the setup used in the outdoor measurements, we found that we could only sustain 5.5 Mbps 
without significant datagram errors in the free-field environment. This appears to be the 
baseline capacity of the wireless fink in our setup. 

The data output from Iperf can be redirected to a file for later analysis. Typical data is 
shown in Figure 3.12: 
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Figure 3.12: Example output from Iperf. 


3.6.2 YellowJacket and IEEE 802.11b 

The Berkeley Varitronics Systems YellowJacket PLUS b [29]is an 802.11b receiver wrapped 
around an HP (Compaq) iPAQ (model 5150 in ours) Personal Digital Assistant (PDA). The 
PLUS version includes a Global Positioning System (GPS) receiver and antenna jack. 

Software on the iPAQ provides an interface to the receivers and makes measurements of 
the 802.11 RF environment along with location data for the position of the YellowJacket 
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(YJ). The data can be logged on the iPAQ and then transferred to a PC (Windows only) via 
the built-in infrared (ir) port on the iPAQ. A BAFO brand ir/USB interface is provided by 
BVS. but we could never get it to work with our Dell Latitude laptop. Since the Windows 
OS gives several options to configure a new communications device, much frustration was 
expended just to establish communications with the YJ. We found that the built-in ir port on 
the Dell laptops was not enabled by the default OS and BIOS settings. With assistance from 
a systems engineer familiar with the BIOS, we finally established ir communications with the 
YJ. But, that link is quite slow (is it 9600 baud?) and there are a number of translations that 
must take place to view YJ data: a translation from ir to file using Microsoft ActiveSync, a 
translation from PDA format to “desktop” format, a translation from BVS format to Excel 
using B\ S Chameleon’ and finally using MS Excel to create a graph. We eventually wrote 
custom MATLAB software to analyze the YJ logs. 

Some of the features of the YJ don’t work as well as we would like, but it is nonetheless 
a unique and useful tool: 

1. The GPS receiver performance is mediocre compared to the handheld units meant for 
hiking or navigation. 

2. The screen snapshot function of the YJ only works for the Spectrum screen. 

3. The Marker function to update the log file does not work - this would have been a way 
to synchronize the YJ and Iperf data. 

4. The Relative Correlation screen only logs 22 data points out of the 44 measured. 

5. The Delay page does not log at all, although that data can be derived. 

6. The unit has an odd behavior with the RSSI display; this signal strength has spikes to 
a low value whenever the data rate increases toward the 11 Mbps rate. 

7. The log files slow their update as the data rate increases. At the management and 
housekeeping rate of 1 Mbps, the log updates around 8 Hz, but when data is streaming 
this slows to 1 Hz or less. 

8. The Relative Correlation display can easily “alias” or wrap- around when the “chip” 
delay exceeds 7 chips. 

9. Finally, the battery performance of the iPAQ is not appropriate for field equipment. 

We learned much about the 802.11b protocol. First, while the advertised speed is 11 
Megabits per second (Mbps), we could never achieve that as a channel throughput rate. In 
fact, it is documented [30] that, due to handshake overhead and synch headers, the maximum 
possible data throughput rate is only 7.1 Mbps for UDP and 5.9 Mbps for TCP links. We 
initially tested, using Iperf UDP, with a transfer rate of 7.35 Mbps (empirically derived) in 
the lab. This was used in the field to stress the link and determine the limits of performance, 
with distance being the only variable. After that test regimen, we lowered our rate to 5.5 


31 



3 METHODS, ASSUMPTIONS, AND PROCEDURES 


Mbps to have some headroom and get reliable performance out of the link. In our field 
testing we constantly ran up against the issue of low power versus distance. We needed 
to increase the TX-RX distance to simulate long multipath or delayed reflections. But, at 
the same time, the loss of RF power (r-squared free space loss) would sometimes put our 
signal too near the noise floor. The use of directional antennas was required. The possibility 
of purchasing an 802.11 RF power booster exists, but the expense was not justified. The 
Cushcraft S24012P panel antennas offer +12 dBi gain due to the pattern. One aspect of 
the 802.11 specification serves to meet an FCC requirement that limits the effective radiated 
power of 802.11b devices - the combination of RF power and antenna gain cannot exceed 4W 
EIRP [31]. One way the manufacturers hold to this rule is by using an antenna connector 
that makes it inconvenient to swap antennas on an Access Point. These “reverse gender” 
SMA connectors keep the AP/antenna combination proprietary. We either replaced the 
SMA on the AP or used a purchased adaptor to connect to the Cushcraft panels. We were 
still under the legal limit after the modification. 

We learned that two AP’s will not connect wirelessly; a “bridge” or “router” is needed 
for this function. We found that the antennas internal to a laptop might be in an unknown 
orientation. We did some limited experimentation to find an optimum laptop orientation, 
which wasn’t usually horizontal on a desktop! In the close proximity of an office environment, 
even the polarization loss from vertical to horizontal may be insignificant. For field use, 
commercial AP’s don’t offer enough RF output, and should be possible to power from DC 
directly; we had to carry bulky DC/AC inverters to run the measly “wall warts” for the 
AP’s. 

In order to use two directional antennas, we needed another standard SMA connection to 
our laptops. The Dell Intel wireless card was not amenable to having an external antenna. 
We tried a Linksys WUSB11 802/USB dongle, but it never achieved error free performance 
with the drivers provided. We borrowed a Linksys WET11 Bridge and connected via ethernet 
(cat5 cable) to the Dell (the Mac is cat5 to the D-link hardwire). Using the SMA gender 
changer, we could then have two wireless connection to the Cushcraft directional antennas 
(The YJ antenna is a standard SMA). In our experience, any change or addition to the 
Dell/Windows laptop was a lesson in frustration for Mac users. 

The results of our field testing confirmed most assumptions about the performance of 
802.11b. We had incremental surprises and detours. 

1. The performance of a commercial AP with its omnidirectional antenna falls off after 
about 60 meters in the free field. 

2. The presence of a reflector about 40m away does not destroy an established 802.11b 
Iperf link. 

3. The presence of a reflector about 90m away, along with the long path/low power 
situation does affect Iperf performance. 

4. The intentional orientation of the RX antenna in a multipath configuration (La Cueva) 
did result in a substantial degradation of Iperf performance. 
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4 Results and Discussion 

4.1 Introduction 

In this section, we describe the results from the validation study and the physical layer sim- 
ulations. The validation study is based outdoor measurements made in the region around 
NMSU. The simulation study is to estimate the characteristics of the physical layer perfor- 
mance for the 802.11 a/'b protocols on the Martian surface based on the DEM data. 

4.2 Validation of Power Delay Profile Simulations 

Power Delays Profiles were simulated at the Tortugas Mountain and Dripping Springs sites 
and a simple RF reflection measurement was also made at these sites. Table 4.1 shows 
the direct path and reflected path RF measurement results. For Tortugas Mountain, the 
reflection is 24 dBm less than the power for the direct path (since the sidelobe power levels 
were below the noise floor, the 24 dBm figure is entirely due to reflected power). 

Table 4.1: PDP validation data. 



Relative Power Level Difference 
Between Direct and Reflection Paths 

Reflection 

Delay 

Tortugas Mountain-Simulation 

20 - 25 dBm 

0.6 //s 

Tortugas Mountain-Measurement 

23 dBm 

1.4 //s 

Dripping Springs-Simulation 

26 dBm 

0.7 fj.s 

Dripping Springs-Measurement 

20.5 dBm 

0.87 /xs 


Using the appropriate DEMs, antenna patterns, radiated power levels, and positions of 
transmit and receive antennas, the power delay profiles were simulated for the local sites. 
Figure 4.1 gives the PDPs for Tortugas Mountain and Dripping Springs. In these figures, 
we see for Tortugas Mountain a strong response occurs at approximately 1.4 fxs and is 
25 dBm down from the direct path and for Dripping Springs a strong response occurs at 
approximately 0.7 /as and is 26 dBm down from the direct path. This agrees quite closely 
with measurement data. 

4.3 Power Delay Profile Measurement Data 

Figure 4.2 shows the PDPs measured using the YellowJacket at the NMSU football stadium 
parking lot (free field). With the access point at a fixed location, the PDP was measured at 
40 m and 80 m distances in the south and east directions (south leg and east leg). The plots 
show the PDP (relative correlation vs. delay in chips) at various points in time. Figure 4.3 
shows the PDPs measured at the NMSU football stadium parking lot in the presence of the 
reflector. With the access point at a fixed location, we measured the PDP at 40 m with and 
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(b) 


Figure 4.1: Power delay profile simulation results from (a) Tortugas Mountain and (b) 
Dripping Springs. 
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without the reflector. We clearly see in Figure 4.3(b) the secondary path delayed by 3 chip 
periods relative to the direct path centered at 0 chips. 

Figure 4.4 shows the PDPs (relative correlation vs. delay in chips) at various points in 
time measured using the YellowJacket at the La Cueva rock formations at Dripping Springs 
(complex multipath). Figure 4.4(a) illustrates the PDP when the antennas were pointed 
directly at each other. With this setup, there is essentially a single, direct path. Figure 
4.4(b) illustrates the PDP when the antennas were pointed at La Cueva. With this setup 
there is clearly, significant multipath with the assumed direct path centered at 0 chips and 
secondary reflections around 1.5 chips and —2 chips (aliased). 

4.4 Outdoor IEEE802.11b Results— Stadium Parking Lot 

Prior to measuring the performance of IEEE802.11b outdoors in a complex multipath en- 
vironment, measurements were made in the presence of a single reflector (three large steel 
bookshelves) in the NMSU football stadium parking lot. In the absence of the reflector the 
parking lot is essentially a free field (see previous subsection). As described in the previ- 
ous section, a baseline measurement was made without the bookshelf (directional antennas 
pointed toward each other) and a measurement was made with the reflector (mainlobes 
pointed at the bookshelf with first sidelobes pointed directly toward each other). The mea- 
surements were of network bandwidth (fastest data rate without significant UDP datagram 
loss for the baseline) and percentage of received packets. Figure 4.5 shows plots of the data 
for a 100 second time period. For the free field we see virtually no packet loss at a sustained 
rate of 5.5 Mbps, however, in the presence of a very simple reflector, we have significant 
packet error rates while attempting the same rate UDP data flow. 

4.5 Outdoor IEEE802.11b Results— Dripping Springs 

In order to access the performance of IEEE802.11b in the presence of multipath such as 
that found in an outdoor environment, several tests were conducted at the La Cueva rock 
formations at Dripping Springs. These tests were described in the previous section. For 
the direct path test, both directional antennas were pointed straight at each other and the 
resulting network bandwidth is given in Figure 4.6(a). This result serves as our baseline 
since there is no multipath involved. For the multipath test, each antenna was pointed at La 
Cueva which resulted in a strong multipath environment. The resulting network bandwidth 
is given in Figure 4.6(b). 

4.6 Simulated WLAN Performance at Different Sites 

Data packets for 802.11a or 802.11b were generated according to the PPDU format shown in 
Fig. 3.7. The simulation software is developed around the mWLAN toolbox from Comm Ac- 
cess Technologies [32], The data packets from the transmitter are sent through a random 
multipath channel generated for the particular transmitter and receiver locations on the 
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Figure 4.2: Power delay profiles of the NMSU football stadium parking lot (free field) as 
measured using the YellowJacket. Plots show the PDP (relative correlation vs. delay in 
chips) at various points in time, (a) South leg, 40 m distance; (b) South leg, 80 m distance; 
(c) East leg, 40 m distance; and (d) East leg, 80 m distance 
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Figure 4.3: Power delay profiles of the NMSU football stadium parking lot with and without 
a reflector present. Plots show the PDP (relative correlation vs. delay in chips) at various 
points in time, (a) without reflector showing a single, direct path centered at 0 chips and 
(b) with a reflector showing a secondary path at 3 chip periods and a direct path centered 
at 0 chips. 
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(a) 


(b) 


Figure 4.4: Power delay profiles of the La Cueva rock formations at Dripping Springs. Plots 
show the PDP (relative correlation vs. delay in chips) at various points in time, (a) with 
directional antennas pointed directly toward each other, direct path centered at 0 chips and 
(b) with directional antennas pointed toward rock formations, direct path centered at 0 chips 
and multiple secondary paths. 
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Figure 4.5: Outdoor IEEE802.11b network bandwidth and percent UDP received packets 
for (a) free field and (b) simple reflection. 


38 




4 RESULTS AND DISCUSSION 



Figure 4.6: Outdoor IEEE802.11b network bandwidth and percent UDP received packets 
for (a) free field and (b) multipath environment. 


Martian surface using the ICS Telecom software. The received packets are processed by the 
corresponding receiver. The 802.11b receiver’s performance is studied with and without a 
RAKE structure. Note that a RAKE receiver coherently combines different multipath con- 
tributions before detection and thus improves performance. For both 802.11a and b, only 
truncated channel impulse responses are estimated at the receiver using the corresponding 
PLCP Preamble. 

4.6.1 Performance versus distance between the transmitter and the receiver 

In order to obtain packet error rate (PER) and bit error rate (BER) results versus dis- 
tance, it is necessary to estimate both the received signal and the input-referred noise for 
an 802.11 receiver on the Martian surface. The RF propagation simulations using ICS Tele- 
com provide an estimate of electric field intensity at the receiving antenna. A first-order 
estimate of receiver noise is based on a noise figure Fr = 7.2 dB for a typical 802.11a re- 
ceiver implementation [33]. Assuming noise figure is measured for a reference temperature 
T 0 = 290 K, the equivalent noise temperature for the Martian receiver may be calculated [34] 
as T R = ( Fr — 1)T 0 = 1522 K. An omnidirectional antenna pattern sees roughly half sky and 
half surface, so we approximate the brightness temperature (A,) as T& = T v j 2 = 250 K / 2, 
where T p is the physical temperature. Further assuming a radiation efficiency rj = 0.9, we find 
an equivalent temperature for the antenna of Ta — r l'Fb + (1 — 7 /)^p = 138 K. Thus, our sim- 
ulations use an equivalent noise temperature for the receiver input of T eq = Ta + Tr = 1560 
K. 

The packet error rates (PER) for various distances (d) between the transmitter and 
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receiver are given in Tables 4.2-4.6. Note that a CRC failure is considered as a packet 
error in 802.11b while any error in the OFDM SIGNAL symbol constitutes a packet error 
in 802.11a. Transmit power is 1 W, and antenna height is 1.5 m above the ground, for both 
802.11a and b. The 802.11b results in the table are obtained without RAKE. The data rates 
for 802.11a and b axe 12 Mbps and 11 Mbps respectively. 

The packet error rate tables show that both 802.11a and b perform well for receivers 
within several hundred meters from the transmitter. In some cases, we find better packet 
error performance at a longer distance ( d ). For example, with 802.11a at Gusevl Site2, PER 
at 500 m appears to be better than the PER at 200 m. Similarly, in the case of Hematite4 
Sitel, the PER at 200 m is better than the PER at 100 m. In the case of Gusevl Site2, 
we observe that while the received power is higher at 200 m, the rms delay spread in this 
case is smaller for d = 500 m, resulting in fewer packet errors. A similar comment can be 
made about the Hematite4 Sitel PER result. In the case of 802.11b, the performance at 100 
m is better than the performance at 50 m for Gusevl Site2 as well as for Hematite4 Sitel. 
We notice a similar phenomenon as observed in the case of 802.11a, that is, although the 
received power is smaller for 100 m, the rms delay spread becomes smaller too. Thus, it 
appears that when sufficient power is transmitted (1 W in this case) the multipath effects 
play a dominant role on the performance of both 802.11a and 802.11b. Finally, we note that 
the results with very low PER values must be used with caution as they are not statistically 
significant due to the small number of packet errors observed from transmitting 20,000 data 
packets. 

The effect of distance on the bit error rate (BER) performance is shown in Fig. 4.7 for 
Gusevl Sitel. The BER result for each distance is an average over four different locations at 
the same distance. Transmitted power is 1 mW for all cases, and 802.11b results are obtained 
using a RAKE receiver. The data rate for 802.11a is 12 Mbps and the 802.11b transmits at 
the rate of 11 Mbps. In the case of 802.11a, the overall BER seems to increase with distance 
except for a strong dip at 500 m. This BER dip is believed to be due to favorable terrain 
conditions at that distance and it agrees well with the PER result in Table 4.2. The BER 
for 802.11b seems to be nearly constant up to distances of 1000 m except for a dip at 500 m 
similar to 802.11a. 

4.6.2 Effect of transmit power on PER 

Although Tables 4. 2-4. 6 show PER results for 1 W of transmit power, it is instructive 
to study the effects of transmit power on the PER. This is investigated via Table 4.7 for 
Gusevl Sitel. The table shows that when the transmit power is small, 802.11b seems to 
be doing better than 802.11a. As the transmit power increases, the performance for both 
802.11a and b tend to flatten out for high transmit power, but 802.11a shows much better 
performance than 802.11b. Note that the rms delay spread for this location is 0.105 ps 
for 802.11a, and it is much less than the available 0.8 ps guard period. Thus, 802.11a can 
handle this delay spread quite well, and its performance keeps improving with the transmit 
power. As the transmit power becomes large, however, the multipath events with delays 


40 


4 RESULTS AND DISCUSSION 


Table 4.2: Packet Error Rate Performance at Gusevl Sitel. A indicates zero packet errors, 

in 20,000 packets. 

d rms delay spread Received power 
(m) (fis) (nW) 


(fis) 


802.11a 802.11b 802.11a 



PER 

802.11a 

802.11b 

0.001)8 

0.0983 

0.0004 

0.0768 

0.0001 

0.0572 


0.0001 


0.0009 



0.0281 


0.0158 


0.9619 


Table 4.3: Packet Error Rate performance at Gusevl Site2. 


d 

(m) 

rms delay spread 

Received power 
(nW) 

PER 


802.11a 

802.11b 

802.11a 

802.11b 

802.11a 

802.11b 

20 

0.146 

0.186 

38.2 

79 

0.0003 

0.115 

50 

0.131 

0.155 

26.7 

56 

0.0004 

0.082 

100 

0.095 

0.126 

25.8 

54 

0.0001 

0.032 

200 

0.713 

0.719 

0.0822 

0.16 

0.099 

0.51 

500 

0.472 

0.476 

0.0114 

0.02 

0.067 

0.53 


Table 4.4: Packet Error Rate performance at Gusevl Site3. A indicates zero packet errors, 

in 20,000 pac kets. 

d rms delay spread 
(m) (n s) 

802.11a I 802.11b 


Received power 
(nW) 

802.11a I 802.11b 



802.11a 




802.11b 


11.2 

9.2 

.00001 

0.0001 

0.742 

0.718 

.000001 

.00003 
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Table 4.5: Packet Error Rate Performance at Hematite4 Sitel. 


d 

(m) 

rms delay spread 

(^s) 

Received power 
(nW) 

PER 

802.11a 

802.11b 

802.11a 

802.11b 

802.11a 

802.11b 

20 

0.741 

0.634 

59.2 

114.27 

0.0262 

0.2113 

50 

0.747 

0.625 

49.1 

94.02 

0.0272 

0.2844 

100 

0.584 

0.564 

47.0 

80.62 

0.0138 

0.1667 

200 

0.289 

0.297 

29.3 

46.07 

0.0026 

0.1196 

500 

0.069 

0.087 


43.24 



1000 

0.696 


.0374 

0.167 




Pable 4 

.6: Packet Error R 


d 

(m) 

rms delay spread 
(**) 

Received power 
(nW) 

PER 

802.11a 

802.11b 



802.11a 

802.11b 

20 

1.031 

0.913 

45.9 

| 

0.0037 

0.1280 

50 

0.755 

0.694 

36.2 


0.0012 

0.0950 

BH 


0.498 

29.2 

55.97 

0.0004 

0.0724 


0.178 

0.228 

28.6 

55.30 

0.0001 

0.0354 

500 



42.3 

88.95 

0.0004 

0.0370 

1000 

0.287 


2.10" 7 

2 . 10 — 6 

1 

1 


Table 4.7: Effect of Transmit Power on PER for Gusevl Sitel at a distance of 100 m from 


Transmit Power 

PER for 802.11a 


1 /rW 

0.985 


10 /rW 

0.381 


100 /iW 

0.0225 


1 mW 

0.0021 


10 mW 

4 x 10~ 4 


100 mW 

2.5 x 10- 4 

0.0555 

1 W 


0.0516 
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Figure 4.7: BER Performance for 802.11a and b at Gusevl Sitel. 

exceeding 0.8 /ns start affecting its performance with adjacent symbol interference. This 
limits the performance improvement. In the case of 802.11b too, the presence of multipath 
propagation effects does not allow performance improvement beyond a certain value. 

4.6.3 BER Performance versus SNR 

The bit error rate (BER) performance results versus SNR are shown in Figs. 4.8-4.17 for 
IEEE802.11a and b. In the case of 802.11a, we notice that lower data rates provide much 
better BER performance giving several dB advantage over higher rates. However, it is 
also to be noted that low r er rates need to transmit longer than higher rate modes in order 
to send the same amount of information. We also see that the curves tend to flatten at 
the higher SNR region as the performance becomes more dominated by the delay spreads. 
Although the rms delay spread is within 0.8 /xs for the cases studied in these figures, there 
are still multipath channels beyond 0.8 /xs producing adjacent symbol interference. The 
BER performance curves for 802.11b show that multipath channels can severely affect their 
performance. Figures 4.13-4.17 show results without a RAKE structure. Another interesting 
observation is that CCK performs better than the other modulations in some cases. 

4.6.4 Effect of using RAKE for 802.11b 

The use of a RAKE receiver can significantly improve the BER and PER performance for 
802.11b. The BER performance improvements can be seen comparing Figs. 4.13 and 4.18 
for Gusevl Sitel. The PER performance improvements are summarized in Table 4.8. The 
table shows that RAKE can provide PER improvement by a factor as high as eight in this 
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Figure 4.11: BER Performance for 802.11a at Hematite4 Sitel 
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Figure 4.12: BER Performance for 802.11a at Hematite5 Sitel. 



Figure 4.13: BER Performance for 802.11b at Gusevl Sitel without a RAKE structure. 
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Figure 4.14: BER Performance for 802.11b at Gusevl Site2 without a RAKE structure. 



Figure 4.15: BER Performance for 802.11b at Gusevl Site3 without a RAKE structure. 
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Figure 4.18: BER Performance for 802.11b at Gusevl Sitel with a RAKE receiver. 

case. The performance improvement seems to be generally smaller at very large distances. 

4.6.5 Performance versus antenna heights 

The antenna heights can affect the performance of both 802.11a and b significantly. An 
increase in the antenna heights can provide better line-of-sight signals over a rocky terrain 
and can increase the received power. However, it can result in more delay spreads as well, 
resulting in decreased performance at the receiver. In the case of 802.11a, we can observe 
from the PER tables that, of the three sites considered, Gusevl Site3 has the least rms delay 
spread at 100 m. Since the received power is too low, the benefit from an increase in the 
received power becomes significant since the rms delay spread remains much smaller than 
the guard interval. Thus, when the antenna heights are raised, this site shows significant 
improvement in performance despite an increase in the rms delay spread value from 0.036 pts 
(corresponding to antenna height of 0.5 m) to 0.058 p. s (corresponding to antenna height of 
2.0 m). The results in 802.11b do not show significant improvements with antenna heights as 
in 802.11a. This may be because the benefit due to more received power is nearly cancelled 
by the loss due to increased delay spreads. Finally, increasing the heights of the antennas 
beyond a certain value may be impractical for mobile rovers. 

4.6.6 Discussion 

There are several interesting observations. 

• The received power for 802.11b is always greater than 802.11a. This makes sense since 
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Table 4.8: Packet Error Rate Performance at three sites in Gusevl for IEEE 802.11b. The 
indicates non-availability of results. The transmit power is 1 W and the antenna heights 
are fixed at 1.5 m above the ground. 


d 

(m) 

Site 1 

Site 2 

Site 3 

Without 

RAKE 

With 

RAKE 

Without 

RAKE 

With 

RAKE 

Without 

RAKE 

With 

RAKE 

20 

0.098 

0.024 

0.115 

0.022 

0.10 

0.043 

50 

0.077 

0.014 

0.082 

0.011 

0.027 

0.008 

100 

0.057 

0.01 

0.032 

0.008 

0.016 

0.005 

200 

0.028 

0.005 

0.51 

0.284 

0.03 

0.006 

500 

0.016 

0.002 

0.53 

0.315 

0.54 

0.234 

1000 

0.962 

0.52 

- 

- 

1.00 

0.682 



Figure 4.19: BER Performance for 802.11a at Gusevl Sitel. The transmit power is 100 /xW, 
and the distance ( d ) between the transmitter and the receiver is 100 m. 
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Figure 4.20: BER Performance for 802.11b at Gusevl Sitel using a RAKE receiver. The 
transmit power is 100 /xW, and the distance (d) between the transmitter and the receiver is 
100 m. 

the transmit frequency for 802.11a is in the 5 GHz band while the transmit frequency 
for 802.11b is in the 2.4 GHz band. 

• For shorter distances, the rms delay spread for 802.11a seems to be smaller than for 
802.11b in the Gusev sites considered. For larger distances, the rms delay spread 
for 802.11a increases and becomes similar to or larger than the rms delay spread for 
802.11b. The behavior seems to be just the opposite at the Hematite sites. 

• The performance of 802.11a and b is affected by received power and the existence 
of multipath propagation conditions. When the received power is too small, we can 
say that the system is operating in the power limited region. An increase in power 
in the power constrained region improves the performance. On the other hand, when 
sufficient power is received, the performance of the system can still be severely degraded 
due to multipath effects. In this case, we can say that the system is in the multipath 
limited (or equivalently bandwidth limited) region. In the multipath limited region, 
the performance of the system does not improve with transmission of any additional 
power. 

• The PER for 802.11a is observed to be smaller than the PER for 802.11b in almost 
all cases for large transmit power (1 W). For smaller transmit power, 802.11b seems 
to perform better than 802.11a. Note that the received power for 802.11b is larger 
than for 802.11a as they use 2.4 GHz and 5 GHz frequency bands respectively. This 
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higher received power greatly helps 802.11b in this power limited region. A meaningful 
comparison between the two, however, should include the effects of packet sizes, over- 
heads, possible improvement due to RAKE in 802.11b, and implementation complexity 
considerations as well. 
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5 Conclusions 

We have continued the effort begun under Year 1 of the grant to begin applying the technol- 
ogy developed for propagation modelling on the surface of Mars using commercial software 
products and available map data sets. The effort in Year 2 had two major phases, vali- 
dation of simulation techniques and IEEE 802.11a/b physical layer modelling. The major 
accomplishments for the second year of this research are as follows: 

1. Accomplishment 1. Using the Year 1 results regarding antenna coverage patterns, 
maximum link distances, effects of surface clutter, and multipath effects, we have 
simulated the physical layers of IEEE 802.11a and IEEE 802.11b wireless networking 
standards in the Martian environment. The simulations were conducted in MATLAB 
using the CommAccess Technologies’ mWLAN toolbox. These results will be used 
as the basis to begin the study of the behavior of the medium access layer for these 
protocols in the simulated Martian environment. 

2. Accomplishment 2. We performed a detailed validation study of the results of the 
computer-based techniques for estimating the physical layer propagation effects against 
those measured in the field. In this validation study, we measured key parameters (RF 
signal strength, delay spread, data rates, and packet error rates) for an outdoor IEEE 
802.11b wireless network. These measurements were conducted at sites near Tortugas 
Mountain and Dripping Springs due their similarity to the Mars surface (free of man- 
made objects, little vegetation, mostly flat with some terrain variation and rocks, 
etc.). These measurements were compared with expected performance based upon the 
DEMs for the sites. We judge the agreement to be very good, especially considering the 
wide variation in measurements that can occur by moving the transmitter or receiver 
antenna by a small distance. 

From this study effort, we conclude that: 

• The results of the physical layer simulations for the Martian surface show that success- 
ful IEEE 802.11a/b-based communications are possible within a few hundred meters 
of the transmit antenna when the transmit power is more than a few milliwatts and 
the antenna heights are fixed at more than 1 meter above the ground. The packet er- 
ror rate performance of 802.11b without a RAKE receiver seems to be more adversely 
affected by the multipath conditions than 802.11a. Further, the lowest data rate mode 
of 802.11a provides the best bit error performance. 

• Using higher power in the communications system does not always help the perfor- 
mance of the system. This result is known for terrestrial environments where vegeta- 
tion and atmospheric phenomena are important. It is also true in sparse environments 
with negligible atmospheric attenuation as well. 

• Transmission power and antenna height can be traded to a certain extent. The desired 
link Quality of Service and data rate may be more of a driver in link design and 
performance than transmission power and antenna height. 
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• Use of a RAKE-type of receiver can significantly improve performance with 802.11b 
protocols. 
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6 Recommendations 

Given our research results, we propose the following recommendations: 

1. Continue investigation of the technology to employ the MATLAB physical layer sim- 
ulations with the OPNET MAC layer simulations. The incorporation of MATLAB 
with OPNET to produce a cosimulation is proving to be a difficult but necessary task. 
However, recent efforts have shown that there appears to be a path to resolving the 
compatibility issues. Once these issues are resolved, then further progress in seeing the 
effects of the physical channel on MAC layer performance will follow. 

2. Investigate MAC layer improvements to make outdoor wireless communications more 
reliable. The 802.11 protocol suite has a rich set of features that can be investigated for 
how they can be improved to make the communications more reliable. Most of these 
features were initially designed for indoor wireless use and not the outdoor environment. 
Now that we can predict the propagation effects, both link budget and multipath, with 
a degree of reliability, we have the means to investigate which protocol features are 
most important for the 802.11 wireless environment. 

3. Investigate the potential for other protocol families to be used in the planetary surface 
environment. For example, it is expected that the European Space Agency will con- 
sider using the 802.16/large subcarrier number OFDM as the baseline protocol surface 
communications in their mission studies. With the technology baseline being developed 
here, we can assess if this protocol suite would make sense when the various propa- 
gation factors are considered. The performance can then be benchmarked against the 
802.11 simulations for evaluation purposes. 
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C Steps for Generating Simulations 

C.l RF Coverage Simulations using HerTZ Mapper 

Initial Setup 

To create a project in HerTZ Mapper you must have three files available 

1. .geo: This file contains the topographic information 

2. .img: This is a image file corresponding to the topographic map 

3. .pal: This file contains the palette information of the image 

Creating a Project 

From the Windows Start menu , navigate to the ATDI software selection within Programs, 
and choose HerTZ Mapper. Once the package has opened, click on the image on the middle 
of the screen. From the File menu, you may now choose Create New File Server to create 
a new project, or Coverage and profile analysis to open an existing one. The procedure is 
explained below. The File server window is shown in Fig. C.l. 

• Enter a file name for project or select an existing.svr file 

• Double-click on Digital Elevation Model, then on Add, and select a .geo file 

• Proceed the same way for the Primary Image, and add first the .img file associated to 
the DEM (then the .pal file) 

• “click Open Coverage” & Profile, and choose OK 


Transmitter setup 

Once you have the files opened in HerTZ Mapper you can click on the DEM or image and 
select “Add site” to place a transmitter. You can also specify the latitude and longitude 
of the location to exactly place a Tx at desired location. When you select Add site a Tx 
parameter window is displayed as one shown in Fig. C.2. The user has to enter all the Tx 
parameters such as frequency, power and antenna height. Once you choose OK a red square 
is displayed on DEM and image with a letter on it as shown in Fig. C.3. This red square 
shows the location of transmitter. 
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Figure C.l: HerTZ Mapper file server window. 
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Figure C.2: HerTZ Mapper station parameters. 



Figure C.3: HerTZ Mapper opening window. 
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Figure C.4: HerTZ Mapper model selection window. 



Figure C.5: HerTZ Mapper ITM parameter window. 
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Propagation Model Setup 

To select a propagation model select Parameters form “File Menu” then “Engineering par 
rameters” and then “Coverage”. Once you click on the coverage a window is selected with 
all the propagation models as shown in Fig. C.4. First you have to select “External dll” 
and then “ITM122”. Once you click OK the ITM parameter window is displayed as shown 
in Fig. C.5. Once you enter all the required values you can close the window to finish the 
propagation model setup. 



Figure C.6: HerTZ Mapper opening window. 


Coverage Calculation and Overlay 

First click on the “Transmitter” and then select the Transmitter and then select “coverage” 
to r un a coverage. Clicking on the Coverage option will instruct the software to simulate 
the site you clicked on, using the parameters and propagation model you have previously set 
up. The computer performs a simulation. Once finished, you are prompted to save the field 
strength file. Click “Save” for this file box (it automatically suggests a name). Minimize or 
close the field strength result window and click on the map. The field strength file can now 
be displayed in full transparency on top of the image by choosing the “Overlay FLD” option 
(in the popup box), and selecting the appropriate field file. The color of the field denotes 
a discrete band of values of field strength. To see the actual value of field, move the mouse 
cursor over the appropriate area and read the value from the mformation bar at the bottom 
of the screen. It appears in the dBV/m bracket. You can modify the colors and field levels 
in the Parameters/Legend box. Fig. C.6 shows the coverage overlayed on the DEM. 
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C.2 Power Delay Profile Simulations using ICS Telecom 

Initial setup 

As mentioned for HerTZ Mapper the ICS Telecom also needs three files (.geo,.img, and .pal) 
to create a project. 

Creating a Project 

From the Windows Start menu, navigate to the ATDI software selection within Programs, 
and choose ICS Telecom. Once the package has opened, From the File menu choose project 
Manager to create a new project, or to open an existing one. 

Once the project Manager window opened click on the 3 dots in front of the DEM, image, 
and pal options to select the .geo, .img, and .pal file locations. Once all 3 files are selected 
click save as to create the project file with extension .PJT. 

Transmitter setup 

Select the point mode from the tool box. Now click on the DEM and then select add Tx/Rx 
from the popup menu. A parameter window will be displayed. Enter the Tx parameters 
leaving out the options not applicable for the situation. Once you click OK a square with a 
letter on it is displayed on DEM to represent the transmitter position. 

Propagation Model setup 

Select File menu Coverage/Network Calculation/Tx/Rx FS calculation. A window is dis- 
played with Rx antenna height field. Now click on more to get a Advanced coverage pa- 
rameters window. Click on Model to get Propagation models window. The following three 
selection have to be made in this window. 

In Models box click on 3 dots to select the ITM (tuned dll) from the Models directory 
within ICS. In climate box enter the corresponding Mars land and sea radius in Km and 
then click on the 3 dots to select the k factor as 1. In Reflection box enter the reflectance 
value and then select 3d coverage only option. 

PDP calculation 

1. Set on the terrain the transmitter network up(Add Tx/Rx option of the point mode 
left click popup menu). 

2. Proceed to the coverage calculation of the network(File menu Coverage/Network cal- 
culation/Tx/Rx Fs coverage). 

3. Display the composite coverage(File Menu Coverage/Network Analysis/Composite cov- 
erage display) and keep the coverage overlaid on the terrain. Keep the Point mode 
activated. 
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4. Click over the site of the Virtual receiver In the popup menu, choose option Multipath. 

5. In the input box, enter the icon number of the transmitter and the receiver ante nna 
heighten meters). 

The simulation is run and will calculate the reflected field strength and the ToA emanat- 
ing from the designated transmitter to the receiver point touched by the cursor. The result 
is shown as a graph displayed in the Multipath box. The blue color represents the reflected 
ray, and the red color indicates the direct ray. In the status bar, the following information 
regarding the position of the cursor on the graph, is given: 

ToA value in /is 

Reception threshold(FSR) in dB/xV/m 

Altitude in meters 

The print button opens the Print setup box to print the graph. The close button closes 
the Multipath box. 
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D Visualization of Data Acquired with the Yellow- 
Jacket 

D.l Log File Transfer 

The YellowJacket will save measurement data on the iPAQ in a single file named with 
the .YJ3 extension. The default name is log.yj3. In order to save multiple sessions or 
measurements, be sure to rename the file and start again, as the YJ will simply append to 
the existing log file. 

The only way to transfer data out of the YJ is to use the infrared port on the top of the 
iPAQ. This must connect to a corresponding infrared port on a Windows-based PC. The 
software used for the transfer is Microsoft Active Sync. 

1. On the PC, start Microsoft Active Sync, then point the YJ at the ir port of the PC. 
On the YJ, use the main menu screen to select System, then Communications, then 
Connect via IR. 


rz 

Microsoft ActiveSync 

B 

I File View loois Help 


* . • ' Details • • . - 

: .. , r? c ^ ,, 




Not connected 


Information Type Status 

— — - — « 


Figure D.7: Snapshot of ActiveSync. 
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2. When the un i ts “connect” , Active Sync will make a noise and ask what type of part- 
nership is desired. This should be a guest connection, not a synched connection. 

3. While Active Sync should be able to allow the user to browse the files on the remote 
unit, it usually crashes. Try the normal “my computer” method to locate the remote 
unit and look for the log.yj3 files in the home folder. Copy the file(s) to the desktop 
or appropriate folder. 

4. Now that the file has been copied (with translations from PDA to desktop format) to 
the PC, it must be translated from YJ format to columnar. Use Berkeley Varitronics 
Systems (BVS) Chameleon software for this. 



Figure D.8: Snapshot of Chameleon. 


5. In Chameleon, follow the steps numbered to the left of the Chameleon screen: 

(a) Browse for the input file. The default screen will search for the .yj3 extension 
(extensions axe not always visible in Windows machines) and the file type: BVS 
YellowJacket III Log Files. 

(b) Browse for the output file (create a file name) with the default extension .out 

(c) Select Output Format: MS Excel; select Field Delimiter: Tab 

(d) Choose the data source from the “Field Selection Available.” The YJ stores the 
data from the screens in these categories. The most useful are the AP/Sta Data 
(Access Point station data) and GPS. 
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(e) Once a source is selected (by placing the selection is the Selected box using the 
arrows), the DATA: Fields become available. Here, the individual selections are 
made and will include all sequential samples of the various measurements. Al- 
ways select a Channel and/or SSID to be sure which station is being measured. 
The Multipath Data are the first twenty two measurements from the Correlation 
screen, from which the Delay Spread can be derived. 

(f) Use Microsoft Excel to view the Chameleon XXX.out. file and put the data in 
spreadsheet format. The Excel wdzard will assist in translating the layout. The 
file can be saved as XXX. xLs and opened in MATLAB (see next section regarding 
out Yellowjacket Toolbox User’s Guide). For viewing in Excel, some adjustment 
of column width will be required. 



Figure D.9: Snapshot of Excel. 


D.2 MATLAB Yellowjacket Toolbox User’s Guide 

Open MATLAB by double clicking on the MATLAB icon as in Fig. D.10. 
As in Fig. figxommand, type at command prompt 

» YJ_program 
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Figure D.10: Snapshot of Windows desktop. 

Select YellowJacket data folder by date and antenna height. Then choose which Excel 
file by distance from AP. See Fig. D.12. 

Now the program is running as in Fig. D.13. NOTE: To view Iperf data for location 
1. Select file Iperf 
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Figure D.ll: MATLAB’s command window. 
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Figure D.12: YellowJacket toolbox program running. 
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Figure D.13: YellowJacket toolbox program r unning . 
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E IEEE 802 . 11 b Field Measurement Test Procedure 

E.l Introduction - Free Field and Simple Reflector Tests 

Purpose: To verify IEEE 802.11b WLAN hardware and RF performance. We intend to 
discover the free field distance that can support data transfer between commercial access 
points (AP) and standard laptop computer transceivers (built in). In addition, a simple 
RADAR reflection configuration will be measured to determine the effect of multipath signals 
on the performance of an 802.11b system. 

E.2 Equipment 

D-Link DI-524 wireless router / access point (AP) 

SSID = “spirit” or “opportunity” 

Apple Macintosh PowerBook G4 laptop 
SSID = “isidis” or “dsp” 

Dell Latitude D600 laptop 
SSID = “gusev” 


Yellow Jacket Plus (with GPS) - also called “YJ” 

802.11b W-LAN Analysis System 
Berkeley Varitronics Systems 
s/n 024127 

Cushcraft S24012P directional panel antenna(s) 

Garmin rinol20 FRS radios / GPS 
Garmin eTrex GPS receiver 

APC SmartUPS 700 

uninterruptible power supply (battery powered AC source) 
and / or 

Xantrex xPower 600 
DC- AC power source 

software: 

“Iperf” network testing utility 
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E.3 Procedure 

Setup 

1. Connect equipment to power source as required: 
laptop “isidis” 

D-Link AP 

YJ charger cradle, if needed 

2. Establish appropriate connections to AP: 
ethernet cable between Iperf “client” and AP 
RF link between Iperf “server” and AP 

open “SSH Secure Shell Client” window on server to remotely control client: password 
as listed 

run Iperf as client, -c 

3. Verify YJ is acquiring AP by channel and SSID 
under YJ options, clear log file and then enable logging 

Iperf data test 

1. Setup Iperf server and client UDP command lines: 

(server) iperf -s -u -i 1 xxxlog.txt 

where xxx is log file descriptor (distance, elevation) such as “s40-lm” 

(client) iperf -c lO.O.O.xx -u -b 7.35m -t 120 

where .xx is the wireless IP address of the server computer 

determine clean data rate and substitute for 7.35m in above command line 

2. While monitoring distance on the GPS, move away from the AP, watching the data 
quality on Iperf. 

3. Take data at 20 meter increments to limit of system function and data quality. See 
Figure E.14 for positions. 

4. Rename the log file on the YJ, labelling each measurement station 


77 


E IEEE 802. 1 IB FIELD MEASUREMENT TEST PROCEDURE 



Reflection test 

Establish the geometric layout depicted in Figure E.15. 

1. Verily wireless direct link between server and client. 

2. With YJ at far corner of layout triangle, look for double correlation peak by manipu- 
lating orientation or antenna angles. 

3. Name the log file on the YJ 

4. With double peak showing on YJ, attempt Iperf data run and log results. 

E.4 Introduction - Multipath Tests 

Purpose: In a terrestrial environment that features massive rock structures capable of causing 
multipath reflections, determine the signal characteristics and data link effects due to the 
terrain. 


E.5 Equipment 

D-Link model DI-524 wireless router / access point (AP) 
SSID = “spirit” or “opportunity” 
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Reflector 



Figure E.15: Outdoor geometry. 
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Linksys Bridge model WET 11 
SSID = “opportunity” (WLAN AP) 

Apple Macintosh PowerBook G4 laptop 
SSID = “isidis” 

Dell Latitude D600 laptop 
SSID = “gusev” 


YellowJacket Plus (with GPS) - also called “YJ” 

802.11b W-LAN Analysis System 
Berkeley Varitronics Systems 
s/n 024127 

Cushcraft S24012P directional panel antenna(s) 

Garmin rinol20 FRS radios / GPS 
Garmin eTrex GPS receiver 

APC SmartUPS 700 

uninterruptible power supply (battery powered AC source) 
Xantrex xPower 600 
DC-AC power source 

software: 

“Ip erf” network testing utility 
“Ping” network connection utility 

E.6 Procedure 

Setup 

1. Connect equipment to power source as required: 
laptop “isidis” 

D-Link AP 

YJ charger cradle, if needed 

2. Verify YJ is acquiring AP by channel and SSID 
under YJ options, clear log file and then enable logging 
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Wireless data test 

1. Establish a wireless connection with omnidirectional antennas, using Ping as an exam- 
ple data test (see Figure E.16). 

2. Raise antenna to two meter height and evaluate data improvement, if any. 

3. For sufficient signal using 50mW AP, connect directional antenna(s) and retest. 

4. Find optimal orientation to establish multipath peaks on YJ display, using rock face 
as reflector. Mark RX spot as GPS waypoint. 

5. Name log file on YJ. 

^ ^ ,r ' Network Utility 

—{ Info Netstat AppleTalk Ping Lookup Traceroute Whois Finger Port Scan }- 


Please enter the network address to ping 


10.0.0.10 (ex. 10.0.2.1 or www.dofnain.com) 


C Send an u nlimited num ber of pings 

® Send only S pings ( pj n g j 


Ping has started ... 

PING 10.0.0.10 (10.0.0.10): 56 data bytes 
64 bytes from 1 0.0.0.1 0: icmp_seq=0 ttl=1 28 time=2.082 ms 
64 bytes from 10.0.0.10: icmp_seq=1 ttl=128 time=1 .983 ms 
64 bytes from 10.0.0.10: icmp_seq=2 ttl=128 time =2. 094 ms 
64 bytes from 1 0.0.0.1 0: icmp_seq=3 ttl=1 28 time=1 .861 ms 
64 bytes from 10.0.0.10:icmp_seq=4tll=128time=1.921 ms 

— 10.0.0.10 ping statistics — 

5 packets transmitted, 5 packets received, 0% packet loss 
round-trip min/avg/max * 1 .861/1 .988)2.094 ms 


/, 


Figure E.16: Sample output from ping test. 


Iperf Multipath test 

1. Setup UPS power to run laptop and Bridge at RX location. 
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2. Establish appropriate connections to AP: 

ethernet cable between Iperf “client” and AP at TX location. 

RF link between Iperf “server” and AP using Linksys Bridge as intermediate element 
(to enable use of directional antenna) at RX location. 

3. Carefully, without moving RX antenna position, connect antenna to Bridge in place of 
Yellow J acket . 

4. Rim Ping to verify wireless link. 

5. Setup Iperf server and client UDP command lines: 

(server) iperf -s -u -i 1 xxxlog.txt 

where vw is loc file deserintnr 

■ ~ — - - o r * 

(client) iperf -c 20.0.0.xx -u -b 5.5m -t 60 

where .xx is the wireless IP address of the server computer 

6. Perform sufficient Iperf runs to establish link behavior under multipath conditions. 
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F Code Listings 

This Appendix lists the major MATLAB code files used in the analysis and simulation work. 

Yellow Jacket Code 

ask_npts.m 

1 function answer = ask_npts (handles) 

2 

3 default = num2str (length (handles. current .data)*. 1) ; 

4 prompt = {[’How many points (’ num2str (length (handles . current _data) ) ’):’])•; 

5 dlg_title = ’Moving Average filter’; 

6 miin_lines= 1; 

7 def = {default} 

8 answer = str2num(char(inputdlg (prompt, ... 

9 dlg_title,num_lines,def))) ; 

10 if isempty(answer) 

11 answer = . l*length (handles. current.data) ; 

12 end 


BW_filter.m 

1 function BW = BW_f ilter(BWO) 

2 

3 for i = 1: length (BWO) 

4 

5 if BWO(i) < 7.5 

6 BW(i) = BWO(i); 

7 end 

8 end 


cellJind.m 

1 function [labels, index, in_ref] = cell_find(s, char_array) 

2 ‘/.This program find a string in a cell array and return the index 

3 

4 

5 •/,================================ 

6 %This is not correct 

7 ‘/.The Data Rate has a extra tab char ac tar 
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9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 


[R C] = size(s); '/.define search criteria 

[Rc Cc] = size(char_array) ; 

index = 0 ; '/.Set intial index 

in_ref = [] ; 

ind_end = 1; 

for j = l:Cc 


for i = 1:C 


x = double (char (s{l, i}) ) ; 
if x(l) == 9 

labels{i> = {char(x(2:length(x)))>; 

else 

labels{i} = s{l,i}; 

end 


'/.If desired string of given headers 

if length(char(labels-(i})) == length (char(char_array(j))> ... 

& char(labels{i» == char(char_array(j)) 

’/.If consecutive headers are the same 

if i“=l & length(char(labels{i}) ) == length (char (labels{i-l}) ) ... 

& char(labels-(i}) == char (labels{i-l}) 
if ind_end = 1; ’/(describe when consecutive headers began 
ind_end = 0; 
initial = i— 1 ; 
len = length (in_ref ) ; 

end 

else 

in_ref {length (in_ref )+l} = [{i} labels{i}] ; 

end 

index(length(index)+l) = i; 

elseif i >3 & length(char(labels{i-l» ) == length(char(labels{i-2») ... 
& char (labels{i-l}) == char(labels{i-2}) 
if ind_end == 0 

in_ref{len> = [{initial} {i-1} labels{i-l>] ; 
ind_end = 1; 

end 

end 
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54 end'/ochar^array search 

55 end%header search 

56 

57 


cellarray.m 

1 function names = cellarray(ssid) 

2 

3 len - length(ssid) ; 

4 

5 for i = l:len 

6 names{i,l> = char (ssid{iHl>{l>) ; 

7 end 


colJincLm 

1 function handles = col_find (handles) 

2 ^Called from YJ_workl on line 96 

3 Y.This program read the colums specified in handles. ref 

4 

5 row = [2 length (handles.dat a) ]; /(Length of data 

6 /(Columns 

7 7.SSID 

8 temp = handles .ref 

9 handles . col_ssid = temp; 

10 /(Multipath start 

11 temp = handles. ref {2} ( 1) ; 

12 handles . col_in = temp{l}(l); 

13 ‘/.Multipath finish 

14 temp = handles. ref {2} (2) ; 

15 handles . col_fn = temp{l}(l); 

16 /(Data Rate 

17 temp = handles. ref 

18 handles. DR = data_f etch (row, temp, handle s . data) ; 

19 /(RSSI 

20 temp = handles .ref 

21 handles .rssi = data_f etch (row, temp, handle s . data) ; 

22 

23 

24 if length (handles .ref ) > 4 
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25 ‘/(Latitude 

26 temp = handles. ref {5}{1} ; 

27 handles .latitude = data.fetch (row, temp, handles .data) ; 

28 end 

29 if length (handles. ref ) > 5 

30 /(Longitude 

31 temp = handles. ref {6}{1}; 

32 handles . longitude = aata_f etch (row, temp, handles . data) ; 

33 end 

34 if length (handles. ref ) > 6 

35 '/.Altitude 

36 temp - handles .ref {7}{1}; 

37 handles. Altitude - data„f etch (row, temp, handles .data) ; 

38 end 


dataJetch.m 

function d = data_f etch (row, col, data) 

[r c] = size (col); 

/(This is for Multiple col umn data 
if c >1 


for i = row(l) :row(2) 

for j = col(l) :col(2) 

d(i-row(l)+l , j) = str2num(char(data{i, j») ; 


10 

end 

11 

end 

12 

else '/.th: 

13 

if isem] 

14 


15 

for 

16 


17 


18 

end 

19 

else 

20 

for 

21 


22 


23 

end 

24 

end 

25 

end 

26 



i = row(l) :row(2) 
y - char (data{i, col}) ; 
d(i*-row(l)+l) = str2num(y(l:5)) ; 


i = row(l) :row(2) 
y = char (datafi, col}) ; 
d(i-row(l)+l) = str2num(y) ; 
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27 


delay .spread . m 

1 function d = delay_spread(MP) 

2 

3 

4 N = 22; 

5 T = 1/22000000; 

6 n=[0:N-l] ; 

7 [numrows.numcols] = size(MP); 

8 d = zeros (numrows , 1) ; 

9 

10 for i = l:numrows 

11 m = sum(MP(i,3:24).*n)/sum(MP(i,3:24)); 

12 d(i) = T*sqrt(sum(MP(i, 3:24) .* ((n-m) . ~2))/sum(MP(i, 3:24))) ; 

13 end 

14 d = d*le9; 


edit.external.plot.m 

1 load handles 

2 

3 PD = handles. MP(: ,4:25); 

4 [r c] = size (PD) 

5 start = 1; 

6 stop = r; 

7 figure (23); 

8 

9 V, Specify Range for plot 

10 % start = 2000 

11 % stop = 2500 

12 

13 

14 

15 ’/. Plot Power Delay Profile 

16 X = -2:5/c:3-(4/c) ; 

17 Y = start: stop; 

18 mesh(X,Y,PD(Y, :)) 

19 ylabel(’ Relative Time’); 

20 xlabel (’ Chips’) ; 
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21 zlabel( ’Relative Correlation* ) ; 

22 title(*When pointed at La CuevaO; 


for_data.m 

1 function lab = for_data(s) 

2 

3 ^Format data 

4 7. clear; clc; close all 

5 7. 7.file name of data 

6 7o filename = , log629_ap_mtlb. out * ; 

7 7. disp (filename) 

8 7. s = importdataCf ilename, *\t *) ; 

9 7. disp ( * *) 

10 

11 clc; 

12 rows = length(s); 

13 7.tab = s{1}(1,8); 

14 load tab_char 7.tab character 

15 for j = l:rows 

16 

17 index = f ind(char (s{j» == tab); 7.find tab character 

18 

19 len = length ( index ) ; 

20 

21 temp= s{j>(l : index(l)-l) ; 7oSave data string (initial) 

22 lab{j,l> - cellstr(temp(l , : )) ; 7.1oad data string (initial) 

23 

24 for i = l:len-l 

25 temp = s{j > (index (i) +1 : index (i+l)-l) ; 

26 lab{j,i+l> - cellstr (temp) ; 

27 end 

28 lab{j,i+2> = cellstr (s{j>(index(i+l) : length (s{j>))) ; 7.final 

29 end 

30 

31 

32 


generalJperf.m 

1 function iperf = general^iperf (file, path) 
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2 7 . 

3 7. SPACE DELIMITER 

4 % 

5 % USE get_iperf_data.m if possible 

6 7 .==========— ====-================================ 

7 X The program reads data files taken by iperf 

8 % and then literally copied to an Excel spreadsheat 

9 7. Process : 

10 % 1) open text file of iperf data 

11 7. 2) Ctl+a 

12 7. 3) Ctl+c 

13 7. 4) In Excel spreadsheet ctl+v 

14 7. 5) File Save 

15 7. Now run this program and the iperf data will be stored as followed 

16 % all in the iperf handle, iperf.: 

17 7. Time in seconds (iperf. time) 

18 7. Transfer in KBytes (iperf .transfer) 

19 7* Bandwidth in Mbit/sec (iperf .BW) 

20 % Packet Error Rate in percent (iperf .PER 

21 7 , 

22 7 . =====^==-======================— ==================== 

23 

24 not_all - 0; 7.condition if all data if imported properly 

25 

26 % [f ile path] = uiget f ile ( } . xls } ) ; 

27 

28 7. import data 

29 data = importdata( [path file]); 

30 7.Space delimeter 

31 sp_del = double (char ( } * ) ) ; 

32 /(condition varables 

33 garbage = 7; gar = 0; 

34 7#Main loop 

35 for i = 8: length (data) 

36 7.1 f data is not empty 

37 if ~isempty(data{i» 

38 %Read only data lines 

39 if char(data{i}(D) = char^CO & char (data{i> (2) ) ~= char( > *) %data line 

40 sp_ind = find(data{i> = sp_del) ; 7.Save delimeter col location 

41 sp_ind = [0 sp.ind] ; 

42 

43 for k = 1: length (sp_ind) 

44 if ~(k = 4 k i < 17) 

45 if sp_ind(k+l)-sp_ind(k) > 1 7.Test for multiple spaces 

46 start = sp_ind(k)+l; 
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47 


stop = sp_ind(k+l)-l ; 

48 


row_data{i-garbage, k-gar} = {data{i} (start : stop) >; 

49 


if (k == 3 & i < 18) | (k “ 2 & i >= 18) 

50 


tmp = double (char (row_data{i -garbage , k-gar}) ) ; 

51 


ind « find (tmp == 45); 

52 


if ~isempty(ind) 

53 


row_data{i-garbage , k-gar} = {char (tmp ( 1 : ind-1) ) } ; 

54 


end 

55 



56 


end 

57 


else 

58 


gar = gar +1; 

59 


end 

60 


else 

61 


gar = gar +1; 

62 


end 

63 



64 


%PER column 

65 


if (k+2) > length(sp_ind) 

66 


start » sp_ind(k+l)+2; 

67 


stop = length(data{i})-2; 

68 


if start > length (data{i} ( : ) ) 

69 


if "isempty (char (row_data{i-garbage , 11}) ) 

70 


per.tmp = [char(row_data{i-garbage, 10}) char(row_data{i-garbage 

71 


row_data{i-garbage , k-gar+1} = {num2str(eval(per.tmp)*100)}; 

72 


else 

73 


row_data{i-garbage} = {’’}; 

74 


end 

75 


elseif data{i} (start) == I (i-garbage) == 1 

76 


row_data{i-garbage, k-gar+1} = {’-l*}; 

77 



78 


else 

79 


row_data{i-garbage, k-gar+1} = {data{i} (start : stop)}; 

80 


end 

81 


gar = 0; 

82 


break 

83 



84 


end%End of PER column 

85 


endXEnd for loop or the row 

86 

else %If row contain labels 

87 


garbage = garbage + 1; 

88 

end 


89 

7. If 

statement for last row of data (PER column) 

90 

if 

(i == length(data)) 

91 


row_data{i-garbage , 12} = row_data{i-garbage,ll}; 
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92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 


end 

else 

not_all = 1; 

end'/, End of is data cell row empty 
end% End of row 

'/. '/.If there is a WARNING: in the last row of row_data disp and delete 
if not_all = 0 k length ( char (row_data\i-garbage, 2} )) -- length ( ' WARNING : 1 ) 

3 ■ U 
str = □ ; 

while “ isempty (row_data{i-garbage , j » 

str = [str char(row_data-Ci-garbage, j» ’ ’] ; 
row_data{i-garbage, j> = {”}; 

j = j + l 

if length (row_data) > j 
break 

end 

end 

disp(str) ; 

end 

'/.WWWWWWW Create iperf handle //////////////// 

7.T ime in seconds 

iperf .time = abs((iperf_cell_array(row_data,2))) ; 

'/.Transfer in KBytes 

iperf .transfer = iperf_cell_array(row_data,4) ; 

'/.Bandwidth in Mbit/sec 

iperf .BW = iperf_cell_array(row_data,6) ; 

'/.Packet Error Rate in percent 

iperf .PER = iperf _cell_array(row_ data, 12) ; 

'/.Make data same length 

if length (iperf .time) "= length (iperf .PER) 

len = min ([length (iperf. PER) length (iperf .time) ... 

length (iperf .BW)] ) ; 
iperf .time * iperf .timed : len) ; 
iperf . transfer = iperf . transfer (1 : len) ; 
iperf . BW = iperf . BW ( 1 : len) ; 
iperf .PER = iperf .PER(1: len) ; 

end 

if not_all = 1 

disp( ’Mat lab could not import all of the data’); 


end 
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get jperfldata.m 

1 

2 % This is called from the menu of YJ_program 

3 

4 [file path] = uigetfile; 

5 ip_fi = file;ip_pa = path; 

6 if fliplr(file(length(file)-(0:2))) == 'txt’ 

7 temp = importdata([path file]); 

8 time = temp.data(: ,1) ; 

9 BW = temp.data(: ,3) ; 

10 PER = temp.data(: ,7) ; 

11 else 

12 temp = general_iperf (file, path) ; 

13 time = temp. time; 

14 BW = temp.BW; 

15 PER = temp. PER; 

16 end 

17 

18 

19 iperf ..numeric = 2; ‘/.Conditon for iperf _gui 

20 

21 save (’ iperf _num_wksp’ , ’time' , ’BW 1 , ’PER’ , ’ iperf .numeric ’ , ... 

22 ’ ip_f i ’ , 1 ip_pa ’ ) 

23 run iperf _gui 


get jperf_data jgui . m 

1 

2 % This is called from the menu of iperf _gui under: 

3 % User Input > iperf numeric 

4 

5 [file path] = uigetfile ( , *.txt’ , ’ .xls’) ; 

6 

7 if fliplr (file (length (file) -(0:2))) == ’txt’ 

8 temp = import data ([path file]); 

9 time = temp.data( : , 1) ; 

10 BW = temp.data(: ,3) ; 

11 PER = temp.data( : ,7) ; 

12 else 

13 temp = general_iperf (file, path) 
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14 time = temp. time; 

15 BW = temp.BW; 

16 PER = temp. PER; 

17 end 

18 

19 iperf .numeric = 3; ’/.Conditon for iperf_gui menu 

20 

21 save (’ iperf _num_wksp 5 , ’time-' , ’BW’ , ’PER'* , ’ iperf .numeric ’ , 

22 ’file', ’path’) 

23 run iperf.gui 


getfolder.m 

1 function ret = getf older (path) 

2 

3 

4 ind = findCpath == ’\’); 

5 

6 ret = [path(l:3) ’...’ path(ind(length(ind)-l) : length (path)-l)] ; 


iperf_cell_array.m 

1 function ret = iperf_cell_array(data, j) 

2 

3 len = length (data) ; 

4 

5 for i = l:len 

6 if "isempty (data-fi , j}) '/A char(data-(i, j>{l» 

7 x = char(data{i, jHl» ; 

8 

9 ’/, condtion for time 

10 k = find(x == ’-’); 

11 

12 if isempty (k) 

13 

14 if length (x) == length (’sec’) ... 

15 & x == ’sec’ 

16 ret(i,l) = data{i , j+l> ; ’/.Transfer column 

17 elseif length(x) ==length( ’KBytes’) ... 

18 b x == ’KBytes’ 

19 ret(i,l) = data-Ci, j+1}; ’/.Bandwidth column 
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20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 


elseif j — 12 ft str2num(x) > 100 
ret(i,l) = data{i, j+1}; 

else 

ret{i,l> = x; ’/.If number 

end 


else 


if k = 1 

ret{i,l} = ’O’; ’/. bad data in file 

else 

ret{i,l> = x(l:k-l); ‘/.bad data import 

end 


end 

if isempty(ret{i, 1}) 
ret{i, 1} = 0; 

end 


end 

end 

ret = str2num( char (ret)) ; 
if j ==2 ft ret(len) == 0 
ret = ret(l:len-l) 

end 


'/, x= * sec ’ 

% if (length (x) == length(’sec’) I length ( 'KBytes’)) ft ... 
% (x == 'sec' I ’KBytes') 

% disp(’yes’); 

% end 


Iperf_data.m 

1 function iperf = Iperf_data(iperf ) 

2 % 

3 */. SPACE DELIMITER 

4 */. 

5 '/, USE get_iperf _data.m if possible 
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6 y================================== ================== 

7 / The program reads data files taken by iperf 

8 / and then literally copied to an Excel spreadsheat 

9 / Process : 

10 % 1) open text file of iperf data 

11 % 2) Ctl+a 

12 % 3) Ctl+c 

13 / 4) In Excel spreadsheet ctl+v 

14 7 . 5) File Save 

15 7 . Now run this program and the iperf data will be stored as followed 

16 7 . all in the iperf handle , iperf . : 

17 7 , Time in seconds (iperf. time) 

18 7 , Transfer in KBytes (iperf .transfer) 

19 7 ' Bandwidth in Mbit/sec (iperf .BW) 

20 7 . Packet Error Rate in percent (iperf .PER 

21 / 

22 / ===================== — ==================== ========= 

23 

24 

25 load iperf _loc */.pa -path and fi -file 

26 path = [pa ’iperf/’]; 

27 ind = find(fi == ’-’) 

28 file - [fi(l:ind-l) ’.xls’]; 

29 

30 

31 not_all = 0; /(condition if all data if imported properly 

32 

33 7 . [file path] = uigetfile(’ .xls’) ; 

34 

35 7 , import data 

36 data = importdata( [path file]); 

37 /.Space delimeter 

38 sp_del * double (char (’ ’)); 

39 /condition varables 

40 garbage * 7;gar = 0; 

41 /Main loop 

42 for i = 8: length (data) 

43 /If data is not empty 

44 if - isempty(data-Ci» 

45 /Read only data lines 

46 if char(data{i>(l)) = char(’[’) & char(data-[i}(2)) ~= char(’ ’) ‘/.data line 

47 sp_ind = f ind (data(i} = sp_del) ; /Save delimeter col location 

48 sp_ind = [0 sp_ind] ; 

49 

50 for k = 1 : length (sp_ ind) 
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51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 


if "(k == 4 k i < 17) 

if sp_ind (k+ 1 ) -sp_ ind (k) > 1 ‘/.Test for multiple spaces 
start = sp_ind(k)+l; 
stop = sp_ind(k+l)-l; 

row_data-Ci-garbage, k-gar} = {data{i} (start : stop) > ; 
if (k = 3 k i < 18) | (k = 2 ft i >= 18) 

tmp = double (char (row_data{i-garbage , k-gar}) ) ; 
ind = find (tmp == 45); 
if 'isempty(ind) 

row_data{i-garbage, k-gar} = {char(tmp(l:ind-l))}; 

end 


end 

else 

gar = gar +1; 

end 

else 

gar = gar +1; 

end 

'/(PER column 

if (k+2) > length (sp_ind) 
start = sp_ind(k+l)+2; 
stop = length (data{i}) -2; 
if start > length (data{i}(:)) 

if "isempty(char(row_data{i-garbage, 11}) ) 

per_tmp = [char(rov_data{i-garbage, 10}) char(row_data{i-garbage 
row_data{i-garbage,k-gar+l} = {num2str(eval(per_tmp)*100)}; 

else 

row_data{i-garbage} = { ’ ’ } ; 

end 

elseif data{i} (start) == | (i-garbage) == 1 

row_data{i-garbage,k-gar+l} = {*-1*}; 

else 

row_data{i-garbage , k-gar+1} = {data{i} (start: stop)}; 

end 

gar = 0; 
break 

end‘/.End of PER column 
end%End for loop or the row 
else %If row contain labels 
garbage = garbage + 1; 

end 
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96 %If statement for last row of data (PER column) 

97 if (i == length (data)) 

98 row_data{i-garbage,12> = row_data{i-garbage,ll>; 

99 end 

100 else 

101 not_all = 1; 

102 end% End of is data cell row empty 

103 end/* End of row 

104 

105 7. y.If there is a WARNING: in the last row of row.data disp and delete 

106 if not_all 0 k length ( char (row_data{i -garbage, 2})) == length ( * WARNING : ’ ) 

107 j = 1; 

108 str = □ ; 

109 while ~isempty (row_data{i-garbage , j}) 

110 str = [str char (row_data{i-garbage , j } ) 3 *] ; 

111 row_data{i-garbage, j} = { ,J >; 

112 j = j+1; 

113 end 

114 disp(str); 

115 end 

116 

117 %\\\\\\\\\\\\\\ Create iperf handle //////////////// 

118 %Time in seconds 

119 iperf .time = abs(round(iperf_cell_array(row_data,2))) ; 

120 ^Transfer in KBytes 

121 iperf .transfer = iperf_cell_array(rov_data,4) ; 

122 ^Bandwidth in Mbit /sec 

123 iperf. BW = iperf_cell_array(row_data,6) ; 

124 %Packet Error Rate in percent 

125 iperf .PER = iperf_cell_array(row_data, 12) ; 

126 TMake data same length 

127 if length (iperf .time) < length (iperf .PER) 

128 len « length (iperf .time) ; 

129 iperf .time * iperf .timed : len) ; 

130 iperf .transfer = iperf .transferd: len) ; 

131 iperf. BW « iperf .BW(1: len); 

132 elseif length (iperf .PER) < length (iperf .time) 

133 len = length(iperf .PER) ; 

134 iperf .time = iperf .timed : len ) ; 

135 iperf .transfer - iperf .transferd : len) ; 

136 iperf. BW = iperf .BW(1: len) ; 

137 end 

138 

139 if not_all — 1 

140 dispdMatlab could not import all of the data’); 
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141 

142 end 

143 

144 '/, handles . iperf .figure = figure 

145 V, set (gcf , ’ Name ’ , handles . filename) ; 

146 % Chares, h_per,h_bw] = plotyy(iperf .time, iperf .PER, iperf. time, iperf .BW) 

147 '/. 

148 '/, axes Chares ( 1 ) ) JCPER parameters 

149 '/. ylabelC ’Packet Error Rate’); 

150 '/, ylimCCO 100]); 

151 % set (h_per , ’ Color ’ , ’ r ’ ) ; 

152 axes (haxes (2) ) %BW parameters 

153 '/, ylabelC ’Bandwidth [Mbits/sec]’) 

154 ylimCCO max C iperf .BW) +20]) 

155 set Ch_bv , ’ LineStyle ’ , ’ — ’ ) ; 


iperf_file.m 

1 

2 y. This is called from the menu of YJ_program under: 

3 '/, File > Iperf 

4 this command finds the iperf data for the current YJ data file 

5 load handles ; 

6 ip_fi = handles. file; 

7 ip_pa = handles. path; 

8 

9 iperf _numeric =1; '/.Conditon for iperf _gui menu 

10 

11 save C ’ iperf _num_wksp ’ , ’ iperf .numeric ’ , ... 

12 ’ip_fi’, ’ip_pa’) 

13 run iperf _gui 


iperf-gui.m 

1 function varargout = iperf _gui Cvarargin) 

2 '/, IPERF.GUI M-file for iperf _gui. fig 

3 '/, IPERF.GUI, by itself, creates a new IPERF.GOI or raises the existing 

4 ’/. singleton* . 

5 y. 

6 A H = IPERF.GUI returns the handle to a new IPERF.GUI or the handle to 

7 '/, the existing singleton*. 
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8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 


7. 

% IPERF_GUI (’ CALLBACK’ .hObject, eventDat a, iperf, .. .) calls the local 

/• function named CALLBACK in IPERF_GUI.M with the given input arguments 

% 

'/, IPERF_GUI ( * Property ’ , ’Value ’,...) creates a new IPERF.GUI or raises the 

% existing singleton*. Starting from the left, property value pairs are 

% applied to the GUI before iperf _gui_0peningFunction gets called. An 

% unrecognized property name or invalid value makes property application 

V. stop. All inputs are passed to iperf _gui_0peningFcn via varargin. 

% 

/• *See GUI Options on GUIDE’S Tools menu. Choose "GUI allows only one 

V. instance to run (singleton)". 

7 . 

7. See also: GUIDE, GUIDATA, GUIHANDLES 


% Edit the above text to modify the response to help iperf_gui 
'/. Last Modified by GUIDE v2.5 01-Nov-2004 16:56:26 


V. Begin initialization code - DO NOT EDIT 
gui_Singleton = 1; 

gui_State = struct (’ gui_Name ’ , mfilename, ... 

’gui_Singleton’ , gui_Singleton, ... 

’ gui_0peningFcn ’ , ®iperf_gui_0peningFcn, 
’gui_0utputFcn’ , ®iperf_gui_0utputFcn, . 
’gui_LayoutFcn’ , □ , . . . 

’ gui_Callback ’ , □ ) ; 

if nargin & isstr(varargin{l}) 

gui_State . gui_Callback = str2func(varargin-(l}-) ; 


if nargout 

[varargout{l:nargout>] = gui_mainf cn(gui_State, varargin{:»; 

else 

gui_mainf cn (gui_State , varargin{ : » ; 

end 

'/, End initialization code - DO NOT EDIT 


% Executes just before iperf _gui is made visible. 

function iperf _gui_0peningFcn(h0bject, eventdata, iperf, varargin) 

'/, This function has no output args, see OutputFcn. 

7, hObject hand! e to figure 

V, eventdata reserved - to be defined in a future version of MATLAB 
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53 

54 

55 

56 

57 

58 

59 

60 
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62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 
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81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 


7. iperf structure with iperf and user data (see GUIDATA) 

'/, varargin command line arguments to iperf_gui (see VARARGIN) 

V.Creat Menu OptionsO 

load iperf _num_wksp 

if iperf _numeric ”= 3 %If called from YJ_workl 
'I, iperf_gui menu 

file = uimenu( ’Label * , ’File ’ ) ; 

uimenu(f ile, ’Label’ , ’Open’ , ’Callback’ , ’get_iperf_data_igui’) ; 
uimenu(f ile, ’Label’ , ’Save Plot’ , ’Callback’ , ’save’) ; 
uimenu(file, ’Label’ , ’Quit ’ , ’Callback’ , ’exit ’ , . . . 

’Separator’ , ’on’ , ’Accelerator’ , ’Q’) ; 

PP = uimenu( ’Label ’, ’Plot Properties’); 
uimenu(PP, ’Label’, ’Grid’, ’Callback’, ’Grid_stat’) ; 
uimenu(PP, ’Label’, ’Clear Graph’, ’Callback’, ’cla’); 

help = uimenu ( ’ Label ’ , ’ Help ’ ) 

uimenu(help, ’Label ’ , ’Help’ , ’Callback’ , ’help_Doc’ ) ; 

end 

if iperf_numeric == 1 
load handles; 

set(gcf , ’Name’ .handles. filename) ; 
iperf .rssi = handles. rssi; 
iperf = Iperf_data(iperf ) ; 
set (iperf . rssi_but , ’Value ’ , 1) ; 
set (iperf .rssi_but, ’Visible’ , ’on’) 

else 

set (gcf , ’ Name ’ , ip_f i) 

iperf .time = time 
iperf .BW = BW 
iperf .PER * PER 

iperf .rssi_val = 0; 

set (iperf .rssi_but, ’Value’ ,0) ; 

set (iperf . rssi_but , ’Visible ’ , ’ off ’ ) ^ 

end I 

delete_file = 0; 

save ( * iperf _num_wksp 7 , 7 delete_f ile 7 ) ; 
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98 

99 

100 iperf.BW = BW_filter (iperf .BW) ; 

101 

102 save iperf 

103 re_plot (iperf ) 

104 

105 '/, Choose default command line output for iperf_gui 

106 ’/.iperf . output = hObject; 

107 

108 V, Update iperf structure 

109 guidata(h0bject, iperf); 

110 

111 7, UIWAIT makes iperf_gui wait for user response (see UIRESUME) 

112 '/, uiwait (iperf .figurel); 

113 

114 

115 Outputs from this function are returned to the command line. 

116 */. function varargout = iperf _gui_0utputFcn(h0bject, eventdata, iperf) 

117 % */, varargout cell array for returning output args (see VARARGOUT); 

118 '/, '/, hObject handle to figure 

119 ’/, 1 eventdata reserved - to be defined in a future version of MATLAB 

120 ‘/. y, iperf structure with iperf and user data (see GUIDATA) 

121 •/. 

122 '/, y. Get default command line output from iperf structure 

123 '/, varargout{l)- = iperf . output ; 

124 

125 

126 Executes on button press in txt_per_but. 

127 function per_but_Callback (hObject, eventdata, iperf) 

128 % L hObject handle to txt_per_but (see GCB0) 

129 eventdata reserved - to be defined in a future version of MATLAB 

130 'I, iperf structure with iperf and user data (see GUIDATA) 

131 

132 re_plot (iperf) 

133 y. Hint: get (hObject , 'Value ’ ) returns toggle state of txt_per_but 

134 

135 

136 y, Executes on button press in txt_rssi_but . 

137 function rssi_but_Callback(hObject, eventdata, iperf) 

138 '/, hObject handle to txt_rssi_but (see GCB0) 

139 ’/, eventdata reserved - to be defined in a future version of MATLAB 

140 '/. iperf structure with iperf and user data (see GUIDATA) 

141 re_plot ( iperf ) 

142 ’/, Hint: get (hObject, ’Value’) returns toggle state of txt_rssi_but 
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143 

144 

145 

146 

147 

148 

149 
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168 

169 

170 

171 
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180 
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187 


7. Executes on button press in bv_but. 

function bw_but ..Callback (hObject, eventdata, iperf) 

7, hObject handle to bv_but (see GCBO) 

V. eventdata reserved - to be defined in a future version of MATLAB 
7, iperf structure with iperf and user data (see GUIDATA) 
re.plot (iperf ) 


7, Hint: get (hObject, ’Value’) returns toggle state of bw_but 



function re_plot (iperf) 
load iperf 
load handles 


Xclear statistics 
set(iperf -mean_rssi, ’String’ , ”) ; 
set (iperf .mean_bw, ’String’ , ; 

set (iperf .mean_per, ’String’ , ” ) ; 
set (iperf .var_rssi, ’String’ , ”) ; 
set (iperf . var_bw, ’String’ , ” ) ; 
set ( iperf . var _per , ’ String ’ , ” ) ; 
opt = 0; 

bw_val = get (iperf . bw_but , ’ Value ’ ) ; 
per_val = get (iperf . per_but , ’ Value ’ ) ; 
rssi_val = get (iperf ,rssi_but, ’Value’); 
subplot (1,1,1) 
if rssi_val — 0 

Xclear txt_rssi statistics 

set (iperf .mean_rssi, ’Visible’ , ’off ’) ; 

set (iperf . var.rssi , ’Visible ’ , ’ off ’ ) ; 

subplot (1,1,1) 

if (get (iperf .bw_but, ’Value’) & get (iperf .per.but, ’Value’)) == 1 
‘/.Plot TXT_PER and BW information 

[haxes , h_per , h_bw] = plotyy ( iperf . time , 100-iperf . PER , ... 

iperf . time ( 1 : length (iperf . BW) ) , iperf . BW) ; 
title (’Datagrams Received [*/.] and BW [Mbits/sec]’); 

'/.Axis Label 

set (get (haxes (1) , ’Ylabel’) , ’String’ , ’Percent’) 
set (get (haxes (2) , ’Ylabel’) , ’String’ , ’Mbit/sec’) 

%Axis limits 
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209 

210 
211 
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214 
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ylim(haxes(l) , [0 105]); 
ylim(liaxes(2) , [0 10] ) ; 

'/.Line Properties 

set (h_per , ’ Color ’ , ’ b ’ ) ; 
set (h_bw , ’ LineStyle — ’ ) ; 

%Set GUI Visulation Parameters 

set (iperf .mean_bv, ’Visible ’ , ’off ’ ) ; 
set (iperf .var_bv, ’Visible’ , ’off’); 

set (iperf .mean_bw, ’Visible’ , ’on’) ; 
set (iperf .var_bv, ’Visible’ , ’on’) ; 

set (iperf .mean.bw, ’String’ ,num2str(mean(iperf .BW))) ; 
set (iperf .var_bv, ’String’ . num2str (var ( iperf . BW) ) ) ; 

set (iperf .mean_per, ’Visible’ , ’on’) ; 
set (iperf . var.per , ’Visible ’ , ’ on ’ ) ; 

set (iperf .mean_per, ’String’ , num2str (100-mean (iperf . PER) ) ) ; 
set (iperf .var.per, ’String’ , num2str(var( iperf .PER))) ; 

time = iperf .time; 

PER - iperf .PER; 

BW = iperf .BW; 
file = ip_fi; 

save ( ’ current_data’ , ’time’ . ’PER’ , ’BW’ , ’file’) ; 

elseif get (iperf .per_but, ’Value’) == 1 k get ( iperf . bw_but , ’ Value ’ ) ==0 
‘/.Plot TXT_PER 

plot(iperf .time, 100-iperf. PER) ; title (’Datagrams Received [*/,]’); 
ylabel ( ’Percent ’ ) ; 
ylim( [0 105]); 

'/.Set GUI Visulation Parameters 
set (iperf .mean.bir, ’Visible’ , ’off ’) ; 
set (iperf .var _bw,’ Visible’ , ’off ’); 

set (iperf .mean_per, ’Visible ’ , ’on’ ) ; 
set (iperf .var.per, ’Visible’ , ’on’) ; 

set (iperf .mean_per, ’String’ ,num2str( 100-mean (iperf .PER))) ; 
set (iperf .var _per, ’String’ ,num2str(var(iperf .PER))) ; 

time = iperf .time; 

PER = iperf .PER; 
file = ip_fi; 

save (’ current _data’ , ’time’ , ’PER’ , ’file’); 
elseif get ( iperf . per_but , ’ Value ’ ) == 0 & get ( iperf . bw_but ,’ Value ’ ) ==1 
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'/.Plot BW 

plot(iperf .time, iperf .BW) ; title (’BW [Mbits/sec] ’); 
axis ([0 lengthCiperf .time) 0 10]); 

‘/.Set GUI Visulation Parameters 
set (iperf .mean_bv, ’Visible ’ , ’on* ) ; 
set (iperf . var_bw, ’Visible ’ , ’on’ ) ; 

set (iperf .mean_bw, ’String’ , num2str (mean ( iperf . BW) ) ) ; 
set (iperf .var_bw, ’String’ , num2str(var( iperf .BW))) ; 

set (iperf .mean_per, ’Visible’ , ’off ’) ; 
set (iperf .var_per, ’Visible’ , ’off ’) ; 

time = iperf .time; 

BW = iperf. BW; 
file = ip_fi; 

save (’ current _data’ , ’time’ , *BW’ , ’file’) ; 

elseif (get (iperf .bw.but, ’Value’) k get (iperf .per _but, ’Value’ )) == o 
subplot (1,1,1) 
plot (0,0) 

t it le(’ Nothing Selected’); 

%Set GUI Visulation Parameters 
set (iperf .mean_per , ’Visible ’ , ’ off ’ ) ; 
set (iperf .var_rssi, ’Visible’ , ’off ’) ; 
set (iperf . var_bv, ’Visible’ , ’off’) ; 
set (iperf .var^per, ’Visible’ , ’off’) ; 

end 

elseif rssi_val == 1 

XSet txt_rssi statistics 

set (iperf .mean_rssi, ’Visible’ , ’on’) ; 

set (iperf .var_rssi, ’Visible’ , ’on’); 

set (iperf .mean_rssi, ’String’ ,num2str (mean (iperf .rssi))) ; 
set ( iperf .var_rssi, ’String’ , num2str(var (iperf .rssi))) ; 

subplot (211) 

if (get (iperf .bw_but, ’Value’) & get (iperf .per _but, ’Value’)) == 1 

Diaxes, h_per,h_bw] = plot yy( iperf .time, 100-iperf. PER, ... 

iperf .time , iperf . BW(1 : length(iperf . time) ) ) ; 
title (’Datagrams Received [7.] and BW [Mbits/sec]’); 

XAxis Label 

set (get (haxes(l) , ’Ylabel’) , ’String’ , ’Percent’) 
set (get (haxes(2) , ’Ylabel’) , ’String’ , ’Mbit/sec’) 

Vjbcis Limits 

ylim(haxes(l) , [0 105]); 
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ylim(haxes(2) , [0 10] ) ; 

'/.Line Properties 

set (h_per , ’ Color ’ , ’ b * ) ; 
set (h_bv, ’LineStyle — ’ ) ; 

‘/.Set GUI Visulation Parameters 
set(iperf .mean_bv, ’Visible* , ’on*) ; 
set(iperf .var_bw, ’Visible’ , ’on’) ; 

set (iperf . mean_bw , ’ String ’ , num2str (mean (iperf . BW) ) ) ; 
set (iperf . var_bw , ’ String ’ , num2str (var ( iperf . BW) ) ) ; 

set (iperf .mean_per, ’Visible’, ’on’) ; 
set (iperf .var _per, ’Visible’ ,’on’); 

set (iperf .mean_per, ’String’ ,num2str(100-mean(iperf .PER))) ; 
set (iperf . var.per , ’String ’ ,num2str (var (iperf . PER) ) ) ; 

time = iperf .time; 

PER = iperf .PER; 

BW = iperf. BW; 
rssi = iperf .rssi; 
file = ip_fi; 

save (’ current _data’ , ’time* , ’PER’ , ’BW* , ’rssi’ , ’file’) ; 

elseif get (iperf .per_but, ’Value’) == 1 & get (iperf .bw_but, ’Value’) ==0 

plot(iperf .time, 100-iperf. PER) ;title( ’Datagrams Received ['/,]’); 
ylabel ( ’Percent ’ ) ; 
ylim([0 105]); 

'/(Set GUI Visulation Parameters 
set (iperf .mean_bw, ’Visible’ , ’off ’) ; 
set (iperf .var _bw, ’Visible’ , ’off’); 

set (iperf .mean_per , ’Visible ’,’on’); 
set (iperf . var_per , ’Visible ’,’on’); 

set (iperf .mean_per, ’String’ , num2str( 100-mean (iperf .PER)) ) ; 
set (iperf . var_per , ’ String ’ , num2str (var ( iperf . PER) ) ) ; 

time = iperf .time; 

PER = iperf .PER; 
rssi = iperf. rssi; 
file = ip_fi; 

save(’current_data’ , ’time’ , ’PER’ , ’rssi’ , ’file’) ; 

elseif get (iperf .per _but, ’Value’) == 0 & get (iperf .bw_but, ’Value’) ==1 
plot (iperf .time, iperf .BW) ; title (’BW [Mbits/sec] ’); 

'/.Set GUI Visulation Parameters 
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set(iperf .mean_bw, ’Visible’ , ’on’) ; 
set(iperf .var_bw, ’Visible’ , ’on’) ; 

set(iperf .mean_bw, ’String’ , num2str (mean (iperf .BW))) ; 
set (iperf . var_bw , ’ String ’ , num2str (var (iperf . BW) ) ) ; 

set (iperf .mean_per, ’Visible’ , ’off ’) ; 
set (iperf . var_per , ’Visible ’ , ’off ’ ) ; 

time = iperf. time; 

BW = iperf .BW; 
rssi = iperf. rssi; 
file * ip_fi; 

save (’ current _data’ , ’time’ , ’rssi’ , ’BW’ , ’file’) ; 

elseif (get ( iperf .bv_but, ’Value’) & get(iperf .per _but, ’Value’ )) == 0 
subplot (1,1,1) 

plot (iperf .rssi) ;title(’RSSI’) ; 

ylabel(’dBm’) 

opt = 1; 

%Set GUI Visulation Parameters 
set (iperf .mean_bw, ’String ’ , ” ) ; 
set (iperf . var_bw , ’ String ’,’’); 

set (iperf .mean_per, ’String’ , ”) ; 
set (iperf .var _per, ’String’ , ”) ; 

end 

if opt == 0 

subplot (212) 

plot (iperf .rssi) ; title ( ’RSSI’) ; 
ylabel(’dBm’) 


end 

end 


y. Executes during object creation, after setting all properties. 

function mean_txt_rssi_CreateFcn(hObject, eventdata, handles) 
hObject handle to mean_txt_rssi (see GCBO) 

*/, eventdata reserved - to be defined in a future version of MATLAB 

'/. handles empty - handles not created until after all CreateFcns called 

'/, Hint: edit controls usually have a white background on Windows. 

See ISPC and COMPUTER. 

if ispc 
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368 set (hObject, 'BackgroundColor' , 'white') ; 

369 else 

370 set (hObject , 'BackgroundColor' ,get(0, 'def aultUicontrolBackgroundColor ' )) ; 

371 end 

372 

373 

374 

375 function mean_txt_rssi_Callback(hGbject , eventdata, handles) 

376 7 hObject handle to mean_txt_rssi (see GCBO) 

377 7 eventdata reserved - to be defined in a future version of MATLAB 

378 7 handles structure with handles and user data (see GUIDATA) 

379 

380 7 Hints: get (hObject, 'String') returns contents of mean_txt_rssi as text 

381 7 str2double (get (hObject, 'String')) returns contents of mean_txt_rssi as a double 

382 

383 

384 7 Executes during object creation, after setting all properties. 

385 function var_txt_rssi_CreateFcn (hObject , eventdata, handles) 

386 7 hObject handle to var_txt_rssi (see GCBO) 

387 7 eventdata reserved - to be defined in a future version of MATLAB 

388 7. handles empty - handles not created until after all CreateFcns called 

389 

390 7, Hint: edit controls usually have a white background on Windows. 

391 7. See ISPC and COMPUTER. 

392 if ispc 

393 set (hOb j ect , ' BackgroundColor ' , ' white ' ) ; 

394 else 

395 set (hObject, 'BackgroundColor' ,get(0, 'def aultUicontrolBackgroundColor ' )) ; 

396 end 

397 

398 

399 

400 function var__txt_rssi_Callback (hObject, eventdata, handles) 

401 7 hObject handle to var_txt_rssi (see GCBO) 

402 7. eventdata reserved - to be defined in a future version of MATLAB 

403 7. handles structure with handles and user data (see GUIDATA) 

404 

405 7 Hints: get (hObject, 'String') returns contents of var_txt_rssi as text 

406 7 str2double (get (hObject, 'String')) returns contents of var_txt_rssi as a double 

407 

408 

409 7 Executes during object creation, after setting all properties. 

410 function var_bw_CreateFcn (hObject, eventdata, handles) 

411 7 hObject handle to var_bw (see GCBO) 

412 7 eventdata reserved - to be defined in a future version of MATLAB 
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413 7 . handles empty - handles not created until after all CreateFcns called 

414 

415 7 . Hint: edit controls usually have a white background on Windows. 

416 7 . See ISPC and COMPUTER. 

417 if ispc 

418 set (hOb j ect , * Backgr oundColor ’ , ’ white ’ ) ; 

419 else 

420 set (hObject , ’Backgr oundColor ’ ,get(0, : def aultUiconrroiBackgroundCoior ’ )) ; 

421 end 

422 

423 

424 

425 function var_bw_Callback (hObject, eventdata, handles) 

426 7 , hObject handle to var_bw (see GCBO) 

427 7 . eventdata reserved - to be defined in a future version of MATLAB 

428 7 . handles structure with handles and user data (see GUIDATA) 

429 

430 7 , Hints: get (hObject, ’String’) returns contents of var_bw as text 

431 str2double (get (hObject, ’String’)) returns contents of var_bw as a double 

432 

433 

434 7 , Executes during object creation, after setting all properties. 

435 function mean_bw_CreateFcn (hObject, eventdata, handles) 

436 7 % hObject handle to mean_bw (see GCBO) 

437 7 % eventdata reserved - to be defined in a future version of MATLAB 

438 7 , handles empty - handles not created until after all CreateFcns called 

439 

440 7 » Hint: edit controls usually have a white background on Windows. 

441 7 . See ISPC and COMPUTER. 

442 if ispc 

443 set (hOb j ect , ’ BackgroundColor ’ , ’ white ’ ) ; 

444 else 

445 set (hObject, ’BackgroundColor’ ,get(0, ’defaultUicontrolBackgroundColor’)) ; 

446 end 

447 

448 

449 

450 function mean_bw_Callback (hObject, eventdata, handles) 

451 7 % hObject handle to mean_bw (see GCBO) 

452 7 . eventdata reserved - to be defined in a future version of MATLAB 

453 7 . handles structure with handles and user data (see GUIDATA) 

454 

455 7 % Hints: get (hObject, ’String’) returns contents of mean_bw as text 

456 7 . str2double (get (hObject, ’String’)) returns contents of mean_bw as a double 

457 
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458 

459 7 * Executes during object creation, after setting all properties. 

460 function var_txt_per_CreateFcn (hObject , eventdata, handles) 

461 7. hObject handle to var_txt_per (see GCBQ) 

462 % eventdata reserved - to be defined in a future version of MATLAB 

463 7 handles empty - handles not created until after all CreateFcns called 

464 

465 7. Hint: edit controls usually have a white background on Windows. 

466 7. See ISPC and COMPUTER. 

467 if ispc 

468 set (hObject, ’ BackgroundColor ’ , ’white’) ; 

469 else 

470 set(hObject, ’BackgroundColor’ ,get(0, ’def aultUicontrolBackgroundColor ’ )) ; 

471 end 

472 

473 

474 

475 function var_txt_per_Callback (hObject, eventdata, handles) 

476 % hObject handle to var_txt_per (see GCBO) 

477 % eventdata reserved - to be defined in a future version of MATLAB 

478 7 * handles structure with handles and user data (see GUIDATA) 

479 

480 7. Hints: get (hObject, ’String’) returns contents of var_txt_per a s text 

481 % str2double (get (hObject, ’String’)) returns contents of var_txt_per as a double 

482 

483 

484 7 % Executes during object creation, after setting all properties. 

485 function mean_txt _per_CreateFcn (hObject , eventdata, handles) 

486 7 % hObject handle to mean_txt_per (see GCBO) 

487 % eventdata reserved - to be defined in a future version of MATLAB 

488 7 * handles empty - handles not created until after all CreateFcns called 

489 

490 7 . Hint: edit controls usually have a white background on Windows. 

491 % See ISPC and COMPUTER. 

492 if ispc 

493 set (hOb j ect , ’ BackgroundColor ’ , ’ white ’ ) ; 

494 else 

495 set (hOb j ect , ’ BackgroundColor ’ , get (0 , * def aultUicontrolBackgroundColor ’ ) ) ; 

496 end 

497 

498 

499 

500 function mean_txt_per_Callback (hObject, eventdata, handles) 

501 % hObject handle to mean_txt_per (see GCBO) 

502 7. eventdata reserved - to be defined in a future version of MATLAB 
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7 . handles structure with handles and user data (see GUIDATA) 

% Hints: get (hObject, ’String’) returns contents of mean_txt_per as text 
7 , str2double (get (hObject, ’String’)) returns contents of mean_txt_per as a double 

% Executes on button press in Ext_fig. 

function Ext_fig_Callback (hObject, eventdata, handles) 

7 % hObject handle to Ext_fig (see GCB0) 

7 . eventdata reserved - to be defined in a future version of MATLAB 
7 . handles structure with handles and user data (see GUIDATA) 
load current_data 

name = figure 

set (name , ’ Name ’ , f ile ) 

set (name , ’ NumberTitle ’ , ’ of f ’ ) 


if exist(’PER’) k exist(’BW’) 

•/.Plot TXT_PER and BW information 

[haxes,h_per,h_bw] = plotyy(time,100-PER, ... 
time (1 : length (BW) ) ,BW) ; 

title ( ’Datagrams Received [/(] and BW [Mbits/sec]’); 
/(Axis Label 

set(get(haxes(l) , ’Ylabel’) , ’String’ , ’Percent’) 
set (get (haxes (2) , ’Ylabel’) , ’String’ , ’Mbit/sec’) 
/(Axis limits 

ylim(haxes(l) , [0 105]); 
ylim(haxes(2) , [0 10] ) ; 

/(Line Properties 

set (h_per , ’ Color ’ , ’ b ’ ) ; 
set (h_bw , ’LineStyle 

legend ( [h_per h_bw] , ’Datagram’ , ’Bandwidth’ ) 
elseif exist (’PER’) 

'/(Plot TXT.PER 

plot (time, 100-PER) ; title (’Datagrams Received [*/,]’); 
title ( ’Datagrams Received [%] ’ ) ; 
ylabel ( ’Percent ’ ) ; 
ylim( [0 105]); 
legend ( ’ Percent ’ ) 
elseif exist (’BW’) 

‘/Plot BW 

plot (time, BW) ; title (’BW [Mbits/sec] ’) ; 
title(’BW [Mbits/sec] ’ ) ; 
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548 legend ( ’ Bandwidth ’ ) ; 

549 

550 end 


iperfLpre.m 

1 load iperf 

2 if ishandle (iperf .figure 1) 

3 closeCiperf .f igurel) 

4 end 

5 

6 run iperf _gui ; 


mouse_data.m 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


/Sebastian Stewart 

— - 

‘/.This program is what allows the user to 
bright click on the delay spread plot and 
/select function to execute 
/=======——=—=================== 

cmenu = uicontextmenu; 

mu = handles . DS_mean* . . . 

ones ( 1 , length (handles . current _data) ) ; 

handles . ds_plot = plot (handles . current _data , ’ UIContextMenu ’ , cmenu) ; 
hold on; 

handles . sds_plot = plot(mv_data, ’r’ , ’UIContextMenu’ .cmenu) ; 
handles. mu.plot = plot (mu, ’g’ ... 

, ’ UIContextMenu ’ , cmenu) ; 
hold off; 

set(cmenu, ’Tag’ , char (handles. ssid_name) ) ; 

'/these are the funciton that execute when item are selected 
cbl = [’this’]; 
cb2 = [’sds_menu (handles) ’] ; 
cb3 = [’uiresume’] ; 

/These cure what get displayed when right click of mouse on object 

iteml = uimenu(cmenu, ’Label’, ’This’, 

’Callback’, cbl); 

item2 = uimenu(cmenu, ’Label’, ’MA Delay Spread’, ... 
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27 ’Callback’, cb2) ; 

28 item3 = uimenuCcmenu, ’Label’, ’resume’, ... 

29 ’Callback’, cb3) ; 


new_file.m 


1 

2 [iperf.file iperf.path] = uigetfileC’*.txt’) ; 

3 setCgcf, ’Name’ .iperf .file) ; 

4 temp = importdataC [iperf .path iperf.file]); 

5 

6 iperf.time = temp.dataC: ,1) ; 

7 iperf.BW = temp.dataC: ,3) ; 

8 iperf.PER = temp.dataC: ,7) ; 

9 

10 load handles ; 

11 iperf .rssi = handles . rssi ; 

12 

13 iperf = Iperf _dataCiperf) 

14 save iperf 

15 re_plot C iperf ) 


re_plot.m 


ssidJindnew.m 

1 function [ssid_names, ssid_ index] - ssid_findnewCdata,col) 

2 

3 %===—==== 

4 ‘/.Assumptions The data is sorted by channel and SSID 

5 ‘/.this program:] 

6 % 1) outputs the index of each ssid 

7 7 , 2) output cell array of the ssid in the data 

8 %=================================================== 

9 

10 ’/, initial = 0 ’/.searching for new ssid name 

11 '/, start = [] ; 

12 
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13 % '/.set initial ssid name 

14 ssid_names = O; 

15 % %set initial ssid index 

16 ssid_ index = [] ; 

17 

18 %Search all rows 

19 for i = 2: length (data) 

20 

21 if char (data{i, col}) 

22 

23 %if strings are same 

24 1 ength( char (data{i, col})) ; 

25 length (char (data{i-l, col}) ) ; 

26 char (data{i , col}) ; 

27 char (data{i-l, col}) ; 

28 

29 if length (char (data{ i , col}) ) length (char (data{i-l, col})) 

30 

31 s s id_index (length ( ss id_index) + 1 ) = i-1; 

32 ssid_index(length(ssid_index)+l) * i; 

33 ssid_names{length(ssid_names)+l, 1} = data(i,col); 

34 

35 elseif length (char (datali, col}) )== length ( char (datali-1, col})) ... 

36 t char(data{i, col}) ~= char (dat a{ i- 1 , col}) 

37 ssid_index(length(ssid_index)+l) * i-1; 

38 ssid_index(length(ssid_index)+l) = i; 

39 ssid_names{length(ssid_names)+l, 1} = data(i,col); 

40 end 

41 end 

42 end 

43 ssid_index(length(ssid_index)+l) = i; 

44 

45 


time_sample_data.m 

1 function handles = time_sample_data(handles) 

2 ‘/.handles used: row handles created: tsmpl 

3 user_str = {'user*}; 

4 

5 tm_per = . 1 : . 1 : 1 

6 ssid_len = handles. row (2) -handles .row (1) 

7 time_pts * unique (ceil (ssid_len.*tm_per)) 
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8 handles. tsmpl = time_pts(l) 

9 

10 7. time_pts = cellstr(time_pts) 

11 y = cellstr(time_pts) ; 

12 time_pts = cat(2,y,user_str) ; 

13 

14 set (handles. dwn_samp, ’String’ ,time_pts) 

15 

16 


YJ_workl.m 

1 function varargout = YJ_workl (varargin) 

2 

3 7.for_data (90,184) 

4 %data_fetch (445) 

5 7.time_sample_data (356, 421) 

6 7delay_spread (371,421,439) 

7 ‘/.col.find (99,190) 

8 */.cell_find (188) 

9 

10 

11 

12 

13 V, YJ_W0RK1 M-file for YJ_workl.fig 

14 7 % YJ_W0RK1, by itself, creates a new YJ_W0RK1 or raises the existing 

15 7 , singleton*. 

16 */. 

17 % H = YJ_W0RK1 returns the handle to a new YJ_W0RK1 or the handle to 

18 7 , the existing singleton*. 

19 */. 

20 % YJ_W0RK1(’ CALLBACK’ .hObject, event Data, handles, .. .) calls the local 

21 7 % function named CALLBACK in YJ_W0RK1.M with the given input arguments . 

22 7 . 

23 7t YJ_W0RK1 ( ’Property’ , ’Value ’,...) creates a new YJ_W0RK1 or raises the 

24 7 * existing singleton*. Starting from the left, property value pairs are 

25 7 * applied to the GUI before YJ_workl_0peningFunction gets called. An 

26 7 , unrecognized property name or invalid value makes property application 

27 7 % stop. All inputs are passed to YJ_workl_0peningFcn via varargin. 

28 7. 

29 7. *See GUI Options on GUIDE’S Tools menu. Choose "GUI allows only one 

30 7 , instance to run (singleton)". 

31 7. 
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32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 


*/. See also: GUIDE, GUIDATA , GUIHANDLES 

*/, Edit the above text to modify the response to help YJ_workl 

’/. Last Modified by GUIDE v2.5 24-Jul-2004 14:29:00 

% Begin initialization code - DO NOT EDIT 
gui_Singleton = 1; 

gui_State = struct (’ gui_Name ’ , mfilename, ... 

’gui_Singleton’ , gui_Singleton, ... 

’ gui_OpeningFcn * , ®YJ_workl_0peningFcn, . . . 
’gui_0utputFcn’ , ®YJ_workl_0utputFcn, ... 

’ gui_LayoutFcn ’ , [] , ... 

’ gui_Callback ’ , [] ) ; 

if nargin & isstr(varargin{l» 

gui_State . gui.Callback = str2f unc (varargin{l» ; 

end 

if nargout 

[varargout-Cl : nargout}] - gui_mainf cn(gui_State, varargin-C : }) ; 

else 

gui_mainf cn(gui_State , varargin!:}) ; 

end 

’/, End initialization code - DO NOT EDIT 

V, Executes just before YJ_workl is made visible. 

function YJ_workl_0peningFcn(h0bject, eventdata, handles, varargin) 
V, This function has no output args, see OutputFcn. 

% hObject handle to figure 

% eventdata reserved - to be defined in a future version of MATLAB 
% handles structure with handles and user data (see GUIDATA) 

*/, varargin command line arguments to YJ_workl (see VARARGIN) 

/(This is not yet working completly 

V.Creat Menu OptionsO 

file = uimenu (’ Label’ , ’File’) ; 

uimenu(f ile, ’Label’ , ’Iperf ’ , ’Callback’ , ’iperf_f ile’ ) ; 
uimenu(file, ’Label’ , ’Save Plot’ , ’Callback’ , ’save’) ; 
uimenu(file, ’Label’ , ’Quit’ , ’Callback’ , ’exit’ , . . . 

’Separator’ , ’on’ , ’Accelerator’ , ’Q’) ; 

PP = uimenu ( ’ Label ’ , ’ Plot Properties ’ ) ; 

uimenu(PP, ’Label’, ’Grid’, ’Callback’, ’ Grid_stat ’ ) ; 
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77 uimenu (PP, ’Label*, ’Clear Graph’, ’Callback’, ’cla’); 

78 

79 User_ input = uimenu( ’Label’ , ’User Input’); 

80 uimenu (User.input , ’ Label ’ , ’ iperf numeric ’ , ’ Callback ’ , ’ get_iperf _data ’ ) ; 

81 

82 help = uimenu ( ’Label’ , ’Help’); 

83 uimenu (help, ’Label’ , ’Help’ , ’Callback’, ’help_Doc’); 

84 

85 * 

86 

87 '/.Read file name and path 

88 [handles. filename handles . path] = uigetfile(’*.out’) ; 

89 %Save file name for iperf call 

90 f i=handles . f ilename ; pa-handles . path ; 

91 save iperf.loc fi pa 

92 %Shov Path on YJ figure 

93 folder = getf older (handles. path); 

94 set (handles . f igure 1 , ’ Name ’ , f older) ; 

95 ^Define file location 

96 file_loc = [handles. path handles. f ilename] ; 

97 %========================================== 

98 %Show file name in filename box 

99 set (handles, file, ’String’ .handles, filename) ; 

100 

101 h = waitbar(0, ’Please wait...’); 

102 

103 '/, /(Import the data as Cell Array 

104 s = importdata(file_loc, ’\t’) ; 

105 

106 waitbar(10/100) ; 

107 

108 '/.If structure then convert to Cell Array 

109 if isstruct(s) 

110 s = struct2cell(s) 

111 end 

112 % /(Discretize the Cell Array by delimiter 

113 handles. data = for_data(s); 

114 

115 waitbar(.5) '/.update vaitbar 

116 

117 %Read labels 

118 [labels header_index handles. ref] = ... 

119 cell_find (handles. data, [{’SSID’} {’Multipath Data’) {’Data Rate’) ... 

120 {’RSSI’) {’Latitude (Dec. Deg. )’) {’Longitude (Dec. Deg.)’) ... 

121 {’Altitude (FT)’}]); 
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122 

123 

124 handles = col_f ind(handles) ; 

125 

126 7 

127 7 handles. ref contain columns location of 

128 % cell array specified above 

129 7 

130 

131 7SSID he ader_ index has ssid column number in first element 

132 vaitbar ( . 75) 

133 

134 TFind the rows of the rows of each ssid 

135 [ssid_names handles . ssid_row] = ssid_findnew (handles .data, handles. col_ssid) ; 

136 

137 7, Display SSid names is ssid_name_ array pop-up menu 

138 ssid_names = cellarray(ssid„names) ; 

139 set (handles . ssid_name_array , * String , , ... 

140 ssid_names) ; 

141 

142 %Get first ssid name displayed in handles. ssid_name_array 

143 handles « ssid_name_array_Callback (handles, ssid^name^array, eventdata, handles); 

144 7Get first Time Point displayed in handles . tsmpl 

145 handles = dwn_samp_Callback (handles. dvn_samp, event data, handles); 

146 

147 plot _Callback (hObject , eventdata, handles) 

148 waitbar(l) 

149 close (h) ; 

150 save handles 

151 

152 7===GUI STUFF= == ”============= ==” =================~===== s == ; 

153 7. I I I I I I I I I I I I I I i 

154 7.\/ V V \/ \/ \/ \/ \/ V \/ V \/ \/ \/ \/ 

155 % Choose default command line output for YJ_workl 

156 handles . output = hObject; 

157 % Update handles structure 

158 guidata (hObject, handles); 

159 7 UIWAIT makes YJLworkl wait for user response (see UIRESUME) 

160 7 uiwait (handles . f igurel) ; 

161 
162 

163 7 Outputs from this function are returned to the command line. 

164 function varargout = YJ_vorkl_OutputFcn (hObject, eventdata, handles) 

165 7 varargout cell array for returning output args (see VARARGOUT) ; 

166 7 hObject handle to figure 
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167 X eventdata reserved - to be defined in a future version of MATLAB 

168 % handles structure with handles and user data (see GUIDATA) 

169 

170 % Get default command line output from handles structure 

171 varargout{l> = handles . output ; 

172 

173 

174 % PUSH BUTTONS 

175 

177 

178 X Executes on button press in browse. 

179 function browse_Callback (hObject, eventdata, handles) 

180 clc; 

181 load handles 

182 7,Create Menu Options 

183 X 

184 % f * uimenu (’Label’ , ’Plot Protperties 9 ) ; 

185 X uimenu (f , * Label 7 , * New Figure 7 , 7 Callback 7 , 7 f igure ’ ) ; 

186 X uimenu (f , * Label * , 1 Save 7 , ’Callback’ , 9 save*) ; 

187 X uimenu (f , ’ Label * , * Quit * , ’ Callback ’ , ’ exit 7 , . . . 

188 X ’ Separator 7 , * on 7 , 7 Accelerator 7 , 7 Q 7 ) ; 

189 X 

190 file = handles. filename; pa = handles. path ; 

191 X save wksp hObject file pa 

192 X clear 

193 X load wksp 

194 X hObject handle to browse (see GCBD) 

195 X eventdata reserved - to be defined in a future version of MATLAB 

196 X handles structure with handles and user data (see GUIDATA) 

197 Xset (handles . ssid_name_array , 7 String 7 , ) ; 

198 XXWWWWW DATA DISPLAY PROGRAM /////////////////// 

199 

200 
201 

202 XRead file name and path 

203 [handles . f ilename handles . path] = uigetf ile ( , * . out 7 ) ; 

204 

205 

206 if handles. path ~= 0 

207 

208 XShow Path on YJ figure 

209 folder = getf older (handles. path) ; 

210 set (handles . f igure 1 , 7 Name 7 , f older) ; 

211 XDefine file location 
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212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 


file_loc = [handles . path handles. f ilename] ; 

‘/.Save file name for iperf call 
f i=handles . filename ; pa=handles . path ; 

save iperf_loc fi pa 

7. = ================================================================ 

'/.Show file name in filename box 

set (handles, file, ’String’ .handles, filename) ; 

h = waitbar(0, ’Please wait 

7 . ’/.Import the data as Cell Array 
s = importdata(file_loc, ’\t’) ; 

waitbar (10/100) ; 

7 . If structure then convert to Cell Array 
if isstruct(s) 

s = struct2cell(s) ; 

end 

’/. /(Discretize the Cell Array by delimiter 
handles. data = f or_data(s) ; 

waitbar (.5) '/.updata waitbar 

/(Read labels 

[labels header_index handles . ref] = ... 

cell_find (handles. data, [{’SSID’} {’Multipath Data’} {’Data Rate’} ... 

{’RSSI’} {’Latitude (Dec. Deg.)’} {’Longitude (Dec. Deg.)’} 

{’Altitude (FT)’}]); 

handles = col_find (handles) ; 

7. 

% handles. ref contain columns location of 

'/. cell array specified above 

•/. 


'ASS ID header_index has ssid column number in first element 
waitbar (.75) 

/(Find the rows of the rows of each ssid 

[ssid_names handles . ssid_row] = ssid_f indnew (handles .data, handles. col_ssid) ; 
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257 V 9 Display SSid names is ssid_name_array pop-up menu 

258 ssid_names = cellarray(ssid_names) ; 

259 set (handles. ssid_name_arr ay, ’String 1 , ... 

260 ssid_names) ; 

261 

262 %Get first ssid name displayed in handles . ssid_name_array 

263 handles * ssid_name_array_Callback(handles.ssid_name_array, eventdata, handles) 

264 /.Get first Time Point displayed in handies. tsmpl 

265 handles = dwn_samp_Callback (handles. dwn_samp, eventdata, handles); 

266 

267 waitbar(l) 

268 close (h) ; 

269 else 

270 handles. filename - file; 

271 handles. path = pa; 

272 end 

273 save handles ; 

274 % 

275 7.=-=====================— =========-======================= 

276 7A\\\\\\\\\\\\\\\\\\\\ Other plotting functions//////////////// 

277 % Executes on button press in surf. 

278 function surf ..Callback (hObject , eventdata, handles) 

279 % hObject handle to surf (see GCBO) 

280 % eventdata reserved - to be defined in a future version of MATLAB 

281 handles structure with handles and user data (see GUIDATA) 

282 [r c] = size (handles. current_data) ; 

283 if c > 1 

284 x = [-2:0.25:3.25] ; 

285 y = [1 :r] ; 

286 surf (x,y, handles. current _data(: , handles. col_in: handles . col_fn) ) ; 

287 xlabel ( 1 Chips 1 ) ; 

288 ylabel ( ’Relative time ’ ) ; 

289 

290 end 

291 

292 % Executes on button press in contour. 

293 function contour^Callback (hObject , eventdata, handles) 

294 % hObject handle to contour (see GCBO) 

295 7 eventdata reserved - to be defined in a future version of MATLAB 

296 % handles structure with handles and user data (see GUIDATA) 

297 [r c] = size (handles, current _data) ; 

298 if c > 1 

299 x = [-2:0.25:3.25] ; 

300 y - [l:r] ; 

301 contour (x,y, handles . current _data( : , handles. col_in: handles. col_fn)) ; 
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302 xlabel ( ' Chips ' ) ; 

303 ylabel ( 'Relative time ' ) ; 

304 end 

305 

306 

307 

308 % Executes on button press in mesh. 

309 function mesh_CaIIback (hObject, eventdata, handles) 

310 % hObject handle to mesh (see GCB0) 

311 % eventdata reserved - to be defined in a future version of MATLAB 

312 % handles structure with handles and user data (see GUIDATA) 

313 [r c] = size (handles. current_data) ; 

314 if c > 1 

315 x * [-2:0.25:3.25] ; 

316 y = [l:r] ; 

317 mesh(x, y, handles. current _data(: , handles. col_in: handles. col_fn)) ; 

318 xlabel ('Chips'); 

319 ylabel ('Relative time') ; 

320 end 

321 

322 

323 % Executes on button press in plot. 

324 function plot_Callback (hObject, eventdata, handles) 

325 % hObject handle to plot (see GCB0) 

326 % eventdata reserved - to be defined in a future version of MATLAB 

327 % handles structure with handles and user data (see GUIDATA) 

328 clc 

329 

330 if handles .file 

331 

332 

333 ts = handles . tsmpl ; 

334 

335 [r c] = size (handles, current _data) ; 

336 

337 %If Delay Spread is Selected 

338 if c “ 1 

339 %N-point Moving average filter for Delay spread 

340 len = length (handles, current _dat a) 

341 npts - ask_npts (handles) 

342 num = ones (1, npts) /npts; 

343 mv_data = filter (num, 1 , handles . current_data) ; 

344 

345 

346 title ( [num2str (npts) ' MA filter']) 
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347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 


grid on; 

ylabel ( ’Nanoseconds ’ ) ; 
xlabel ( ’Relative Time ’ ) ; 

'/.Mouse select Property 

run mouse.data 
%If Multipath is Selected 
elseif c > 1 

h = waitbar(0, ’Please wait...’); 

'/(Get Status of Radio Button 

val = get (handles . dB , ’Value ’ ) ; 
x = [-2:0.25:3.25]; %Chip vector 

%Create down sample index 
dsi = l:r/ts:r; 

y = [round(dsi)] ; '/.Time Sample vector 

'/.Check if Radio buton 

if val = 1 %In dB 

z - 10*logi0 (handles. current _data(y, ... 

handles . col_in : handles . col_f n) ) ; 
waitbar(.6) 

else %In decimal 

z = handles. current _data(y, ... 

handles . col_in : handles . col_f n) ; 
waitbar ( . 6) 

end 

axes (handles . axes 1 ) 

mesh (x, y, z) ; title ( [num2str(ts) ’ samples’]) 

xlabel ( ’ Chips ’ ) ; 

ylabel ( ’Relative time’); 

waitbar (.9) 
close (h) ; 

else 

title (’Could not display data’); 

end 

end 

save handles 

XWWWWWWWWWW POPUP MENU ///////////////////////////// 

•/.wwwwwwwwwwwww/////////////////////////////////// 

7 .\\\\\\\\\\\\\\\\\\\\\\\\\\/////// //////////////////////////// 

7 , Executes during object creation, after setting all properties. 

function select _data_CreateFcn(hObject, eventdata, handles) 

7 , hObject handle to select_data (see GCB0) 
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392 % eventdata reserved - to be defined in a future version of MATLAB 

393 7, handles empty - handles not created until after all CreateFcns called 

394 

395 guidata(hObject , handles) ; 

396 

397 % Hint : popupmenu controls usually have a white background on Windows . 

398 7. See ISPC and COMPUTER. 

399 if ispc 

400 set (hOb j ect , ’ BackgroundColor ’ , * white ’ ) ; 

401 else 

402 set (hObject, ’ BackgroundColor ’ ,get(0, ’defaultUicontrolBackgroundColor’ )) ; 

403 end 

404 

405 

406 % Executes on selection change in select .data. 

407 function handles = select.data.Callback (hObject, eventdata, handles) 

408 7# hObject handle to select.data (see GCB0) 

409 7. eventdata reserved - to be defined in a future version of MATLAB 

410 % handles structure with handles and user data (see GUIDATA) 

411 

412 val = get (hObject, ’Value*) ; 

413 str = get (hOb j ect , ’ String ’ ) ; 

414 switch str{val>; 

415 case 7 Multipath 7 

416 [r c] = size (handles. current .data) 

417 if r > 10 

418 time.sample.data(handles) ; 

419 end 

420 handles . current. data = handles. MP; 

421 %Clear Delay Spread info 

422 set (handles, mean , 7 String 7 , [] ) ; 

423 set (handles . std.dev , 7 String ’,[]); 

424 case ’Delay Spread’ 

425 set (handles . dwn.samp , ’ String ’ , □ ) ; 

426 handles . current .data = delay.spread (handles. MP) ; 

427 7.Calc mean and standard deviation 

428 handles . DS.mean = mean (handles , current.dat a) 

429 set (handles, mean, ’String’ .handles. DS.mean) ; 

430 handles. DS.sigma = std (handles, current. data) 

431 set (handles . std.dev, ’String’ , handles .DS.sigma) ; 

432 end 

433 guidata(h0bj ect .handles) ; 

434 % Hints: contents = get (hOb j ect , ’ Str ing ’ ) returns select.dat a contents as cell array 

435 % contentslget (hObject, ’Value’)} returns selected item from select.data 

436 
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437 %==================SSID NAME== : ====================== 

438 %=^=“==“===============POPUP============================ 

439 •/=== s=s=s = = = == == === = ======= == ==== = ===================== 

440 

441 •/' Executes during object creation, after setting all properties. 

442 function ssid_name_array_CreateFcn (hObject, eventdata, handles) 

443 % hObject handle to ssid_name_array (see GCBO) 

444 */ 9 eventdata reserved - to be defined in a future version of MATLAB 

445 % handles empty - handles not created until after all CreateFcns called 

446 guidat a (hObject .handles) ; 

447 % Hint: edit controls usually have a white background on Windows. 

448 % See ISPC and COMPUTER. 

449 

450 if ispc 

451 set (hOb j ect , 3 BackgroundColor 3 , 3 white 3 ) ; 

452 else 

453 set (hOb j ect , 3 BackgroundColor 3 , get (0 , 3 def aultUicontrolBackgroundColor * ) ) ; 

454 end 

455 

456 ‘/.THIS IS THE USER SELECT SSID FUNCTION 

457 ^^^ ==== == = == = ^ == ^ ====== = ^^„^ = ^ ======= ^ ================ 

458 % I | I | | | | i | I | | | | | | | | | | 

459 X\/ \ f \/ \/ \/ \/ \/ \/ V \/ \/ \/ \/ \/ V \/ \/ \/ \/ \/ 

460 function handles = ssid_name_array_Callback (hObject, eventdata, handles) 

461 % hObject handle to ssid_name_array (see GCBO) 

462 % eventdata reserved - to be defined in a future version of MATLAB 

463 % handles structure with handles and user data (see GUIDATA) 

464 

465 val = get (hOb j e ct , 3 Value * ) ; 

466 str = get (hObject, 3 String 7 ); 

467 if length (str) <2 k val >1 

468 val = 1; 

469 end 

470 handles . ssid_name = str (val); 

471 

472 /(handles = col_find(handles) 

473 

474 7, Use this code until fixed ssid_rov from function 

475 */ 9 ssid_f indnew 

476 len = length (handles, ss id_r ow) ; 

477 if mod (len, 2) 

478 handles. row = handles. ssid_row(2*val: (2*val+l)) ; 

479 else 

480 handles. row = handles. ssid_row(2*val-l : (2*val)) ; 

481 end 
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482 

483 >CSet the time sample data select options 

484 

485 handles = time .sample .data (handles) ; 

486 

487 

488 

489 


490 % Multipath data 

491 col = [handles . col.in handles . col.fn] ; 

492 handles. MP = dat a.f etch (handles, row, col, handles .data) ; 

493 *L — Set current data 

494 if get (handles, select _dat a, ’Value*) == 2 

495 handles. current.data = delay.spread (handles .MP) ; 

496 ftCalc mean and standard deviation 

497 handles. DS.mean = mean (handles, current .data) ; 

498 set (handles . mean , * String ’ , handles . DS .mean) ; 

499 handles . DS.sigma = std (handles . current .data) ; 

500 set (handles . std.dev, ’String* , handles .DS.sigma) ; 

501 else 

502 handles . current .data = handles. MP; 


503 end 

504 

505 

506 ^handles . ssidname = get (handles . ssid.name.array , * String * ) ; 

507 guidata(h0bject, handles) ; 

508 % Hints: get (hObject, 'String* ) returns contents of ssid.name.array as text 

509 % str2double (get (hObject , ’String’)) returns contents of ssid.name.array as a double 

510 

511 ^======= — ======== === ==fiLE NAME=== :: = s ========= 

512 X— — 

513 •/========================================= 

514 

515 */. Executes during object creation, after setting all properties. 

516 function file.CreateFcn (hObject , eventdata, handles) 

517 '/. hObject handle to file (see GCB0) 

518 % eventdata reserved - to be defined in a future version of MATLAB 

519 % handles empty - handles not created until after all CreateFcns called 

520 guidat a (hObject , handles) 

521 

522 y % Hint: edit controls usually have a white background on Windows. 

523 % See ISPC and COMPUTER. 

524 if ispc 

525 set (hObject, ’BackgroundColor ’ , ’white’); 

526 else 
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527 set (hObject, 'BackgroundColor' ,get(0, 'def aultUicontrolBackgroundColor ' ) ) ; 

528 end 

529 

530 

531 

532 function file.Callback (hObject, eventdata, bandies) 

533 X hObject handle to file (see GCBO) 

534 % eventdata reserved - to be defined in a future version of MATLAB 

535 X handles structure with handles and user data (see GUIDATA) 

536 ha n d l es .file * get (handles . file , * String ' ) ; 

537 guidata(hObject, handles) 

538 X Hints: get (hObject , ' String’) returns contents of file as text 

539 X str2double (get (hObject , ’String ' ) ) returns contents of file as a double 

540 

541 

542 %====— =======—=EXTERNAL FIGURE=-===-=-=-=====-======-======== 

543 •/•========— 

544 X Executes on button press in ext.fig. 

545 function ext _fig.Callback (hObject, eventdata, handles) 

546 7, hObject handle to ext.fig (see GCBO) 

547 eventdata reserved - to be defined in a future version of MATLAB 

548 X handles structure with handles and user data (see GUIDATA) 

549 name = figure 

550 file = handles .filename; 

551 ssid = handles. ssid.name{l>; 

552 set (name, 'Name* , handles, filename) 

553 set (name, 'NumberTitle ' , 'off ') 

554 

555 val = get (handles. dwn.samp, 'Value') ; 

556 str = get (handles . dwn_samp , 'String'); 

557 ^Convert to integer 

558 

559 ts = handles . tsmpl ; 

560 

561 [r c] * size (handles, current .data) ; 

562 if c — 1 

563 plot (handles . current .data) 

564 hold on; 

565 plot (handles. DS.mean*ones(l, length (handles, current .data)) , 'g') ; 

566 hold off; 

567 title ([' Delay Spread of ' f ile(l: length (f ile)-9)] ) ; 

568 elseif c > 1 

569 %Get Status of Radio Button 

570 val = get (handles. dB, 'Value') 

571 x = [-2:0.25:3.25]; %Chip vector 
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573 

574 

575 

576 

577 

578 

579 

580 

581 

582 

583 

584 

585 

586 

587 

588 

589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 

600 
601 
602 

603 

604 

605 

606 

607 

608 

609 

610 
611 
612 

613 

614 

615 

616 


‘/.Create down sample index 
dsi = l:r/ts:r; 

y = [round(dsi)] ; ‘/.Time Sample vector 

‘/.Check if Radio buton 

if Veil = 1 ‘/.In dB 

z = 10*logl0(handles. current _data(y, ... 
handles . col_in: handles . col_fn) ) ; 
else %In decimal 

z = handles. current _dat a (y, ... 

handles . col_in : handles . col_fn) ; 

end 

mesh(x,y,z) ; 

title([’ Power Delay Profile of ’ f ile(l :length(f ile)-9)] ) ; 

else 

title ([’Could not display data for ’ ssid]); 

end 

xlabel ( ’ Chips ’ ) ; 

ylabelC ’Relative time’); 

zlabel( [’Relative Correlation of ’ ssid]); 

% Executes during object creation, after setting all properties. 

function mean_CreateFcn(hObject, eventdata, handles) 

V, hObject handle to mean (see GCB0) 

'/, eventdata reserved - to be defined in a future version of MATLAB 

handles empty - handles not created until after all CreateFcns called 

’/. Hint: edit controls usually have a white background on Windows. 

7. See ISPC and COMPUTER, 

if ispc 

set (hOb j ect , ’ BackgroundColor ’ , ’ white ’ ) ; 

else 

set(h0bject, ’BackgroundColor’ ,get(0, ’defaultUicontrolBackgroundColor’)) ; 

end 


function mean_ Callback (hObject, eventdata, handles) 
y, hObject handle to mean (see GCB0) 

eventdata reserved - to be defined in a future version of MATLAB 
'/, handles structure with handles and user data (see GUIDATA) 

*/, Hints: get (hObject, ’String’) returns contents of mean as text 

'/. str2double (get (hObject, ’String’)) returns contents of mean as a double 


127 


F CODELISTINGS 


617 

618 X Executes during object creation, after setting all properties. 

619 function std_dev_CreateFcn (hObject, eventdata, handles) 

620 7, hObject handle to std^dev (see GCBO) 

621 X eventdata reserved - to be defined in a future version of MATLAB 

622 X handles empty - handles not created until after all CreateFcns called 

623 

624 % Hint: edit controls usually have a white background on Windows. 

625 X See ISPC and COMPUTER. 

626 if ispc 

627 set (hObject, 'BackgroundColor* , ’white*) ; 

628 else 

629 set (hObject , ’BackgroundColor * ,get (0, *def aultUicontrolBackgroundColor ’ ) ) ; 

630 end 

631 

632 

633 

634 function std_dev_Callback (hObject, eventdata, handles) 

635 X hObject handle to std_dev (see GCBO) 

636 X eventdata reserved - to be defined in a future version of MATLAB 

637 X handles structure with handles and user data (see GUIDATA) 

638 

639 X Hints: get (hObject , ’ String 1 ) returns contents of std_dev as text 

640 X str2double (get (hObject, ’String’)) returns contents of std_dev as a double 

641 

642 

643 X Executes on button press in dB. 

644 function dB_Callback(hObject, eventdata, handles) 

645 X hObject handle to dB (see GCBO) 

646 X eventdata reserved - to be defined in a future version of MATLAB 

647 X handles structure with handles and user data (see GUIDATA) 

648 val = get (handles. select__data, * Value' ) ; 

649 if val == 1 

650 plot ^Callback (hObject, eventdata, handles) 

651 end 

652 X Hint: get (hObject, y Value’ ) returns toggle state of dB 

653 

654 

655 X Executes during object creation, after setting all properties. 

656 function dwn_samp_CreateFcn (hObject, eventdata, handles) 

657 X hObject handle to dwn_samp (see GCBO) 

658 X eventdata reserved - to be defined in a future version of MATLAB 

659 X handles empty - handles not created until after all CreateFcns called 

660 

661 X Hint: popupmenu controls usually have a white background on Windows. 
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662 */. See ISPC and COMPUTER. 

663 if ispc 

664 set(hObject , 'BackgroundColor' , ' white 1 ) ; 

665 else 

666 set (hOb j ect , ' BackgroundColor ' , get (0 , 9 def aultUicontrolBackgroundColor 9 ) ) ; 

667 end 

668 

669 

670 % Executes on selection change in dwn.samp. 

671 function handles = dwn_samp_Callback(hObject , eventdata, handles) 

672 % hObject handle to dwn_samp (see GCBO) 

673 eventdata reserved - to be defined in a future version of MATLAB 

674 % handles structure with handles and user data (see GUIDATA) 

675 % 

676 yCThis functin get the displayed value on Time Points 

677 % handles . tsmpl 

678 7. 

679 val * get (hOb j ect , 9 Value * ) ; 

680 str = get (hOb j ect , 9 String * ) ; 

681 

682 opt « 1; 

683 

684 switch str {val}; 

685 case ’user' 

686 %Ask user to input number of samples desired 

687 tot = handles. row (2) -handles. row(l ) ; 

688 answer = tot+1; 

689 while answer <= 1 I I answer >= tot 

690 

691 prompt - {[* Enter number between 1 : 9 num2str(tot)] }; 

692 dlg_title = 'Time Points'; 

693 num_lines= 1; 

694 answer = str2num( char (input dig (prompt ,dlg_title,num_lines))) ; 

695 

696 if answer <= 1 I I answer >= tot 

697 errordlg([ 'Please enter number between 1 and ' 

698 num2str (tot)] ) ; 

699 waitf or (gcf ) ; 

700 end 

701 end 

702 

703 handles. tsmpl * answer; 

704 opt = 0; 

705 end 

706 
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707 if opt == 1 

708 handles . tsmpl = str2num(char(str(val))) ; 

709 end 

710 y.plot_Callback(hObject , eventdata, handles) 

711 guidata (hOb j ect , handles) 

712 

713 

714 7, Hints: contents = get (hObject, ’String’) returns dwn_samp contents as cell array 

715 7, contents {get (hObject, ’Value’) > returns selected item from dun.samp 

716 

717 

IEEE 802 . 11 a Simulation Code 

simJERJLla.m 

1 % SIM_ER_11A.M 

2 5* This is the top-level script for IEEE 802.11a PHY layer simulations, 

3 5C using the mWLAN Toolbox by CommAccess Technologies. 

4 */. 

5 % To run a set of simulations, just modify the ’simName’ variable below, 

6 X so that ’RunParams_<simName> .m’ is the name of the script file 

7 X which sets all the simulation parameters for all the simulation runs 

8 Ji to be done in this set. 

9 % 

10 X Some terminology to make sense of the nested iterations: 

11 7, Calling ’sim_ER_lla.m’ once executes a ’simulation set’. 

12 % A simulation set entails N ’simulation runs. ’ 

13 X Each simulation run simulates a bunch of packets, one by one. 

14 X The number of packets in a run is upper-bounded by the ’MaxPkts’ 

15 X parameter in the ’RunParams_X.m’ file. But if ’enough’ errors 

16 % occur earlier, then the run will stop before sim’ing all MaxPkts 

17 % packets. ’Enough’ is defined by the ’MinErrPkts’ parameter, 

18 % and corresponds to the number of pkts with errors, and not just 

19 X the number of bit errors, so that one pkt with 1000 bit errors 

20 % will not end the run. 

21 51 The ’vsEbNo’ flag, set in the ’RunParams_X.m’ file, determines 

22 % how noise is calculated. If ’vsEbNo’, then noise is normalized 

23 5* to produce BER-vs-EbNo curves. If not, then noise is based on 

24 5* an absolute-scale estimate of the noise on Mars, and the signal 

25 5i power is the absolute power ’txPwr’ (Watts) from ’RunParams_X.m’ . 

26 51 The ’skipDemod’ flag allows faster simulation, if you only want 

27 5* Pkt Error Rate, based on correct reception of the pkt header. 
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28 X If ’skipDemod’, then only the SIGNAL symbol in the header is 

29 X demodnlat ed/decoded. The data is ignored. This, together with 

30 7, specifying a very short ’pktLen’, results in much better sim 

31 X times, when only PER is desired. 

32 7. 

33 % This script produces a (very large) diary file ’simName.txt ’ . 

34 7. After each run in the set, a result table will be dumped into the 

35 % diary file. The format of the result table is described in the 

36 X comments for the dumpResultTableO function. The result tables 

37 X are cumulative throughout the diary file, so that the table at 

38 X the very end of the diary file summarizes all the results for 

39 X all the runs in the entire diary file. This final table should 

40 X be extracted into a managably-sized ’ Result s^simName.m’ file 

41 X (containing only the one table), and a ’SumAndPlot_simName .m’ 

42 X or ’SumAndTab.simName.m’ script can be customized to produce 

43 X whatever table/plot is desired. 

44 X 

45 X For discussion about the calculation of noise power on Mars, read 

46 X ’abacusl/nasa_scp/Simulations/80211a/PowerDistanceEtc/pwr.pdf * . 

47 X and the subsection 

48 X , 4. PERFORMANCE AT DIFFERENT SITES’ 

49 X ’Performance versus distance between the transmitter and the receiver’ 

50 X in the 2005 IEEE Aero Conference paper: 

51 X 9 abacus l/nasa_scp/Documents/IEEE_Aero_Conf _2005/Perf EvalllFinal/docf inal . pdf ’ . 

52 X 

53 X 2004, Gaylon R Lovelace, New Mexico State University 

54 X 

55 

56 clear all, format compact, close all; 

57 randn(’state’ ,sum(100*clock)) ; 

58 

59 simName = ’pwrCheck’; 

60 

61 diary ( [ simName ’.txt’ ] ); datestr(now) 

62 

63 skipDemod = 0; X default to normal demod 

64 vsEbNo = 0; X default to SNR determined from ICS Telecom 

65 eval( [ ’RunParams_’ simName ] ); X sets an Nx9 cell array with all the simulation paramet* 

66 runParamNums = cell2mat (runParams ( : ,1:8)) ; 

67 nRuns = size( runParams, 1 ); 

68 

69 magicNums = setMagicNums; X sets a bunch of 802.11a specific magic numbers 

70 

71 pktCums = zeros (nRuns, 1) ; bitCums = pktCums; errBitCums = pktCums; nErrs = 0; 

72 berrPktCums = pktCums; perrPktCums = pktCums; serrPktCums = pktCums; aerrPktCums = pktCui 
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rmsDelays = pktCums ; rxPowers = pktCums; estEbNoDBs = pktCums; 


for runNdx = [ 1 : nRuns ] 


Mbps = 
tx_rateMode = 
tx_codeRate = 
tx.NJBPSC = 
tx_N_CBPS = 
EbNoDB 

distance = 
tx_PvrScalar= 
tx_nBytes = 
maxPkts = 
minErrPkts = 
chnl Type = 
pdpFile = 


cell2mat ( runParams (runNdx, 2) ) ; 

Mbps2rateMode( Mbps ); 

magicNums .RATE_dependent (magicNums . xMODE(tx_rateMode+l) ,3) ; 

magicNums . RATK_dependent (magicNums . xMODE(tx_rateMode+l ) , 4) ; 

magicNums .RATE_dependent (magicNums .xMODE(tx_rateMode+l) ,5) ; 

cell2mat( runParams (runNdx, 3) ); 

cell2mat( runParams (runNdx, 4) ); 

cell2mat( runParams (runNdx, 5) ); 

cell2mat( runParams (runNdx, 6) ); 

cell2mat ( runParams (runNdx, 7) ) ; 

cell2mat( runParams (runNdx, 8) ); 

char( runParams (runNdx, 9) ); 

char( runParams (runNdx, 10) ); 


[ dateStr, V_m ] = readPDP( [ J . . \PDP\ ’ pdpFile ’ .csv’ ), chnlType, 10e-9 ); 
V_m = downsamplePDP ( V_m, 5, 500 ) .* sqrt (tx_PwrScalar) ; 


% Constants for Mars Propagation 
c = 3e8; 

lambda = c / 5.25e9; 7. IEEE 802.11a ==> 5.25 GHz 

Gr = 1.0; 7, Isotropic Rx Antenna 

kBoltz = 1.38e-23; 

Teq = 1560; % Equivalent Noise Temp Kelvins 

Convert from E-Field (V/m) to Rxr Input Voltage RMS: 
rxVrms = V_m/10 * ( sqrt(Gr) * lambda / (pi*sqrt (480)) ); 
if (vsEbNo) 

rxVrms » rxVrms’ / sqrt( abs( rxVrms ’ * rxVrms ) ); % normalize 

end 

7, Ensure at least 64 samples in h(t) 

if (length(rxVrms)<64); rxVrms = [ rxVrms ( : ) ’ , zeros(l,64) ]; rxVrms = rxVrms(l:6‘ 
if (vsEbNo) 

*/. tx_wave is normd so Es = 1.000 — > reqd_No = 1 / EsNoLinear; 

EbNoLinear = 10 .“ ( EbNoDB / 10 ); 

awgnSigma = sqrt ( 1 / ( EbNoLinear * tx_N_CBPS * tx_codeRate ) ) ; 

awgnSigma = awgnSigma * sqrt (80); 7, because is done below, for tx_wave 

else 

V, Compute noise power on Mars 
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awgnSigma = sqrt ( kBoltz * Teq * 20e6 / 2 ) ; 

end 

rmsDelays (runNdx) = rmsDelayCalc( rxVrms, l/20e6, 1 ); 
rxPowers (runNdx) = ( rxVrms ( : ) * * rxVrms ( : ) ) ; 
rxVrms = rxVrms ( : ) ; 

leadPwr = ( rxVrmsCl: 16) ’ * rxVrms(l:16) ); 
tailPwr = ( rxVrms (17: end) '• * rxVrms (17 rend) ); 
totalPwr = ( rxVrms(:)’ * rxVrms(r) ); 
disp( ’leadPwr tailPwr totalPwr’ ); 
disp( [ leadPwr tailPwr totalPwr ] ); 

for repNdx = [ 1 : maxPkts ] 

disp( sprintf ( ’runNdx=*/.d; AntHt=7.gm; EbNoDB=’/,g; Mbps=7.d; Rep#7.d; 7.s ’ , run] 

tx_PSDU = bitgen( 8*tx_nBytes ) ; 

tx.wave = plcp_lla( tx.rateMode, tx.PSDU ); 


7 . Renormalize tx_wave to SampleEnergy=l . 0 , rather than 0FDMSymbolEnergy=l . 0 
tx_wave = tx_wave * sqrt (80) ; 
h_t = rxVrms; 


7 . AWGN Channel is just one tap. 

7Ji_t = zeros (1,64); h_t(l) “ 1.0; 

7 , Fading Model: h(t) = h(t) * Rayleigh ampl * uniform phase 

h_t = ( sqrt(0.5) * [1 j] * randn (2, length (h_t)) ) .* h_t(:)’; 

rx_wave = conv( tx_wave, h_t ); 

5lrx_wave = tx_wave; 

7 , For debug purposes, estimate Eb/No 

estEbNoLinear = ( var(rx_wave) / awgnSigma“2 ) / ( tx_N_CBPS * tx.codeRate / 80 ); 
awgnoise = awgnSigma * sqrt (1/2) * [1 j] * ( randn ( 2 , length (rx_wave ) ) ); 
rx_wave = rx_wave + awgnoise; 


h_t_hat = magicNums . LSE_mtrx * [rx_wave(l: (magicNums.preambLen)) . ’] ; 

7. Equalization at rxr will require transform of whichever impulse response : 
H_k_hat = fft( h_t_hat(l :64) ); 
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'/. Up to this point, everything is in MKS units (Volts, meters. Watts, etc) 

% The receiver simulation below has no gain control. Except for premult ipl; 
V. tx_wave by sqrt(80) (to make it apples-to-apples with the noise), everyth: 
% done above is un-done by H_k_hat in the receiver. Scale the rx’d signal « 
% together to ’set the receiver preamp gain’ . 

rx_wave = rx_wave / sqrt(80); 7, Only for ’real power’ multipath via ICS' 

signalStart = 320 +16+1; 

[ rx_rateMode , rx_nSyms, rx_nBytes, signalErrFlags, signalErrMsg ] = ... 

processSIGNAL( rx_wave (signalStart : (signalStart +63) ) , magicNums, H_k_hai 

rx_rateMode = tx_rateMode; ’/, for BER purposes, ignore the SIGNAL symbol contents 
rx_nBytes = tx_nBytes ; 

if (~skipDemod) 

[ serviceBits, dataBits ] = processData( rx_wave, rx_rateMode, rx.nBytes, m; 
end 

if CskipDemod) , nErrs = sum( dataBits ~= tx_PSDU ); end; 
pktCums (runNdx) = pktCums (runNdx) + 1; 

bitCums (runNdx) = bitCums (runNdx) + 8*tx_nBytes ; 

errBitCums (runNdx) = errBitCums (runNdx) + nErrs; 

berrPktCums (runNdx) = berr PktCums (runNdx) + (nErrs >0) ; 

perrPktCums (runNdx) = perrPktCums (runNdx) + bitand( signalErrFlags, 1 ); 
serrPktCums (runNdx) = s err PktCums (runNdx) + ( signalErrFlags > 1 ) ; 
aerr PktCums (runNdx) = aerr PktCums (runNdx) + ( (nErrs>0) II (signalErrFlags: 
estEbNoDBs (runNdx) = estEbNoDBs (runNdx) + estEbNoLinear; '/. not in dBs until late] 

if( length(signalErrMsg) > 0 ) 

disp( [ ’ ’ signalErrMsg ] ) 

end 

statStr = sprintf ( ’ This pkt : 7,d/7.d bits . ’ , nErrs , 8*tx_nBytes ) ; 

statStr • [ statStr sprintf ( ’Cumulative: 7.d/7,d errPkts k 7.d/7,d bits.’, aen 
disp( statStr ); 

if ( aerr PktCums (runNdx) >= minErrPkts ) 

'/, Have plenty of errors for significance. Skip remaining packets ii 
break; 

end 

end repNdx 
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208 estEbNoDBs (runNdx) = 10 * loglO ( estEbNoDBs (runNdx) / pktCums (runNdx) ) ; 

209 

210 ‘/.save ( sprintf( ,, /,s_7s_curve702d.mat ’ , simName , chnlType, runParamNums (runNdx , 1) ) 

211 

212 resultMatrix = [ errBitCums, be rr PktCums , perrPktCums , serrPktCums, aerrPktCums, . . 

213 bitCums, pktCums, rmsDelays, rxPowers, estEbNoDBs ]; 

214 resultTable = dumpResultTable( resultMatrix, runParams ); 

215 disp( resultTable ); 

216 

217 end '/, runNdx 

218 

219 diary off ; 

220 


Mbps2rateMode.m 

1 function rateMode = Mbps2rateMode( Mbps ) 

2 7MBPS2RATEM0DE Convert 802.11a Mbit/s to a mWLAN rateMode value 

3 7. 48 54 12 18 24 36 6 9 Mb/s 

4 7 01234567 rateMode 

5 7. 

6 7 Subfunction for SIM_ER_11A.M 

7 7 

8 7 2004, Gaylon R Lovelace, New Mexico State University 

9 7 
10 

11 table = [ 

12 48 54 12 18 24 36 6 9 7 Mb/s 

13 0 1 2 3 4 5 6 7 7 rateMode 

14 ]; 

15 ndx = find( Mbps = table (1,:) ); 

16 rateMode = table (2, ndx) ; 


setMagicNums.m 

1 function [magicNums] = setMagicNumsO 

2 7SETMAGICNUMS Returns a struct full of magic numbers for use in 802.11a Rcvr 

3 7. Subfunction for SIM_ER_11A.M 

4 7 

5 7 2004, Gaylon R Lovelace, New Mexico State University 

6 7 
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7 

8 load( ’ tx_wave_preamble ’ ); tx_wave_preamble = tx_wave_preamble ( : ) ; 

9 chill Len = 64; preambLen = length (tx_wave_preamble) ; 

10 A = convmtx( tx_wave_preamble , chnlLen ); 

11 A = A( [1 : preambLen] , : ) ; 

12 '/, don’t use the zero-prefixed rows at bottom; they overlap with the SIGNAL symbol. 

13 


14 

magicNums = struct ( . 

. 


15 

, pnVector > , [ ... 



16 

1111-1-1 

-1 1 -1 -1 

-1-1 1 .. 

17 

h* 

i 

1 

1 

1-111 

-1 1 1 .. 

18 

1111-11 

11-11 

1-1-1 

19 

1 1 l-l 1-1 

-1 -1 1 -1 

1-1-1 

20 

1-1-1 1 1 1 

1 1-1-1 

1 1 -1 .. 

21 

-1 1-1 1-1 1 

1 -1 -1 -1 

1 1 -1 . . 

22 

-1 -1 -1 1 -1 -1 

1-111 

1 1 -1 . . 

23 

1-1 1-1 1-1 

-1 -1 -1 -1 

1-1 1 . . 

24 

1-1 1-1 1 1 

1-1-1 1 

-1 -1-1 . . 

25 

1 1 1 -1 -1 -l 

-1 -1 -1 -1] 

i • * • 

26 

’u64’ , [0:1:63], ... 



27 

’ i_symbol ’ , [8 22 44 58] , ... 


28 

’i_data’ , [1:5 7:19 21:26 28:33 35 

:47 49:53] , 

29 

’v_symbol’, [1-111], ... 


30 

’vl’, [28:1:53], ... 



31 

’v2’ , [2:1:27], ... 



32 

’v3’ , [1:1:26], ... 



33 

’v4’ , [39:1:64], ... 



34 

’RATE_dependent’ , [ . 

. 


35 

6 11/21 48 24 



36 

9 1 3/4 1 48 36 



37 

12 2 1/2 2 96 48 



38 

18 2 3/4 2 96 72 



39 

24 3 1/2 4 192 96 



40 

36 3 3/4 4 192 144; . . . 


41 

48 4 2/3 6 288 192; ... 


42 

54 4 3/4 6 288 216] .... 


43 

’xMODE’ , [783456 

12], ... 


44 

’am’, [l./sqrt([l 1 2 2 10 10 42 42])], ... 


45 ’preambLen’ , preambLen, . . . 

46 ’ LSE_mtrx ’ , inv ( A ’ *A) *A ’ ... 

47 ); 

48 

49 ’D_hat’, [ zeros (1,53) ], ... 

50 ’S_hat’, [ zeros(l,48) ), ... 

51 % ’kpn’ , [1] , ... 


136 


F CODELISTINGS 


52 y. ’u80’ , [80], ... 

53 

54 

55 

56 %, ... 

57 y.’tx_vave_preamble’ , tx_wave_preamble 


readPDP.m 

1 function [ dateStr, V_m ] = readPDP( fileName, grepTxt, Ts ) 

2 XREADPDP Read a delay profile from an ICS Telecom CSV file 

3 % [ dateStr, V_m ] = readPDP( fileName, grepTxt , Ts ) 

4 % opens the named file, and reads the comma separated data from 

5 V, all lines which begin with the text string GREPTXT. 

6 '/, DATESTR returns the timestamp from the first line of the file. 

7 % V_M returns the E-Field Intensity profile, converted from dBuV/m 

8 % into linear units of Volts/meter. Zeros are inserted for all 

9 '/, missing time intervals, including those prior to the first 

10 y, tap in the data. The elements of V_M correspond to Times-of- 

11 % Arrival of [ 0 : Ts : maxTime ] where maxTime is the maximum ToA 

12 in the file data. 

13 y. 

14 'I, An error will result if any of the Times -of -Arrival given in 

15 the file are not an integer multiple of Ts. 

16 y. 

17 The CSV file should match the following six-column format: 

18 */. 5/10/2004 10:17,,,,, 

19 y. 

20 % Subscriber, #,X,Y,FS dBV/m, ToA s 

21 % rx20ml, 6, 102319. 2031, 92928, 102, 0.07 

22 % rx20ml, 6, 102319. 2031, 92928, 99, 0.1 

23 rx20ml, 6, 102319. 2031, 92928, 100, 0.11 

24 % rxlOOml, 4, 102319. 2031, 92855. 20313, 102, 0.31 

25 '/. rxlOOml, 4, 10231 9. 2031, 92855. 20313, 93, 0.32 

26 % rxlOOml, 4, 102319. 2031, 92855. 20313, 99, 0.33 

27 % rxlOOml, 4, 102319. 2031, 92855. 20313, 100, 0.35 

28 y. 

29 For this example, READPDP( fileName, ’rxlOOml’, 10e-9 ) will return: 

30 % DATESTR = ’5/10/2004 10:17’; 

31 */. V_M = 36x1 array = [ zeros(l,31), 10*(102/20), 10‘(93/20), 10“(99/20), 0, 10~(100/20) 

32 */. T0A_S = 36x1 array = [ 0 : 10e-9 : 350e-9 ] ’ ; 

33 y. 

34 y. 
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’/. Subfunction for SIM_ER_11A.M 

7 . 

% 2004, Gaylon R Lovelace, New Mexico State University 

*/. 

fid = f open( f ileName ) ; 
dateStr = fgetl( fid ); 
dateStr = dateStr (i: (end-5)) ; 

% fgetl( fid ); fgetl( fid ); % skip blank line b column headers 


state = 2; 
k = 0; 

data = repmat ( [ -1 -1 ], 8, 1 ); */, create ‘empty’ data array 

while ( state > 0 ) 

oneLine = fgetlC fid ); 
if ( oneLine == -1 ) 
state = 0; 

elseif ( isempty( oneLine ) ) 

% do nothing; skip blank lines 

else 


end 


end 


c = find( oneLine == ); 

txt = oneLine(l : (c(l)-l)) ; 
if ( strcmpi( txt, grepTxt ) ) 
k = k + 1; 
state = 1; 

data(k,l) = str2num( oneLine ((c(4)+l) : (c(5)-l)) ); 
data(k,2) = str2num( oneLine ((c (5) +1) rend) ); 
elseif ( state = 1 ) 

7.7, did have ‘matching’ lines, but now are past them 
state =0; 

end 


fclose( fid ); 


data = data(l:k,:); 


t_over_T = 1 + round( ( data(:,2) * le-6 ) / Ts ) ; 
V_m = zeros ( max(t_over_T) , 1 ); 

V_m(t_over_T) = 10. *( ( data(:,l) / 20 ) - 6 ) ; 


return 
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7,7, Example of how to use this f ctn 
lambda = 3e8 / 5.25e9; 
kBoltz = 1.38e23; 

Teq = 1560; 

tx_wave = plcp_lla( rateNode, PSDU ); 
tx_wave = tx_wave * sqrt(80); 
pdp = readPDP( fileName, tag, 10e-9 ); 
pdp = pdp (min (find (pdp) ): end) ; 

7, h(t) = pdp * Rayleigh ampl * uniform phase 
h_t = ( sqrt (0.5) * [ 1 j ] * randn (2, length (pdp) ) ) 
h_t = h_t * ( lambda / (pi*sqrt (480) ) ); 
rx_wave = conv( tx_wave, h_t ); 
awgnSigma = sqrt( kBoltz * Teq * 20e6 / 2 ); 
awgnoise = awgnSigma * sqrt(l/2) * [l;j] * ( randn(2,length(tx_wave)) ); 
rx_wave = rx_wave + awgnoise; 


7. IEEE 802.11a => 5.25 GHz 


7, Equivalent Noise Temp Kelvins 


'/. Read V/m profile from file. 
7, Discard leading zeros 


pdp; 


downsamplePDP.m 

1 function V_m = downsamplePDP( raw_V_m, M, N ) 

2 7.D0WNSAMPLEPDP Filter k downsample a delay profile from ICS Telecom 

3 7% V_m = readPDP( raw_V_m, M, N ) takes a PDP from readPDPO, and 

4 % lowpass filters it with an Nth order filter, and downsamples it 

5 % by a factor of M. V_m is normalized to have total energy equal 

6 7% to that of raw_V_m. 

7 7. 

8 7. Subfunction for SIM_ER_11A.M 

9 7. 

10 '/. 2004, Gaylon R Lovelace, New Mexico State University 

11 7. 

12 

13 myFilt = firl( N, 1/M ); 

14 

15 first * round( min( find( raw_V_m ) ) + (N/2) ); 

16 last = round( first + length( raw_V_m ) ); 

17 V_m = filter( myFilt, 1, [ raw_V_m; zeros(f irst+N/2, 1) ] ); 

18 

19 V_m = V_m(first:M:last) ; 

20 

21 V_m = V_m * sqrt ( ( raw_V_m * * raw_V_m ) / ( V_m ’ * V_m ) ) ; 

22 
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processData. m 

function [ serviceBits, dataBits ] = processDataC rx.wave, rx_rateMode, rx.nBytes, magicNumi 
‘iPROCESSDATA Process the 802.11a data symbols 

7 . RX_WAVE is the entire received baseband signal, including preamble k SIGNAL symbol. 

% The actual data symbols start at sample 401. 

7, 320 samples = skip over the PLCP Preamble (16us * 20MHz); 

7, 80 samples skips over the SIGNAL symbol (4.0us * 20MHz) 

% 1 = MATLAB index is one-based b we want first sample AFTER the SIGNAL symbol 

7 . MAGICNUMS = setMagicNumsO; 

% H_k = 64-sample FFT of multipath channel impulse response 

7t RX_RATEMQDE, RX_NBYTES should come from the SIGNAL symbol, or from knowledge 
7. of the transmitter. 

7. RX_RATEM0DE should be 0=48Mbps, 1=54, 2=12, 3=18, 4=24, 5=36, 6=6, 7=9Mbps 

7. These are just bits R1-R4 of the SIGNAL field 

7t magicNums .xMODE can sort these into ascending bitrate 

7 . 

% Subfunction for SIM_ER_11A.M 

7 . 

7t 2004, Gaylon R Lovelace, New Mexico State University 

7 . 


rx_codeRate = magicNums. RATE_dependent (magicNums. xMODE (rx_rateMode+l) ,3) 
rx_N_BPSC = magicNums. RATE_dependent (magicNums. xM0DE(rx_rateMode+l) ,4) 

rx_N_CBPS = magicNums. RATE_dependent (magicNums. xM0DE(rx_rateMode+l) ,5) 

rx_AMref = magicNums. am( magicNums. xM0DE(rx_rateMode+l) ) 

rx_nSym = ceil ( ( 16 + 8*rx_nBytes + 6 ) / rx_codeRate / rx_N_CBPS ) ; 


demod_seq= [] ; sof t_seq= [] ; stepsize=l; 
kpn = 1; 

samplesPerSymbol = 80; 

npt = 320 + 16 + 1; 7, skip preamble k SIGNAL symbol 

for k=l :rx_nSym 

D_hat=zeros(l,53) ; S_hat=zeros(l,48) ; 
kpn = kpn + 1; 

npt = npt + samplesPerSymbol; 
nn = npt + magicNums. u64; 
z = rx_wave (nn) ; 
y = fft(z) ./ H_k( : ) . • ; 

D_hat (magicNums . vl) = y (magicNums. v2) ; 

D_hat (magicNums . v3) = y (magicNums . v4) ; 

S_hat = D_hat (magicNums . i_data) ; 

W = sum( y (magicNums . i_symbol) .* magicNums. v_symbol ) * ( magicNums. pnVectorOq 
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44 Qbits=4 ; % number of bits for the quantizer 

45 [ Brx soft ] = DeMod_lla( rx_N_BPSC, Qbits, S_hat/rx_AMref ); 

46 demod_seq= [demod_seq Brx] ; y.size(demod_seq) 

47 sof t_seq= [sof t_seq soft] ; 

48 end 

49 

50 Rcv.bits = Demod_post_lla( rx_codeRate, rx_N_BPSC, rx_N_CBPS, demod_seq, soft_seq ) 

51 

52 serviceBits = Rcv_bits(l:16) ; 

53 dataBits = Rcv_bits(17: (16+8*rx_nBytes)); 

54 


processSIGNAL.m 

1 function [ rx.rateMode, rx_nSyms, rx_nBytes, signalErrFlags , signalErrMsg ] = processSIGNAL 

2 ’/PROCESSSIGNAL Process the 802.11a SIGNAL symbol 

3 V, SIGNALsym should probably be 64 simulation samples, beginning at 

4 7 , sample number 320+16+1. 

5 % 320 samples = skip over the PLCP Preamble (16us * 20MHz) ; 

6 ’/. 16 samples skips over the cyclic prefix of the SIGNAL symbol (0.8us * 20MHz) 

7 % 1 = MATLAB index is one-based ft we want first sample AFTER the cyclic prefix 

8 '/, MAGICNUMS = setMagicNumsO ; 

9 7, H_k = 64-sample FFT of multipath channel impulse response 

10 7, TX_RATEM0DE, TX_NBYTES are compared against contents of SIGNAL sym, to find errs 

11 % TX_RATEM0DE should be 0=48Mbps, 1=54, 2=12, 3=18, 4=24, 5=36, 6=6, 7=9Mbps 

12 7. These are just bits R1-R4 of the SIGNAL field 

13 7 , magicNums . xMODE can sort these into ascending bitrate 

14 % SIGNALERRFLAGS is an integer carrying a 3-bit field: 

15 7 , 0x01 => parity error in SIGNAL symbol 

16 % 0x02 => rx_rateMode "= tx_rateMode 

17 7 . 0x04 => rx_nBytes "= tx_nBytes 

18 % SIGNALERRMSG is a text string conveying the errors flagged by SIGNALERRFLAGS 

19 ’/. 

20 */. Subfunction for SIM_ER_11A.M 

21 ’/. 

22 % 2004, Gaylon R Lovelace, New Mexico State University 

23 ’/. 

24 

25 y = fft (SIGNALsym) ./ H_k(:).’; 

26 D_hat = zeros (1,64); 

27 D_hat (magicNums. vl) = y(magicNums.v2) ; 

28 D_hat (magicNums .v3) = y (magicNums. v4) ; 

29 S_hat = D_hat (magicNums . i_data) ; 


'/, time/freq mapping per 802.11a Fig 109 
*/, dunno why it’s done in two steps 
'/. skip elements from the NULL subcarrs 
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30 S_signal = round (real (S_hat) ) ; 

31 [ rx_rateMode rx.nSyms rx_nBytes rx_prtchkpass ] = parse_signal_lla( S signal ); 

32 

33 tx_codeRate = magicNums.RATE_dependent(magicNums.xMODE(tx_rateMode+l) ,3); 

34 tx_N_BPSC = magicNums.RATE_dependent(magicNums.xMODE(tx_rateMode+l) ,4) ; 

35 tx_N_CBPS = magicNums.RATE_dependent(magicNums.xMODE(tx_rateMode+l) ,5) ; 

36 tx_nSyms = ceil( ( 16 + 8*tx_nBytes + 6 ) / tx_codeRate / tx_N_CBPS ); 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 errParity = ( rx_prtchkpass ~= 1 ) ; 

48 errRateMode = ( rx.rateMode “= tx_rateMode ); 

49 errNBytes = ( rx_nBytes tx_nBytes ) ; 

50 signalErrFlags = errParity + 2*errRateMode + 4*errNBytes; 

51 

52 signalErrMsg = ’ ’ ; 

53 s i gn al Value sErr = ( ( rx_rateHode "= tx_rateMode ) I I ( rx_nSyms "= tx_nSyms ) I I ( 

54 if ( errParity ) 

55 signalErrMsg = ’Parity. 

56 end 

57 if ( errRateMode ) 

58 signalErrMsg = 

59 end 

60 if ( errNBytes ) 

61 signalErrMsg = 

62 end 

63 if ( signalErrFlags ) 

64 signalErrMsg = [ ’SIGNAL symbol error: ’ signalErrMsg ] ; 

65 end 

66 


[ signalErrMsg ’ rateMode: tx=’ num2str(tx_rateMode) ’ rx= 


[ signalErrMsg ’ nBytes: tx=’ num2str(tx_nBytes) ’ rx= ’ nui 


V, Just confirm that mWLAN : parse_signal_lla() uses the same nBytes — >nSymbols 
% formula that I just wrote above. If nBytes ft rateMode are both right, then 
’/, nSymbols should also agree, 
if ( ( rx_rateMode = tx_rateMode ) ftft ... 

( rx_nSyms *= ceil( ( 16 + 8*rx_nBytes + 6 ) / tx.codeRate / tx_N_CBPS ) ) 
error ( [ ’nSymbols calc mismatch: rx_nSyms=’ num2str(rx_nSyms) . . . 

’; rx_nBytes= ’ num2str(rx_nBytes) ’; rx_rateMode= ’ num2si 

end 


rmsDelay Calc . m 

1 fimction [ rmsDelay, meanExcess ] = rmsDelayCailc ( profile, Ts, isVoltage ) 

2 '/.RMSDELAYCALC Calculates multipath delay stats for Delay Profile 

3 [ rmsDelay, meanExcess ] = rmsDelayCalc ( profile, Ts, isVoltage ) 
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4 X takes a delay profile and calculates the mean excess delay and 

5 % the RMS delay spread. 

6 X If ISVOLTAGE is FALSE or omitted, then the delay profile is 

7 % taken to be power vs time. If ISVOLTAGE is TRUE, then the 

8 X profile is taken as voltage (or field strength, etc) vs time. 

9 X In either case, the profile is assumed to be on a linear scale. 


10 

11 

X rather than dB, and uniformly sampled at rate 1/Ts, 

V 

12 

X Subfunction for SIM_ER_11A.M 


13 

X 


14 

X 2004, Gaylon R Lovelace, New Mexico 

State University 

15 

X 


16 



17 

firstNonzero = min( find( profile ) ) 

• 

18 

taus =[0:1: length (profile )-l ] ’ 

.* Ts; 

19 

tau_0 = taus (firstNonzero) ; 


20 

taus = taus - tau_0; 


21 



22 

if ( (nargin<3) 1 1 (isVoltage~=0) ) 


23 

pdp « profile. "2; 


24 

else 


25 

pdp = profile; 


26 

end 


27 



28 

meanExcess = sum( pdp(:) .* taus ) 

/ sum( pdp( : ) ); 

29 

tauSqr„bar = sum( pdp(:) . * taus. "2 ) 

/ sum( pdp(: ) ); 

30 

rmsDelay = sqrt( tauSqr_bar - meanExcess "2 ); 

31 




dumpResultTable.m 

1 function resultTable = dumpResult Table ( resultMatrix, runParams ) 

2 XDUHPRESULTTABLE Produce a vide NxX string matrix with results from a set of sim^n runs 

3 X A simulation set consists of N simulation runs. 

4 X Each simulation run consists of many data packets. 

5 X Each run may in a set may do a different number of pkts. 

6 X Each run ends up producing ten scalar results, 

7 X so the N runs, together, produce the N-x-10 results matrix, 

8 X which is described column-by-column as: 

9 X resultMatrix - [ . . . 

10 X errBitCums, ... X 01 X Total number of bad data bits 

11 X berrPktCums, ... X 02 X Number of pkts with >= 1 bad data bit 

12 X perrPktCums, — X 03 X Number of pkts with bad SIGNAL parity 
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13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 


7. 

7 . 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 


serrPktCums , 
aerrPktCums , 
bitCums , 
pktCums, 
rmsDelays , 
rxPowers , 
estEbNoDBs , 


7. 04 7. Number of pkts with bad SIGNAL contents 
7 05 7 Number of pkts with any of above 3 errs 
7. 06 7. Total number of data bits transferred 
7 07 7 Total number of pkts transferred 
7 08 7 RMS Delay Spread for the channel used 
7 09 7 Avg rx/tx power ratio for the channel used 
7 10 7 Avg rx’d Eb/No for the channel & rateMode used 


All the ’instructions’ for a simulation set are contained in 

a big N-x-9 cell array: (must be cell array so last two columns can be strings) 

runParams = { ... 

7 01 7 Allows e.g. 24 runs to plot as 3 BER curves each 8 points 
... 7 02 7 Used to spec which 802.11a rateMode 

7 03 7 Independent variable for plots. Set zero if plotting vs dist 
7 04 7. Independent variable for plots. Set zero if plotting vs EbNo 
7 05 7 User data bytes per 802.11a packet 

7 06 7 Scales the ICS Telecom PDP (assumed to be for 1.0 W tx pwr) . 
7 07 7 Stop after M pkts, even if found few or no errors 
7 08 7 Stop early, if find this many pkts with >= 1 error. 

7 09 7 ’Subscriber’ string ID from ICS Telecom PDP * . csv file 
7 10 7 Filename of ICS Telecom PDP * ■ csv file 


CurveNum, 

Mbps, 

EbNoDB , 

Distance, 

pktLen , 

txPwr, 

MaxPkts, 

MinErrPkts , 

Channel, 

pdpFile 




beCums pbeCums ppeCums pseCums paeCums 


Subfunction for SIM_ER_11A.M 


6 

nBits 


7 

nPkts 


8 9 

BER rmsDelay 


10 

rxPower estl 


2004, Gaylon R Lovelace, New Mexico State University 


TresultMatrix = [ errBitCums, berrPktCums, perrPktCums , serrPktCums, aerrPktCums, bitCums, ] 
nRuns = size( resultMatrix, 1 ); 

if ( nRuns "= size( runParams, 1 ) ); error ( ’RESULTMATRIX, RUNPARAMS be same height.’ ); 
runParamNums = cell2mat ( runParams (.*,[1:8])); 

resultMatrix(: ,8) = resultMatrix (: ,8) * le9; 7 Display rmsDelaySpread in nanoSecs 
7 don’t 7 resultMatrix (: ,9) = resultMatrix(: ,9) * le9; 7 Display rxPower in nanoWatts 
nonzeroBits = find( resultMatrix (: ,6) ); 

bers = resultMatrix( : ,6) ; 7 BERs default to zero (for runs not completed, etc) 

bers (nonzeroBits) = resultMatrix (nonzeroBits , 1) ./ resultMatrix (nonzeroBits, 6) ; 

resultTable = [ ’7 ’ datestr(now) ]; 

resultTable = strvcat( resultTable, ’7 1 2 3 4 5 6 

resultTable = strvcat( resultTable, ’7 beCums pbeCums ppeCums pseCums paeCums nBits nPl 
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58 f onnatStr = > ‘/.8dy.8d , /.8d , /.8d , /.8d*/.9dy.8d , /.10 . 3g*/,9 . If ‘/.10 . 3 g y,8 . If %10dX6dX8 . 2f‘/.9 . 4g*/.9 . 4 g y.9dy.9d*/.9dy.l3! 

59 for k = [ 1 : nRuns ] 

60 tmpStr = [ charfrunParamsCk.g)) } ]; 

61 thisLine = sprintf ( formatStr, resultMatrixfk, [1:7]) , bers(k), resultMatrix(k, [8: 10.' 

62 resultTable = strvcatC resultTable, thisLine ); 

63 end 

64 


RunParams_BERDistAvg.m 

This last file is one sample parameter file. Similar files are used to define all the different 
simulation sets we did. This one in particular sets up a long set of simulations to find 
BER-vs-(Tx-Rx Distance) averaged over several spots at the Gusevl, Sitel location. 


1 runParams = { 


2 

y. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

3 

% CurveNum 

Mbps 

EbNoDB 

Distance 

txPwr 

pktLen 

MaxPkts 

MinErrPkts Channel 

pdpF: 

4 


1 

12 

0 

1000 

0.001 

128 

500 

10000 

'^lOGOm!* 

’DistA 1 

5 


1 

12 

0 

500 

0.001 

128 

500 

10000 

>rx500ml> 

^DistA" 

6 


1 

12 

0 

200 

0.001 

128 

500 

10000 

’rx200ul’ 

'DistA" 

7 


1 

12 

0 

100 

0.001 

128 

500 

10000 

’rxlOOml 7 

'DistA^ 

8 


1 

12 

0 

50 

0.001 

128 

500 

10000 

’rxBOml * 

’DistA^ 

9 

10 

•/. 

1 

12 

0 

20 

0.001 

128 

500 

10000 

, rx20ml > 

’DistA ? 

11 


2 

12 

0 

1000 

0.001 

128 

500 

10000 

, rxl000m2 > 

’DistA 1 

12 


2 

12 

0 

500 

0.001 

128 

500 

10000 

’rx50Qm2* 

’DistA 1 

13 


2 

12 

0 

200 

0.001 

128 

500 

10000 

* rx200m2 3 

3 List A' 

14 


2 

12 

0 

100 

0.001 

128 

500 

10000 

3 rxlOQw2 3 

’DistA^ 

15 


2 

12 

0 

50 

0.001 

128 

500 

10000 

3 rx50m2 3 

’DistA" 

16 

17 

y. 

2 

12 

0 

20 

0.001 

128 

500 

10000 

3 rx20m2 J 

’DistA 1 

18 


3 

12 

0 

1000 

0.001 

128 

500 

10000 

'rxlOOOmS' 

’DistA 1 

19 


3 

12 

0 

500 

0.001 

128 

500 

10000 

’rxBOOmS’ 

’DistA 1 

20 


3 

12 

0 

200 

0.001 

128 

500 

10000 

, rx200m3* 

’DistA^ 

21 


3 

12 

0 

100 

0.001 

128 

500 

10000 

'rxlOOmS* 

’DistA 1 

22 


3 

12 

0 

50 

0.001 

128 

500 

10000 

3 rx50m3 9 

’DistA^ 

23 

24 

% 

3 

12 

0 

20 

0.001 

128 

500 

10000 

3 rx20m3 * 

*DistA\ 

25 


4 

12 

0 

1000 

0.001 

128 

500 

10000 

'rxlOOOn^ 

’DistA 1 

26 


4 

12 

0 

500 

0.001 

128 

500 

10000 

’rx50Qm4’ 

^DistA 1 

27 


4 

12 

0 

200 

0.001 

128 

500 

10000 

’rx20Qn4’ 

'DistA 1 

28 


4 

12 

0 

100 

0.001 

128 

500 

10000 

’rxlOQn^ 

’DistA^ 

29 


4 

12 

0 

50 

0.001 

128 

500 

10000 

, rx50m4 7 

’DistA 1 
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30 4 12 0 20 0.001 128 500 10000 , rx20m4’ 'DistA' 

31 >; 

32 skipDemod = 0; demodulate data to find BER 

33 vsEbNo = 0; '/, SNR determined from ICS Telecom 

34 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 


IEEE 802 . 11 b Simulation Code 

MultipathEbNoSim.m 

/(Program to simulate BER as a function of Eb/No using Martian Multipath 
‘/.Data 

7. IEEE 802. lla/b/g Toolbox, v.3.02, March 2004 

7 , All rights reserved, (c) CommAccess Technologies, 2001-2004. 

7 % Modified by Dr. Deva K. Borah, NMSU 
7 . Modified by Anirudh Daga, NMSU 

clear; format compact; close all ; clc ; 

opt_mode=l;7, opt_mode = 1 for 1 Mbps, 2 for 2 Mbps, 3 for 5.5 Mbps, 4 for 11 Mbps 
rake = 0; 

barker= [1 -111-1111-1-1 -1]; 

CCK_PBCCtx=0;7. 0 for CCK, 1 for PBCC 

cckdemodtype=l; 70 for suboptimal CCK demod, 1 for near-optimal 
coherent=0 ; 7, 1 for coherent demod, 0 for otherwise 

p=18;7, p = the number of sample per channel symbol 
chiprate=ll; % Mega cps 
deltaT*=l/ chiprate/p ; 
nsample=p; 

‘/.Pulse shaping sequence 

Nspan=5 ; 

ns=p; 

alpha=0 . 35 ; 

hTx=SRRC (alpha, ns , Nspan) ; 

hRx=hTx/nsample; % This scaling operation is simply for the purpose of plotting. 


xxxmmxxmmxmxmmmmmmmmmmxmmmmmxmmxxx 


v. 


ICS TELECOM POWER DELAY PROFILE INCORPORATION 


X 


;xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx%xxxxxmxx 


tl = ’C:\DocumentsPapers\NASA Documents\PDPSimulationData0ct2004\Gusev3_Sitel_100m.txt’; 
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34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 


7. t2 = ’C:\DocumentsPapers\NASA Documents\PDPSimulationData0ct2004\Gusev3_Site2_100m.txt’ 
7 , t3 = 1 C : \DocmnentsPapers\NASA Documents\PDPSimulationData0ct2004\Gusevl_Sitel_100m.txt’ 
files = strvcat(tl); 
for kk = 1: size (files, 1) 

PowerDelayProfile = dlmread (files (kk, :)) ; 

ImpulseResponse = PDP2IMP (PowerDelayProfile) * ; 

ImpulseResponse = Missing_Time_Insert ion (ImpulseResponse) ; 



ImpulseResponse = (ImpulseResponse (: ,l)/norm(ImpulseResponse(: ,1))) ’ ; 
hfl = resample(ImpulseResponse,99,50,l) ; 
hfl = hf 1/norm (hfl) ; 

IRLength = length (ImpulseResponse) ; 


if (opt_mode>4) , opt_mode=4; end 
if (opt_mode<l) , opt_mode=l; end ‘/.protection 

7 , start waveform generation 

L_S=1; 7 . 1 for the long preamble, 0 for the short preamble 

if (opt_mode==l) Rate_modetx=10 ; end; 
if (opt_mode=2) Rate_modetx=20 ; end; 
if (opt_mode==3) Rate_modetx=55 ; end; 
if ( opt _mode=4 ) Rate_modetx=110; end; 


if (opt_mode==3) 

[codevec55,bitvalvec55,phchodd55,phchevn55,phvec55]=cck55demodsupport ; 

end; 


if (opt_mode==4) 

[codevec , bitvalvec , phchodd , phchevn , phvec] =cckl ldemodsupport ; 

end; 


if (CCK_PBCCtx == 1) 
ph0=pi ; 

end 7 , for (PBCC) BER simulation purpose 

nByte=100; 
snrdb = 1000; 
for xx = 1: length (snrdb) 
tic 

bitcounter=0; 
wrongcounter=0 ; 
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79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 
123 


for yy = 1:1000 
tic 

files (kk, : ) 
nBittx=nByte*8 ; 

PSDU=BitGen(nBittx) ; 

PSDULength * length (PSDU) ; 
b23=[0 CCK_PBCCtx] ; 

/.Generate the 3 c ramble! spread data 

[IQsequence , xBit , phO] = PLCP_llb(opt_mode,L_S,PSDU,b23) ; 
if (rake =1) 

TxdSpreadPreamble = IQsequence(l: 128*11) ; 

end 

nIQ=length(IQsequence) ; 
if (b23(2)— 1) 

tap_llb=[l 011011111110 1]; 

VitrB_ini(tap_llb) 

end 

'/.Generate Rayleigh fading 

rayleigh_f ading = (l/sqrt(2))*(randn(size(hf 1)) + j*randn(size(hf 1))) ; 
hf = hfl.*rayleigh_f ading; 

CIR = conv(conv(hf ,hTx) ,hRx) ; 

CIR = CIR(10*ns+l :ns : length(CIR) -10*ns) ; 

'/. generate Tx waveform 

S=zeros (1 ,ns*nIQ) ; 

S ( 1 : ns : ns*nIQ)=IQsequence ; 

S=conv(S,hTx) ; 

Sl=conv(S,hf ) ; 

'/» Add noise to signal samples 
nx=length(Sl) ; 

if (opt_mode==l) noisestd=sqrt (p*ll/ (2*10“ (snrdb(xx) /10) ) ) ; end; 
if (opt_mode==2) noisestd=sqrt(p*ll/(4*10“(snrdb(xx)/10))) ; end; 
if (opt_mode==3) noisestd=sqrt(p/(10“(snrdb(xx)/10))) ; end; 
if (opt_mode==4) noisestd=sqrt(p/(2*10~ (snrdb(xx)/10))) ; end; 

wi=randn ( 1 , nx) *noisestd ; 
wq=randn(l ,nx)*noisestd; 
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124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 


if (opt_mode — 1) SI = SI + wi; else SI = Sl+(wi+j*wq) ; end; 

’/, Receiver section starts 

Sr=conv(Sl ,hRx) ; 

ExPar.llbl 

rcvl = Sr(10*ns+l:ns:length(Sr)-10*ns) ; 
nr=length(rcvi) ; 

if (rake == 1) 

RxdSpreadPreamble = rcvl (1: length (TxdSpreadPreamble)+length (rcvl) -nlQ) ; 
hhat = ChannelEst _NoFor (RxdSpreadPreamble, TxdSpreadPreamble) ; 

[maxval , maxind] =max (hhat ) ; 
k_header_end = (144+48) *11+1; 
elseif (rake — 0) 

[maxval, maxind] - max(CIR) ; 
k_header_end=( 144+48) *11+1+ (maxind-1) ; 

end 


rcv3=rcvl (k_header_end: length (rcvl) ) ; 
if ( (Rate_modetx==55)&(cckdemodtype==l) ) 
if (rake == 1) 

d_rcv = CCK55DemodRake_NoFor(rcv3,nBittx,codevec55,bitvalvec55, . . . 
phchodd55 , phchevn55 , phvec55 , phO , PSDULength , hhat ) ; 
elseif (rake == 0) 

d_rcv = cck55demod(rcv3,nBittx,codevec55,bitvalvec55,phchodd55, . . . 
phchevn55,phvec55,ph0) ; 

end 

elseif ( (Rat e_modetx“ 110)4 ( cckdemodtype== 1 ) ) 
if (rake = 1) 

d_rcv=CCKHDemodRake_NoFor(rcv3,nBittx,codevec,bitvalvec, . . . 
phchodd , phchevn , phvec , phO , PSDULength , hhat ) ; 
elseif (rake == 0) 

d_rcv=ccklldemod(rcv3,nBittx,codevec,bitvalvec, phchodd, . . . 
phchevn, phvec, phO) ; 

end 

else 


if (rake — 1) 
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169 d_r cv=DenKxLllbv4 (Rat e^modetx , nBittx, phO , CCK_PBCCtx , coherent , 

1 70 r cv3 , hhat , PSDULength) ; 

171 elseif (rake == 0) 

172 d_rcv=Demod_l lbv2 (Rate_modetx , nBittx , phO , CCK_PBCCtx , coherent , rcv3) ; 

173 end 

174 

175 end; 

176 

177 dd_rcv=Descrambler_llb(d_rcv(l: 1 : nBittx)) ; 

178 Err=(dd_rcv(l:nBittx)~=PSDU) ; 

179 ErrorKnt=sum(Err) ; 

180 bit counter * bit counter + nBittx; 

181 wrongcounter - wrongcounter + ErrorKnt 

182 ber = wrongcounter /bit counter 

183 yy 

184 snrdb(xx) 

185 if (snrdb(xx)<10 kk wrongcounter >= 20000) 

186 break; 

187 end 

188 if (snrdb(xx)>10 kk wrongcounter >= 10000) 

189 break; 

190 end 

191 toe 

192 end 

193 finber(xx) - ber; 

194 toe 

195 end 

196 BitErrRate(kk f : ) = finber; 

197 end 


mainprog_l lb jntdist .m 

1 % Program to simulate BER as a function of distance using Martian multipath data 

2 

3 7* IEEE 802.11a/b/g Toolbox, v.3.02, March 2004 

4 % All rights reserved, (c) CommAccess Technologies, 2001-2004. 

5 

6 7. Modified by Dr. Deva K. Borah, NMSU 

7 7, Modified by Anirudh Daga, NMSU 

8 

9 clear all; format compact; close alljclc; 

10 barker=[l -111-1111-1-1 -1]; 

11 rake = 1; 
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12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 
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32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 


opt_mode=4 ; 7. opt_mode = 1 for 1 Mbps, 2 for 2 Mbps, 3 for 5.5 Mbps, 4 for 11 Mbps 
CCK_PBCCtx=0 ; '/, 0 for CCK, 1 for PBCC 

cckdemodtype=l ; %0 for suboptimal CCK demod, 1 for near-optimal 
coherent=0 ; '/, 1 for coherent demod, 0 for otherwise 


p=18;% p = the number of sample per channel symbol 
chiprate=ll; % Mega cps 
deltaT=l/ chiprate/p; 
nsample=p ; 

’/, Pulse shaping sequence 


Nspan=5 ; 


ns=p; 


alpha=0 . 35 ; 

hTx=SRRC ( alpha , ns , Nspan) ; 
hRx=hTx/nsample ; 




7. ICS TELECOM POWER DELAY PROFILE INCORPORATION 7. 


xxmmxxxxmmxmmmxmmxmxxxmmmmmmmmmxmmxx 


PowerDelayProfile = dlmreadC’ C:\DocumentsPapers\NASA Documents\ . . . 

PDPSimulationData0ct2004\Gusevl_Sitel_100m.txt ’ ) ; 
[ImpulseResponse , TruePower] = PDP2IMPDistance (PowerDelayProfile) ; 
ImpulseResponse = Missing_Time_Insertion (ImpulseResponse) ; 

ImpulseResponse = ImpulseResponse ( : ,1) ’ ; 


mmmmmmmmmmmmmmmr/mmmmmmmmmn 

% OVERSAMPLING THE MULTIPATH CHANNEL IMPULSE RESPONSE % 

mmrammrammmffimramrarammmramrammnmmm 


hfl = resample (ImpulseResponse, 99, 50, 1) ; 
hfl = (norm ( ImpulseResponse ) /norm (hfl)) . *hf 1 ; 
IRLength = length(ImpulseResponse) ; 


if (opt_mode>4) , opt_mode=4; end 
if (opt_mode<l) , opt_mode=l ; end % protection 


if 

if 

if 

if 


(opt_mode==l) 

(opt_mode“2) 

(opt_mode==3) 

(opt_mode==4) 


Rate_modetx=10; end; 
Rate_modetx-20; end; 
Rate_modetx=55 ; end; 
Rate_modetx=l 10 ; end ; 


if (opt_mode==3) 

[codevec55 , bitvalvec55 , phchodd55 , phchevn55 ,phvec55] =cck55demodsupport ; 

end; 


if (opt_mode==4) 
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57 [codevec , bitvalvec , phchodd , phchevn , phvec] =cckl ldemodsupport ; 

58 end; 

59 

60 if (CCK_PBCCtx == 1) , 

61 phO=pi ; 

62 end % for (PBCC) BER simulation purpose 

63 

64 nByte=500 ; 

65 bitcounter=0 ; 

66 wrongcounter=0 ; 

67 

68 for yy = 1:5000 

69 randn(’ state’ , sum (100* clock)) 

70 tic 

71 % start: waveform generation 

72 L_S=1; X 1 for the long preamble, 0 for the short preamble 

73 nBittx=nByte*8 ; 

74 PSDU=BitGen(nBittx) ; 

75 PSDULength = length (PSDU) ; 

76 b23=[0 CCK_PBCCtx] ; 

77 

78 [IQsequence, xBit,phO]=PLCP_llb(opt_mode,L_S,PSDU,b23) ; 

79 TxdSpreadPreamble = IQsequence (1 : 128*11) ; 

80 nIQ=length (IQsequence) ; 

81 

82 if (b23(2)==l) 

83 tap_llb=[l 0110111111101]; 

84 VitrB_ini (tap_l lb) 

85 end 

86 

87 '/.Generate Rayleigh fading 

88 

89 rayleigh_fading = (l/sqrt(2))*(randn(size(hfl)) + j*randn(size(hf 1))) ; 

90 hf = hf 1. *rayleigh_fading; 

91 CIR = conv(conv(hf ,hTx) ,hRx) ; 

92 CIR * CIR(10*ns+l:ns:length(CIR)-10*ns); 

93 

94 '/, generate Tx waveform 

95 

96 S=zeros(l,ns*nIQ) ; 

97 S(l:ns:ns*nIQ)=IQsequence; 

98 S=conv(S,hTx) ; 

99 S=sqrt(le-3)*S; 

100 Sl=conv(S,hf) ; 

101 nx=length(Sl) ; 
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% Add noise to signal samples 
sigma = 4.8663e-7; 

wi=randn(l ,nx)*sigma; wq=randn(l ,nx)* sigma; 

if (opt_mode == 1) S2 = SI + wi; else S2 = SI + (wi+j*wq) ; 

'/.Receiver section starts 
Sr=conv(S2,hRx) ; 

rcvl = Sr(10*ns+l:ns:length(Sr)-10*ns); 
if (rake == 1) 

RxdSpreadPre amble = rcvl(l:length(TxdSpreadPreamble)+length(rcvl)-nIQ-l) ; 
hhat = ChannelEst_NoFor (RxdSpreadPreamble .TxdSpreadPreamble) ; 

[maxval , maxind] =m ax (hhat) ; 
elseif (rake == 0) 

[maxval, maxind] = max(abs(CIR)) ; 

end 

if (rake == 0) 

k_header_end= ( 144+48) * 11+1+ (maxind- 1) ; 
elseif (rake == 1) 

k_header_end= ( 144+48) +11+1; 

end 

nr=length(rcvl) ; 

rcv3=rcvl (k_header_end : length (rcvl ) ) ; 

if ( (Rate_modetx==55) & (cckdemodtype==l) ) 
if (rake == 0) 

[d_rcv] = cck55demod(rcv3,nBittx,codevec55,bitvalvec55, . . . 
phchodd55 , phchevn55 ,phvec55 , phO) ; 
elseif (rake =1) 

[d_rcv] = CCK55DemodRake_NoFor(rcv3,nBittx, codevec55 ,bitvalvec55 , . . . 
phchodd55 , phchevn55 ,phvec55 , phO , PSDULength , hhat ) ; 

end 


elseif ( ( Rat e_modet x== 1 1 0 ) & ( c ckdemodt ype== 1 ) ) 
if (rake == 0) 

[d_rcv]=ccklldemod(rcv3,nBittx,codevec,bitvalvec,phchodd, . . . 
phchevn,phvec,phO) ; 
elseif (rake =1) 

[d_rcv]=CCKHDemodRake_NoFor(rcv3,nBittx,codevec .bitvalvec, . . . 
phchodd , phche vn , phve c , phO , PSDULength , hhat ) ; 

end 

else 

if (reike == 0) 

d_rcv=Demod_l lbv2 (Rate_modetx .nBittx , phO , CCK_PBCCtx , coherent , rcv3) ; 
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147 elseif (rake == 1) 

148 d_rcv=Demod_llbv4(Rate_modetx,nBittx,phO,CCK_PBCCtx, coherent, . . . 

149 rcv3 , hhat , PSDOLength) ; 

150 end 

151 end; 

152 

153 dd_rcv=Descrambler_l lb (d_rcv (1:1: nBittx) ) ; 

154 Err= (dd_rcv(l: nBittx) ~=PSDU) ; 

155 ErrorKnt=snm(Err) ; 

156 bitcounter = bitcounter + nBittx; 

157 wrongcounter = wrongcounter + ErrorKnt 

158 ber = wrongcounter/bitcounter 

159 yy 

160 

161 if (wrongcounter >= 10000) 

162 break; 

163 end 

164 

165 toe 

166 end 

167 


PERDistanceRake_NoFor.m 

1 clear all; format compact; close all;clc; 

2 

3 

4 '/. IEEE 802. lla/b/g Toolbox, v.3.02, March 2004 

5 V, All rights reserved, (c) CommAccess Technologies, 2001-2004. 

6 

7 */. Modified by Anirudh Daga, NMSU 

8 

9 7. 

10 opt_mode=l;% opt .mode = 1 for 1 Mbps, 2 for 2 Mbps, 3 for 5.5 Mbps, 4 for 11 Mbps 

11 rake = 0; 

12 CCK_PBCCtx=0;‘/, 0 for CCK, 1 for PBCC 

13 cckdemodtype=l; 7,0 for suboptimal CCK demod, 1 for near-optimal 

14 coherent=0 ; 7, 1 for coherent demod, 0 for otherwise 

15 

16 p=18; 7, p = the number of sample per channel symbol 

17 chiprate=ll; % Mega cps 

18 deltaT=l/chiprate/p; 

19 nsample=p; 
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*/. Pulse shaping sequence 


Nspan=5 ; 


ns=p; 


alpha=0 . 35 ; 

hTx=SRRC ( alpha , ns , Nspan) ; 
hRx=hTx/ ns ample ; 


xxmxxxxmxxxmxxxxxmmxxxxxxxmxxxmxmxxxxmxmxmxxxxxxxmxm 


7 . ICS TELECOM POWER DELAY PROFILE INCORPORATION '/. 

xxxmmxmxmxmxmxmmxxmxmmxmmmxmmmmmxxmm 


PoverDelayProf ile = dlmread(’ C:\DocumentsPapers\NASA Documents\ . . . 

PDPSimulationData0ct2004\Gusevl_Sitel_1000m.txt ’ ) ; 
[ImpulseResponse , TruePover] = PDP2IMPDistance (PowerDelayProf ile) ; 
TruePower 

ImpulseResponse = Missing_Time_Insertion(ImpulseResponse) ; 
ImpulseResponse - ImpulseResponse (: , 1) ’ ; 


xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 


% 


OVERSAMPLING THE MULTIPATH CHANNEL IMPULSE RESPONSE 


7 . 


't*nl%hHl%nl*nf*hhhhhhHhl»hhhnhhhhhhhhhhhhhhhhhhhhUhU/*hUhl*i»hf»Ul%HhUl9/\ 

hf 1 = resample (ImpulseResponse ,99,50,1); 
hfl = (norm(ImpulseResponse)/norm(hf 1)) .*hfl; 

IRLength = length(ImpulseResponse) ; 


if (opt_mode>4) , opt_mode=4; end 
if (opt_mode<l) , opt_mode=l; end '/(protection 


if 

if 

if 

if 


(opt_mode==l) 

(opt_mode=2) 

(opt_mode==3) 

(opt_mode=4) 


Rate_modetx=10 ; end ; 
Rat e_modetx=20 ; end ; 
Rate_modetx=55 ; end; 
Rate_modetx=110 ; end; 


if (opt_mode=3) 

[codevec55 , bitvalvec55 , phchodd55 , phchevn55 , phvec55] =cck55demodsupport ; 

end; 


if (opt_mode=4) 

[codevec , bitvalvec , phchodd , phchevn , phvec] =cckl ldemodsupport ; 

end; 

if (CCK.PBCCtx == 1), 
ph0=pi ; 

end 7, for (PBCC) BER simulation purpose 
nByte=50 ; 
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65 packet counter=0; 

66 vrongcounter=0 ; 

67 

68 for yy = 1 : 1000 

69 '/, start waveform generation 

70 L_S=1; '/, 1 for the long preamble, 0 for the short preamble 

71 nBittx=nByte*8 ; 

72 PSDU=BitGen(nBittx) ; 

73 PSDULength = length (PSDU) ; 

74 b23= [0 CCK.PBCCtx]; 

75 

76 [IQsequence , xBit , phO] =PLCP_1 lb (opt_mode , L_S , PSDU , b23) ; 

77 TxdSpreadPreamble = IQsequence (1 : 128*11) ; 

78 nlQ-length(IQsequence) ; 

79 

80 if (b23(2)— 1) 

81 tap_lib=[l 0110111111101]; 

82 VitrB_ini (tap_ 1 lb) 

83 end 

84 

85 ‘/.Generate Rayleigh fading 

86 

87 rayleigh_fading = (l/sqrt(2))*(randn(size(hf 1)) + j*randn(size(hf 1))) ; 

88 hf = hf l.*rayleigh_fading; 

89 CIR = conv(conv(hf ,hTx) ,hRx) ; 

90 CIR = CIR(10*ns+l:ns:length(CIR)-10*ns); 

91 % 

92 '/. generate Tx waveform 

93 

94 S=zeros(l,ns*nIQ) ; 

95 S(l:ns:ns*nIQ)=IQsequence; 

96 S=conv(S,hTx) ; 

97 Sl=conv(S,hf ) ; 

98 nx=length(Sl) ; 

99 

100 

101 '/. Add noise to signzd samples 

102 sigma = 4.8663e-7; 

103 wi=randn(l,nx)*sigma; wq=randn(l,nx)* sigma; 

104 

105 if (opt_mode == 1) S2 = SI + wi; else S2=Sl+(wi+j*wq) ; 

106 

107 /^Receiver section starts 

108 Sr=conv(S2,hRx) ; 

109 
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111 
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113 

114 
| 115 

I 116 

i 117 

I 118 

119 
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147 
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150 
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152 
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154 


ExPar.llbl 

rcvl = Sr(10*ns+1 :ns:length(Sr)-10*ns) ; 
if (rake — 1) 

RxdSpreadPreamble = rcvl(l:length(TxdSpreadPreamble)+length(rcvl)-nIQ) ; 
hhat = ChannelEst _NoFor (RxdSpreadPreamble, TxdSpreadPreamble) ; 

[maxvai .maxind] = max (hhat); 
elseif (rake — 0) 

[maxvai , maxind] =max ( CIR) ; 

end 

if (rake == 0) 

rcvl * rcvl (1:2112); 
elseif (rake == 1) 

rcv2 = rcvl (1:2112 + length(hhat) -1); 

temp = repmat((l:length(hhat)) ’ ,1,2112) + repmat(0:2112-l,length(hhat) ,1) ; 
rcv3 * conj(repmat(hhat(:),l,2112)).*rcv2(temp); 
rcvl = sum(rcv3) ; 

end 

nr=length(rcvl) ; 

ExPar_llb2v2 

z=(T=l); 

b=0; 

signal* (z(b+l:b+8)==l) ; 

Rat e_mode=sum ( signal . * (2 (7: -1:0))); 

service* (z (b+9 : b+16) =1) ; 

CCK_PBCCrx*service (4) ; 
leng=(z(b+17:b+32)*=l) ; 

CRC* (z (b+33 : b+48) =1) ; 7. [PPDU( 177 : 192) ] 

CRC=ones (1,16) -CRC ; 

T= [signal service leng CRC]; 
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155 che cksum=CRCck_ 1 lb (T) 

156 

157 if (checksum ~=0) 

158 wrongcounter = wrongcounter + 1; 

159 end 

160 wrongcounter 

161 packet counter = packet counter + 1 

162 per = wrongcounter /packet counter 

163 if wrongcounter > 1000; break; end; 

164 yy 

165 end 

166 


PDP2IMRm 

1 ^Function to convert Martian PDP into channel impulse response 

2 %By Anirudh Daga, NMSU 

3 

4 function [ImpulseResponse] = PDP2IMP (PowerDelayProf ile) 

5 

6 PowerDelayProfileC: ,1) = PowerDelayProfileC: ,1) - 144.82; 

7 PowerDelayProfileC: ,1) = (le-3) .*(10. ~(. 1* (PowerDelayProf ile (: ,1)))) ; 

8 PowerDelayProf ile (: ,2) = PowerDelayProfileC: ,2) - PowerDelayProf ile (1 ,2) ; 

9 PowerDelayProf ile = PowerDelayProf ile 1 ; imp = zeros (size (PowerDelayProf ile)) ; 

10 imp(l , : ) = sqrt (PowerDelayProf ile (1 , : )) ; imp (2, : ) = PowerDelayProf ile (2, :) ; 

11 ImpulseResponse = imp; 


PDP2IMPDistance.m 

1 ^Function to obtain channel impulse response and true power from Martian PDPs at a given 

2 ^distance 

3 %By Anirudh Daga, NMSU 

4 

5 function [ImpulseResponse, TruePower] - PDP2IMPDistance (PowerDelayProf ile) 

6 

7 PowerDelayProf ile ( : , 1) = PowerDelayProf ile(: , 1) - 144.82; 

8 PowerDelayProfileC : ,1) = (le-3) .*(10. * ( . 1* (PowerDelayProf ile ( : , 1) ) ) ) ; 

9 TruePower = sum(PowerDelayProfile(: , 1)) ; 

10 PowerDelayProf ile = PowerDelayProf ile ’ ; 

11 imp » zeros (size (PowerDelayProf ile)) ; 

12 imp(l,:) = sqrt (PowerDelayProf ile (1, :)) ; 
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13 imp (2 , : ) = PowerDelayProf ile (2 , : ) ; 

14 ImpulseResponse = imp’; 


Missing_Time_Iiisertion.m 

1 XMissing time insertion to fill in zeros 

2 %By Anirudh Daga, NMSU 

3 

4 function [PowerDelayProfile] = Missing_Time_Insertion(pdp) 

5 

6 pdp(: ,2) = rounddec(pdp(: ,2) ,2) ; 

7 b( : ,2) = [pdp(l,2) :0.01:pdp(end,2)] ’ ; 

8 b(:,2) * rounddec(b(: ,2) ,2) ; 

9 [dum , idx] = ismember (pdp( : ,2) ,b( : ,2) ) ; 

10 b(idx,l) =pdp(:,l); 

11 PowerDelayProfile = b; 


rmsdelayspread.m 

1 7 . Program to calculate RMS Delay Spread from PDP 

2 clear; close all; dc; format long; 

3 

4 pdp = dlmread(’ C:\DocumentsPapers\NASA DocumentsV . . 

5 PDPSimulationData0ct2004\AntennaHeight\Gusevl. . . 

6 Site3\Gusevl_Site3_2m_llb.txt’); 

7 pdp(: ,1) = pdp(: ,1)-144.82; 

8 pdp(: , 1 ) = (le-3).*(10.“(.l*(pdp(:,l)))); 

9 pdp = pdp . ’ ; 

10 pdp(2, : ) = (le-6)*(pdp(2, :)) ; 

11 tatibar = sum(pdp(l, :) .*pdp(2, : ))/sum(pdp(l, :)) ; 

12 tausqrbar = sum(pdp(l, :) .*(((pdp(2, :)) .‘2)))/sum(pdp(l, :)) ; 

13 rmsdsp * sqrt (tausqrbar - (taubar“2) ) 


ChannelEst_NoFor.m 

1 */. Channel Estimation 

2 */. By Anirudh Daga, NMSU 

3 

4 

5 function [hhat] = ChannelEst_NoFor(RxdSpreadPreamble,TxdSpreadPreamble) 
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6 

7 n = length (TxdSpreadPreamble) ; 

8 phi = zeros(l.n); 

9 

10 for j= 0:n-l 

11 ind = 1 : n-j ; 

12 phi(j+l) = sum (TxdSpreadPreamble (ind) .* TxdSpreadPreamble (ind + j)); 

13 end 

14 

15 

16 autocorr = toeplitz(phi) ; 

17 

18 IRLengthChipRate = length(RxdSpreadPreamble)-length(TxdSpreadPreamble)+l; 

19 

20 temp = repmat ( (1 : IRLengthChipRate) ’ , 1 , length (TxdSpreadPreamble) ) . . . 

21 + repmat (0 : length (TxdSpreadPreamble) -1 , IRLengthChipRate , 1 ) ; 

22 path = RxdSpreadPreamble (temp) ; 

23 corr = (path*TxdSpr e adPre amble .’).’; 

24 

25 Rbb = autocorr ( 1 : IRLengthChipRate , 1 : IRLengthChipRate) ; 

26 hhat = (inv (Rbb)* corr’) ’ ; 

27 


cckl ldemodsupport.m 

1 */3y Dr. Deva K. Borah, NMSU 

2 /(Please do not distribute without the permission of the author 

3 function Ccodevec,bitvalvec,phchodd,phchevn,phvec]=ccklldemodsupport(void) ; 

4 totalbitvecs=2“8; 

5 veccntr=zeros(l,8) ; 

6 veccntr(l)=-l; 

7 for hypono=l :totalbitvecs 

8 veccntr(l)®veccntr(l)+l; 

9 if (veccntr(l)>=2) 

10 veccntr(l)=0; flowing®!.; 

11 for jj=2:8 

12 if (flowing=l) 

13 veccntr(jj)=veccntr(jj)+l; 

14 if ( veccntr ( j j ) >=2) veccntr(j j)=0; flowing=l; else flowing=0; end; 

15 end; ‘/if flowing==l 

16 end; )ifor jj 

17 end; 7, if veccntr(l)>=2 

18 bitvalvec(hypono, : )=veccntr; 
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19 

20 /(Generate QPSK symbols 

21 for ii~l:4 

22 bitval=2*bitval vec (hypono, 2* (ii-l)+l)+bitvalvec (hypono, 2*ii) ; 

23 switch bitval 

24 case 0, phvec (hypono ,ii)=0; 

25 case 1, phvec (hypono, ii)=pi/2; 

26 case 2, phvec (hypono, ii)=pi; 

27 case 3, phvec (hypono, ii)=3*pi/2; 

28 end; /(switch 

29 end; 7.for ii 

30 codevec (hypono ,1:8)= [exp ( j * (phvec (hypono , 1) +phvec (hypono , 2) +phvec (hypono , 3) + . . . 

31 phvec (hypono , 4) ) ) , exp ( j * (phvec (hypono , 1 ) +phvec (hypono , 3) + . . . 

32 phvec (hypono, 4))) , exp(j* (phvec (hypono, 1 )+phvec (hypono, 2) + . . . 

33 phvec (hypono , 4) ) ) , -l*exp(j*(phvec(hypono,l)+phvec(hypono,4))) , . . . 

34 exp ( j * (phvec (hypono , 1) +phvec (hypono , 2) +phvec (hypono , 3) ) ) , . . . 

35 exp ( j * (phvec (hypono , 1 ) +phvec (hypono , 3) ) ) , “l*exp(j*(phvec(hypono,l)+. . . 

36 phvec (hypono , 2) ) ) , exp ( j *phvec (hypono , 1 ) )] ; 

37 end; */,for bitno 

38 

39 ^Define even and odd phase changes 

40 phchevn(l)=0; 

41 phchevn(2)=pi/2; 

42 phchevn (3)=3*pi/2 ; 

43 phchevn (4) =pi; 

44 phchodd=phchevn+pi ; 

45 


cckl ldemod.m 

1 % CCK - 11 Mbps near optimal demodulation 

2 

3 73 y Dr. Deva K. Borah, NMSU 

4 § 4Please do not distribute without the permission of the author 

5 

6 function [d_rcv] =cckl ldemod (rcv3 , nBittx , codevec , bit valvec , phchodd , phchevn , phvec , phO) ; 

7 

8 

9 lblocks=nBittx/8; 

1 0 pr e vf ir st ph=0 ; 

11 rcv3 = rcv3(l:nBittx); 

12 rcv3=rcv3 . ♦exp (- j *phO) ; 

13 codevec = round(codevec) ; 
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14 for ii=l:lblocks 

15 testvec=rcv3((ii-l)*8+l :ii*8) ; 

16 testmat=repmat (testvec , 256 , 1) ; 

17 testmatstatus=conj (testmat) . *codevec ; 

18 [maxval , maxind] =max( sum (real (testmatstatus) . * ) ) ; 

19 detbits((ii-l)*8+3:ii*8)=bitvalvec(maxind,3:8) ; 

20 f ir stphtx=phvec (maxind , 1 ) ; 

2 1 phchange-f irstpht x-pr e vf ir s tph ; 

22 if (mod(ii,2)=0) 

23 [minval ,minind] =min (abs (exp ( j *phchodd) -exp ( j *phchange) ) . ~2) ; 

24 else 

25 [minval, minind]=min(abs(exp(j*phchevn)-exp(j*phchange)) . ~2) ; 

26 end; ‘/.if 

27 

28 secbit=mod(minind-l , 2) ; 

29 f irstbit=floor((minind-l)/2) ; 

30 detbits ((ii-1) *8+1: (ii-l)*8+2)=[f irstbit,secbit] ; 

31 pre vf ir stph=f ir stphtx ; 

32 

33 end; ‘/.for ii 

34 d_r c v=de t bi t s ; 


cck55demodsupport.m 

1 %By Dr. Deva K. Borah, NMSU 

2 ‘/.Modified by Anirudh Daga, NMSU 

3 ‘/.Please do not distribute without the permission of the author 

4 

5 function [codevec , bitvalvec ,phchodd , phchevn, phvec] =cck55demodsupport (void) ; 

6 

7 totalbitvecs=2~4; 

8 bitvalvec = de2bi([0:totalbitvecs-l]) ; 

9 for hypono=l : totalbitvecs 

10 ‘/.Generate QPSK symbols 

11 bit val=2*bitvalvec (hypono, l)+bitvalvec (hypono, 2) ; 

12 switch bitval 

13 case 0, phvec (hypono, 1)=0; 

14 case 1, phvec (hypono, l)=pi/2; 

15 case 2, phvec(hypono, l)=3*pi/2; 

16 case 3, phve c (hypono, 1) -pi; 

17 end; ‘/.switch 

18 phvec (hypono, 2) = pi*bitvalvec(hypono,3) + pi/2; 

19 phvec (hypono, 3) = 0; 
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20 phvec (hypono , 4) = pi*bitvalvec (hypono, 4) ; 

21 codevec (hypono ,1:8) = [exp( j * (phvec (hypono , 1) +phvec (hypono , 2) +phvec (hypono , 3) + . . . 

22 phvec (hypono , 4) ) ) , exp ( j * (phvec (hypono , 1 ) +phvec (hypono , 3) + . . . 

23 phvec (hypono ,4))) , exp ( j * (phvec (hypono , 1 ) +phvec (hypono , 2) + . . . 

24 phvec (hypono, 4))) , -l*exp ( j * (phvec (hypono , 1) +phvec (hypono , 4) ) ) , ... 

25 exp ( j * (phvec (hypono , 1 ) +phvec (hypono , 2) +phvec (hypono , 3) ) ) , ... 

26 exp( j* (phvec (hypono, l)+phvec (hypono, 3))) , -l*exp(j*(phvec (hypono, 1)+ 

27 phvec (hypono, 2))) , exp(j*phvec(hypono,l))3 ; 

28 end; 'Li or bitno 

29 

30 /(Define even and odd phase changes 

31 phchevn ( 1 ) =0 ; 

32 phchevn (2 )=pi/2; 

33 phchevn (3)=3*pi/2; 

34 phchevn (4)=pi ; 

35 phchodd=phchevn+pi ; 

36 


cck55demod.m 

1 7,CCK - 5.5 Mbps near optimal demodulation 

2 

3 73y Dr. Deva K. Borah, NMSU 

4 ^Modified by Anirudh Daga, NMSU 

5 ’/.P lease do not distribute without the permission of the author 

6 

7 function [d_rcv]=cck55demod(rcv3,nBittx, codevec, bitvalvec.phchodd, phchevn, phvec, phO) ; 

8 
9 

10 lblocks=nBittx/4 ; 

1 1 prevf ir stph=0 ; 

12 

13 rcv3 = rcv3(l:2*nBittx); 

14 rcv3=rcv3 . *exp(-j*ph0) ; 

15 for ii*l:lblocks 

16 testvec=rcv3( (ii-l)*8+l : ii*8) ; 

17 testmat=repmat (testvec ,16,1); 

18 testmatstatus=conj (testmat) . *codevec ; 

19 £maxval,maxind]=max(sum(real(testmatstatus) . ’ )) ; 

20 detbits((ii-l)*4+3:ii*4)=bitvalvec(maxind,3:4) ; 

21 f irstphtx=phvec (maxind , 1 ) ; 

22 phchange=f irstphtx-prevf irstph; 

23 if (mod(ii,2)==0) 
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24 [minval ,minind] =min (abs (exp ( j *phchodd) -exp ( j *ph change) ) . ~2) ; 

25 else 

26 Dttinval,minind3=min(abs(exp(j*phchevn)-exp( j*phchange) ) . “2) ; 

27 end; ‘/.if 

28 

29 secbit=mod (minind-1 , 2) ; 

30 f irstbit=f loor ( (minind-1) /2) ; 

31 detbits ((ii-1) *4+1: (ii-1) *4+2)= [f irstbit . secbit] ; 

32 prevfirstph=firstphtx; 

33 

34 end; %for ii 

35 d_rcv=detbits ; 


CCKllDemodRake_NoFor.m 

1 */.CCK - 11 Mbps near optimal demodulation (with RAKE) 

2 

3 %By Dr. Deva K. Borah, NMSU 

4 '/.Modified by Anirudh Daga.NMSU 

5 ‘/.Please do not distribute without the permission of the author 

6 

7 function [d_rcv] =CCK1 lDemodRake_NoFor (rcv3 , nBittx , codevec , bitvalvec , . . . 

8 phchodd,phchevn,phvec,phO,PSDULength,hhat) ; 

9 
10 

11 lblocks=nBittx/8; 

12 prevf irstph=0; 

13 

14 rev = rcv3(l:PSDULength+length(hhat)-l); 

15 temp = repmat((l: length (hhat) )’,l,PSDULength) + repmat(0:PSDULength-l, length (hhat) ,1) ; 

16 rcv2 = conj (repmat (hhat ( : ) , 1 , PSDULength) ) . *rcv (temp) ; 

17 rcv3 * sum(rcv2); 

18 rcv3 = rcv3 . *exp(-j *ph0) ; 

19 

20 codevec = round (codevec) ; 

21 for ii=l:lblocks 

22 testvec*rcv3 ( (ii-1) *8+1 : ii*8) ; 

23 testmat=repmat(testvec, 256,1) ; 

24 testmatstatus=conj (testmat) . *codevec ; 

25 [maxval , maxind] =Tnax (sum (read (testmat status ) . ’ ) ) ; 

26 detbits((ii-l)*8+3:ii*8)=bitvalvec(maxind,3:8) ; 

27 f irstphtx=phvec (maxind , 1 ) ; 

28 phchange=firstphtx-prevf irstph; 
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29 if (mod(ii,2)=0) 

30 [minval , minind] =min (abs (exp ( j *phchodd) -exp ( j *phchange) ) . ~2) ; 

31 else 

32 [minval .minind] =min (abs (exp( j *phchevn) -exp ( j *phchange) ) . ~2) ; 

33 end; */,if 

34 

35 secbit=mod (minind- 1 , 2) ; 

36 f irstbit=f loor ( (minind-1) /2) ; 

37 detbits ((ii-1) *8+1: (ii-l)*8+2)=[f irstbit,secbit] ; 

38 pre vf irstph=f irstphtx ; 

39 

40 end; 7,for ii 

41 d_rcv=detbits ; 


CCK55DemodeRake_NoFor.m 

1 '/.CCK - 5.5 Mbps near optimal demodulation (with RAKE) 

2 

3 %By Dr. Deva K. Borah, NMSU 

4 ‘/(Modified by Anirudh Daga, NMSU 

5 ‘/.Please do not distribute without the permission of the author 

6 

7 function [d_rcvj =CCK55DemodRake_NoFor (rcv3 ,nBittx , codevec , bitvalvec , . . . 

8 phchodd , phchevn , phvec , phO , PSDULength , hhat ) ; 

9 

10 

11 lblocks=nBittx/4; 

12 prevf irstph=0 ; 

13 

14 rev = rcv3(l:PSDULength*2+length(hhat)-l) ; 

15 

16 temp = repmat((l: length (hhat)) » ,l,PSDULength*2) . . . 

17 + repanat (0 : PSDULength*2 - 1 .length (hhat ) ,1) ; 

18 rcv2 = conj(repmat(hhat(:),l,PSDULength*2)).*rcv(temp); 

19 rcv3 ■ sum(rcv2); 

20 rcv3 = rcv3.*exp(-j*ph0) ; 

21 

22 for ii=l:lblocks 

23 testvec=rcv3((ii-l)*8+l:ii*8) ; 

24 testmat=repmat(testvec,16,l); 

25 testmatstatus=conj (testmat) . *codevec ; 

26 [marval .maxind] =max (sum (real (testmatstatus) . ’ ) ) ; 

27 detbits((ii-l)*4+3:ii*4)=bitvalvec(maxind,3:4) ; 
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28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 


f irstphtx=phvec (maxind , 1 ) ; 
phchange=f irstphtx-prevf irstph ; 
if (mod(ii,2)==0) 

[minval , minind] =min ( abs ( exp ( j *phchodd) -exp ( j *phchange ) ) . ~2) ; 

else 

[minval, minind] ==min(abs(exp(j*phchevn)-exp(j*phchange)) .‘2) ; 
end; '/.if 

secbit=mod (minind- 1 , 2) ; 
firstbit=floor((minind-l)/2) ; 

detbits ((ii-1) *4+1: (ii-l)*4+2)=[firstbit,secbit] ; 
pre vf irstph=f irstphtx ; 

end; %for ii 
d_rcv=detbits ; 


rounddec.m 

1 function y = rounddec(x, n) 

2 '/.RQUNDDEC Round to a specified number of decimals. 

3 7. 

4 % Y = R0UNDDEC(X, N) rounds the elements of X to N decimals. 

5 % 

6 % For instance, rounddec(10*sqrt(2) + i*pi/10, 4) returns 

7 % 14.1421 + 0.3142i 

8 7. 

9 % See also: ROUND, FIX, FLOOR, CEIL, ROUNDDIG, TRUNCDEC , TRUNCDIG . 

10 

11 Author: Peter J. Acklam 

12 % Time-stamp: 2004-09-22 20:06:46 +0200 

13 % E-mail: pjacklamfionline.no 

14 % URL: http://home.online.no/~pjacklam 

15 

16 '/, Check number of input arguments. 

17 error(nargchk(2, 2, nargin)); 

18 

19 7, Quick exit if either argument is empty. 

20 if isempty(x) I | isempty(n) 

21 y = []; 

22 return 

23 end 

24 

25 7, Get size of input arguments . 
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i 


26 

size_x = size(x); 


27 

size_n = size(n); 


28 

scalar _x - all(size_x == 1); % True if x is a 

scalar . 

29 

scalar_n = all (size_n == 1); % True if n is a 

scalar. 

30 



31 

% Check size of input arguments. 


32 

if ~scalar_x kk '’scalar _n kk "isequal (size_x , size_n) 


33 

error ([’When both arguments are non-scalars they must have’ 

34 

’ the same size’]); 


35 

end 


36 



37 

f = 10. ‘n; 


38 

y = round (x .* f) ./ f; 
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List of Symbols, Abbreviations, and Acronyms 

£[*] - Expectation operator 

h - Antenna height 

P(n ) - received power at index n 

fj-s - microseconds (10 -6 seconds) 

a T - rms delay 

T*; - Time delay 

ft - Mean time delay 

AP - Access Point 

BPSK - Binary Phase Shift Keying 

BER - Bit Error Rate 

BVS - Berkeley Varitronics Systems 

CCK - Complementary Code Keying 

CRC - Cyclic Redundancy Code 

dBi - deciBels referenced to an isotropic (unity gain) antenna 

dBm - deciBels referenced to 1 milliwatt 

DBPSK - Differential Binary Phase Shift Keying 

DEM - Digital Elevation Model 

DQPSK - Differential Quadrature Phase Shift Keying 

DSSS - Direct Sequence Spread Spectrum 

EIRP - Effective Isotropic Radiated Power 

GHz - GigaHertz (10 9 Hz) 

GI - Guard Interval 

GPS - Global Positioning System 

IEEE - Institute of Electrical and Electronic Engineers 

IP - Internet Protocol 

ITM - Irregular Terrain Model 

LAN - Local Area Network 

m - meter 

MAC - Media Access Control 

MGS - Mars Global Surveyer 

mW - milliWatts (10 -3 Watts) 

ms - milliseconds (10 -3 seconds) 
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NMSU 

- 

New Mexico State University 

ns 

- 

nanoseconds (10~ 9 seconds) 

OFDM 

- 

Orthogonal Frequency Division Multiplexing 

PBCC 

- 

Packet Binary Convolutional Code 

PHY 

- 

Physical Layer 

PDA 

- 

Personal Digital Assistant 

PER 

- 

Packet Error Rate 

PDP 


Power Delay Profile 

PLCP 

- 

Physical Layer Convergence Procedure 

PPDU 

- 

Physical layer Protocol Data Unit 

PSDU 

- 

Physical Sublayer service Data Unit 

QAM 

- 

Quadrature Amplitude Modulation 

QoS 

- 

Quality of Service 

QPSK 


Quadrature Phase Shift Keying 

RF 

- 

Radio Frequency 

RSSI 

- 

Received Signal Strength Indicator 

RTT 

- 

Round Trip Time 

RX 

- 

Receiver 

SFD 

- 

Start Frame Delimiter 

TCP 

- 

Transmission Control Protocol 

TX 

- 

Transmitter 

UDP 

- 

User Datagram Protocol 

W 

- 

Watt 

WEP 

~ 

Wired Equivalent Privacy 

WLAN 

- 

Wireless Local Area Network 
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